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SUMMARY 


The  use  of  mathematical  models  of  the  human  pilot  in  analyses  of  the  pilot/vehicle  system  has  b’/ought 
a  new  dimension  to  the  engineering  treatment  of  flying  qualities,  stability  and  control,  pilot /vehicle 
integration,  and  display  system  considerations.  As  an  introduction  to  such  models,  elementary  concepts 
and  specific  physical  examples  are  used  to  set  the  stage  for  a  step-by-step  development  of  what  is  known 
about  the  human  pilot  as  a  dynamic  control  component.  In  the  process,  quasi-linear  models  for  single-loop 
systems  with  visual  stimuli  and  multiloop  systems  with  visual  stimuli  are  presented  and  then  extended  to 
cover  multiloop,  multi-modality  situations.  Empirical  connections  between  the  pilot  dynamics  and  pile* 
ratings  are  also  considered. 

Some  of  the  most  important  nonlinear  features  of  human  pilot  behavior  in  adapting  to  changes  in  the 
character  of  the  stimuli  are  described  and  tied  to  the  quasi-linear  models  via  the  Successive  Organization 
of  Pore  notion  (SOP)  theory,  which  is  reviewed  and  elaborated.  Dual-mode  control  models  needed  to  describe 
the  p  „’s  behavior  in  response  to  sudden  transients  are  present  :d,  along  with  pursuit  and  compensatory 
elements  of  the  SOP  continuum. 

The  current  status  of  mathematical  pilot  models  is  shown  to  cover  random,  random-appearing,  and 
transient  inputs  for  single-  and  multi-loop  system  configurations.  An  extensive  bibliography  of  applica¬ 
tions  and  a  nummary  of  analysis  problems  which  have  been  addressed  is  included,  as  is  a  short  general 
status  summary  and  critique  of  existing  models  jn  the  fora  of  a  listing  of  shortcomings  and  problem  areas. 
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CHAPTER  I 


INTRODUCTION 


A.  PURPOSE  AND  OUTLINE!  OF  THE  REPORT 

From  the  authors’  standpoint,  the  genesis  of  this  report  was  a  request  from  the  AGARD  Aerospace  Medical 
Panel  to  prepare  a  concise  and  critical  review  of  mathematical  modeling  of  the  human  pilot  as  a  control 
element  (as  distinct  from  the  somewhat  parallel  technology  in  passive  biodynamic  modeling).  It  was  hoped 
by  the  "mel  that  this  would  highlight  the  most  relevant  and  reliable  conclusions  about  human  operator 
mathematical  modeling  to  be  drawn  from  the  available  literature. 

Our  approach  to  this  tack  is  constructed  upon  two  bases.  The  first  is  an  indication  of  what  is  known 
about  the  human  pilot  as  a  dynamic  control  component}  and  the  second  is,  in  the  light  of  this  knowledge, 

'•bat  further  is  needed,  or  where  do  we  go  from  here?  There  is  an  adequate  variety  of  substantially  current 
surveys  of  the  literature  and  models  readily  available  to  the  specialist.  Therefore,  it  is  not  our  inten¬ 
tion,  to  contribute  another  survey.  We  will,  instead,  attempt  to  distill  out  concepts  and  models  which 
represent  the  essence  of  current  generalised  empirical  know, edge.  These  are  inherently  syntheses  of  models 
and  empirical  data  based  on  the  efforts  of  many  investigators.  Our  emphasis  will  be  on  quasi-linear 
descriptions  for  the  most  part,  since  these  are  the  most  complete  and  have  been  most  extensively  applied. 

On  the  current  frontier  are  nonlinear  and  non-stationary  situations,  so  these  will  also  be  treated  to  the 
extent  possible.  The  next  sections  of  this  introduction  give  brief  backgrounds  of  the  nistory,  intended 
nature,  and  fundamental  concepts  of  human  pilot  mathematical  descriptions.  These  remarks  set  the  stage 
for  the  remainder  of  this  report. 

Chapter  IX  presents  some  elementary  concepts  for  the  quantitative  description  of  triman  dynamic  response 
in  closed-loop  systems.  The  intent  in  this  chapter  is  to  introduce  the  basic  concepts  and  phenomena  in  a 
Primarily  descriptive  and  intuitively  understandable  manner.  By  means  of  concrete  examples,  we  develop 
mathematical  descriptions  (integro-difl’erencial  equations  and  transfer  functions),  the  concept  of  the  rem¬ 
nant,  and  other  facets  of  human  operator  mathematical  descriptions.  The  physical  approach  of  this  presenta¬ 
tion  is  emphasized  by  reliance  on  time  hist  ri.es  to  establish  the  concepts  and  approximate  mathematical 
relationships.  The  lata  for  the  several  specific  cases  presented  are  then  synthesi-ed  into  a  simple  overall 
model  which  characterises  the  dynamics  of  the  human  and  control  system.  The  mathematical  implications  of 
this  model  are  then  explored  in  detail  to  indicate  the  types  of  estimates  and  predictions  one  could  make 
about  control  systems  which  behave  in  the  fashion  exhibited  by  the  simple  model.  Sufficie- 1  detail  is 
provided  to  permit  the  simple  "back  of  the  envelope"  calculations  so  dear  to  the  heart  of  practical 
scientists  and  engineers. 

The  descriptions  of  Chapter  II  are  for  one  pilot  subject  <»nd  the  model  building  there  depends  on  *emi- 
qualitative  measurements.  Chapter  III  describes  the  more  general  state  of  affairs  which  obtains  when  many 
subjects  are  present  and  when  a  more  respectable  mathematical  theory  is  used  as  the  basis  for  measurement. 
Thus,  the  four  key  types  of  variable  which  affect  the  pilot's  action  are  described  in  detail,  and  the 
nature  of  quasi-linear  measurements  is  introduced. 

Chapters  IV  and  V  summarize  the  current  status  of  juasi-iinear  models  for  compensatory  situations.  We 
start  in  Chapter  IV  with  single-loop  systems  for  visual  stimuli.  The  preser nation  includes  an  indication 
of  the  experimental  data  available  and  interpretation  of  these  in  terms  of  the  most  elaborate  current 
models.  Chapter  V  follows  with  the  extension  of  the  single-loop  quasi-linear  compensatory  model  to  mu..-i- 
loop  visual  input  and  to  multiloop,  multi-modality  situations.  Some  empirical  connections  between  pilot 
dynamics  and  pilot  ratings  are  considered  in  ooth  Chapters  IV  and  V. 

Chapter  VI  examines  the  endpoints  of  key  features  of  human  operator  nonlinear  behavior,  that  is,  the 
adaptive  changes  in  behavior  induced  by  changes  in  the  character  of  the  stimulus  presentee.  Pursuit  and 
"precognitive"  behavior  is  introduced  and  tied  together  with  compensatory  behavior  via  the  Successive 
Organization  of  Perception  (SOP)  hypothesis,  which  is  reviewed  and  elaboi*ated.  Step  inputs  and  nearly 
periodic  forcing  functions  are  among  'hose  considered,  and  the  dual-node  control  model  needed  to  describe 
the  pilot's  behavior  in  response  to  sudden  transients  is  featured  as  one  element  in  the  SOP  continuum. 

The  status  quo  summarized  by  Chapters  II- VI  provides  the  prologue  to  Chapters  VII  and  VIII.  Although 
this  report  is  not  intended  1c  provide  a  survey  of  applications  techniques  and  examples,  part  of  an  ade¬ 
quate  status  report  is  an  indication  of  application^.  Therefore,  Chapter  VII  briefly  summarizes,  with  the 
aid  of  an  extensive  bibliography,  the  types  of  pilot/vehicle  analysis  problems  which  have  been  addressed. 

The  critique  portion  of  the  report  follows  in  Chapter  VIII.  This  takes  the  form  of  a  listing  of 
deficiencies  in  the  existing  models  which  thereby  define  shortcomings  and  problem  areas.  Some  attempt  is 
made  to  make  this  critique  constructive  by  noting,  where  available,  analytic  or  other  approaches  which 
might  be  adopted  to  alleviate  the  shortcomings. 

B.  HISTORICAL  BACKGROUND 


Piloted  aircraft  have  always  required  a  satisfactory  match  of  the  aircraft  characteristics  and  the 
controller  properties  of  the  human  pilot.  This  was  explicitly  recognized  in  the  December  2J,  i^Of,  Signal 
Corps  Specification,  No.  1*86,  for  a  heavier-than-air  flyir.g  machine,  as  "....it  must  be  steered  in  all 
directions  without  difficulty  and  at  all  times  under  perfect  control  and  equilibrium."  (Ref.  1).  But, 
need  for  good  man-machine  integration  was  thoroughly  appreciated  even  earlier  by  the  Wright  Brothers,  for 
instance,  in  a  talk  before  the  Western  Society  of  Engineers  in  1901,  two  years  before  their  first  powered 
flight,  Wilbur  Wright  said  (Ref.  2): 
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"Men  already  know  how  to  construct  wings  or  aeroplanes,  which  when  driven  through  the 
air  at  sufficient  speed,  will  not  only  sustain  the  weight  of  the  wings  themselves,  but 
also  that  of  the  engine,  and  of  the  engineer  as  well .  Men  also  know  how  to  build 
engines  and  screws  of  sufficient  lightness  and  power  to  drive  these  planes  at  sustain¬ 
ing  speed. . .  .Inability  to  balance  ar.d  steer  still  confronts  students  of  the  flying 
problem. .. .When  this  one  feature  has  been  worked  out,  the  age  of  flying  machines  will 
have  arrived,  for  all  other  difficulties  a~e  of  minor  importance." 

Although  essential  for  success,  an  agreeable  marriage  between  the  dynamic  properties  of  the  inanimate 
aircraft  and  the  animate  pilot  is  not  spontaneously  achieved  in  the  design  process.  Indeed,  the  provision 
of  proper  aircraft  flying  qualities  has  often  poseu  serious  problems  which  the  designer  must  solve.  Until 
fairly  recently  these  solutions  relied  very  heavily  on  intuitive  eut-and-try  procedures.  Over  the  years 
this  ad  hoc  empirical  approach  fostered  man}*  of  the  udventures  and  uncertainties  of  flight  testing! 

The  desire  to  handle  aircraft  stability  and  control  problems  in  a  more  analy ' cal  fashion  was  recognised 
long  ago.  As  an  illustration,  before  World  War  II  Koppen  stated  (Ref.  *): 

"Since  the  controlled  motion  of  an  airplane  is  a  combination  of  airplane  and  pilot 
characteristics,  it  is  necessary  to  know  something  about  both  airplane  and  pilot 
characteristics  before  a  satisfactory  .jol  of  airplane  design  can  be  done." 

But  the  central  difficulty  in  understanding  controlled  motion  was  recognised  earlier  still.  For  example, 

W.  Crowley  and  Sylvia  3 '•'an  remarked  in  a  19*0  Aeronautical  Research  Committee  report  (Ref.  ■'*): 

"A  mathematical  investigation  of  the  controlled  motion  is  rendered  almost  impossible 
on  account  of  the  adaptability  of  the  pilot.  Thus  if  it  is  found  that  the  pilot 
operates  the  controls  of  a  certain  machine  according  to  certain  lave,  and  so  obtains 
the  oe. t  performance,  it  cannot  be  assumed  that  the  same  pilot  would  applv  the  same 
laws  to  another  machine.  Ke  would  subconsciously,  if  not  intentionally,  change  his 
methods  to  suit  the  new  conditions,  and  the  various  laws  possible  to  a  pilot  are  too 
numerous  for  a  general  analysis." 

Actually,  matters  are  even  worse  than  Crowley  and  Skar.  recognised;  for  while  much  of  the  pilot's  dynamic 
behavior  is  governed  by  the  aircraft  dynamics,  many  additional  factors  also  affect  his  properties. 

Psychologists  too  pointed  to  variability  as  the  hallmark  of  human  behavior,  and  organised  much  of  the 
theory  useful  in  experimental  psychology  so  as  to  cope  with  intersub.jeet  and  intrasubject  differences. 

But,  the  biggest  problem  confronting  anyone  who  wished  to  study  the  control1 «d  motion  of  aircraft  — 
what  we  would  today  call  pilot  Vehicle  analysis  —  was  not  the  human's  adaptability  and  variability  but 
the  absence  of  an  underlying  quantitative  theory  on  which  to  erect  a  structured  approach  to  the  manual  con¬ 
trol  *f  aircraft.  Wha'  was  needed  was  a  theory  for  feedback  control  systems,  and  this  became  available 
in  a  Sufficiently  mature  state  in  the  early  19LOS.  Witn  feedback  control  theory  available  as  a  paradigm, 
on:,  a  need  had  to  be  at  hand  for  attempts  to  be  made  to  overcome  the  a  priori  pessimism.  And  at  the 
te -inning  of  the  Cs  there  was  no  end  to  the  needs!  Thus,  necessity  was  the  mother  of  human  dynamics 
rs‘-ear.  h.  Complex  weapons  were  evolving  which  could  only  function  in  concert  .  ith  human  operators,  so  ail 
urgent  requirement  for  engineering  data  arose. 

The  pioneer  in  human  operator  dynamic  measurements  was  Arnold  Tuotin,  the  electrics*,  engineer,  in 
England  during  World  War  II.  Tustin  extended  the  required  feedback  control  theory  framework  cy  introducing 
the  concept  of  "describing  function"  and  "remnant"  measures  and  quasi-linear  systems  in  General  ( Re f .  t. 

He  then  applied  the  concepts  to  actual  human  operations.  In  reporting  on  nis  ct-’iie"  of  manual  control  of 
a  power-driven  gun  'Refs.  •  and  '  ): 

"The  ooje-'t  of  the  series  of  tests  ..as  to  investigate  the  nature  of  the  layer's 
response  in  a  number  of  particular  cases  and  to  attempt  to  find  the  laws  of  rel1' - 
tion'hips  of  movement  to  error.  In  particular,  it  was  hoped  that  this  relation¬ 
ship  might  be  found  ( within  the  range  of  practical  req  .iro’.ents)  to  be  approximately 
linear  end  so  permit  the  well  developed  theory  of  'linear  servomechanisms'  to  be 
applied  to  ;  anuai  control  in  tne  same  way  as  it  3  a;,  ieu  to  automatic  following." 

Alto  luring  the  3- .ond  World  War  and  responding  to  the  same  urgent  needs  but  acting  independently  of 
Tv-  .  .  A.  Cobciyk  ana  R.  3.  Phillips  at  the  HIT  Radiation  Laboratory  (Ref.  r'l  and  H.  K.  Weiss  (Refs.  9 
and  10)  a'  Aberdeen  Proving  Ground  presumed  quasi-lino'**’  operation  in  a  series  of  studies  on  aided  trussing 
of  puns. 

■  fter  the  war  t'.ese  seminal  efforts,  and  the  nope  for  a  more  rational  approach  to  the  design  of  aircraft, 
del  Leo  Chattier  of  the  U31I  Bureau  of  Aeronautics  and  Charles  Westbrock  of  the  U3  air  Force  Aircraft  Labora¬ 
tory  i.  sponsor  sore  small-scale  ’ontraet  research  efforts  aimed  a.  determining  the  dynamic  characteristics 
of  human  pilots,  lor  the  Ha** y,  lie  Goodyear  Aircraft  group  of  R.  Meade,  I.,  biarantidcs,  and  A.  Cacioppu, 
headed  i>y  R.  "ayn»,  developed  excellent  analog  comp  iter  representations  for  pilots  for  two  specif  e  task- 
variable  situations  (Refs.  11-1-1.  The  Air-Force-sponsored  activity  at  cranuiin  Institute  (Ref.  1  ), 
Princeton  University  (Refs.  1*-1°),  and  Control  Specialists,  dnc.  (Refs.  !•  and  lp'  chose  to  exploit  cross¬ 
correlation  and  cross-spectral  techniques  to  establish  human  pilot  dynamics .  Parallel  university  research 
-as  mderway  at  MIT  with  a  remarkable  Master's  thesis  by  Lindsey  Russell  (Ref.  fO)  and  an  extensive  and 
useful  doctoral  dissertation  by  Jerome  Elkind  (Ref.  l). 

The  end  of  the  pioneering  era  in  experimentation  can  be  •onveniently  put  with  the  publication  of 
"Dynamic  Response  of  Human  Operators"  (°of.  19)  in  October,  19*.  and  its  Journal  summary  in  Ref.  7 '  .  This 
volume  codified  and  ’orrelated  the  available  lmnan  response  data,  developed  predictive  models  compatible 
with  these  data,  and  prescribed  preferred  forms  for  he  operator  which  permitted  specification  01  "ideal" 
characteristics  for  the  controlled  element  compatible  *ith  the  man  in  the  manual  control  system. 


tSL-t .  S-a-rw* 


Since  19?7  an  enormous  number  of  measurements  have  been  made  for  single-loop,  i.e.,  single-input/ 
single-output  manual  control  systems  of  all  kinds.  Using  mere  and  more  refined  measurement  and  data  reduc¬ 
tion  techniques,  many  organisations  and  individuals  have  contributed  in  the  U  S.  and  abroad.  Many  of  these 
will  be  mentioned  later  when  we  attempt  to  define  where  we  stand  today.  In  terms  of  the  historical  record, 
it  is  probably  fair  to  say  that  now  much  more  is  known  about  the  input-output  characteristics  of  some  kinds 
of  single -loop  systems  than  remains  to  be  learned.  The  single-loop  experience  also  provided  a  firm  founda¬ 
tion  for  subsequent  extensions. 

Going  from  single-loop  to  multiloop  situations  is  an  enormous  step  —  so  large,  in  fact,  that  it  is 
still  in  progress.  The  first  attempts  wore  multi-modality  experiments  in  which  a  number  of  fundamental 
dynamic  response  measurements  were  made  in  aircraft  and  moving-base  simulators  to  determine  the  effects  of 
linear  and  rotary  motion  cues  on  the  pilot's  dynamics.  Ordinarily,  these  studies  compared  fixeu-  versus 
moving-base  situations  on  the  basis  of  effective  visual  input  pilot  describing  functions.  In  two  or  three 
instances  it  has  been  possible  to  separate  uniquely  the  motion  and  visual  transfer  characteristics  by  using 
independent  forcing  functions.  As  will  be  seen  later,  all  this  amounts  to  an  excellent  start,  although 
much  of  importance  is  on  the  docket  for  the  future. 

As  one  might  expect  from  this  brief  historical  account,  three  decades  of  joint  engineering  and  experi¬ 
mental  psychological  attempts  to  characterise  mathematically  the  dynamic  characteristics  of  human  pilots 
has  resulted  in  an  enormous  collection  of  literature.  Because  most  of  the  work,  in  che  United  States  at 
least,  has  been  sponsored  by  the  Government,  notably  the  Air  Force  Flight  Dynamics  Laboratory,  the  National 
Aeronautics  and  Space  Administration,  and  the  Naval  Air  Systems  Command,  or  has  been  the  result  of  univer¬ 
sity  research  activities,  the  complete  documentation  is  predominantly  in  Government  reports  and  university 
theses.  Stemming  from  this  enormous  base,  like  the  visible  tip  of  an  iceberg,  is  a  journal  literature 
which  is  also  extensive.  Much  of  the  original  work  is  reflected  in  the  IEEE  Transactions  (Human  Factors  in 
Electronics;  Man-Machine  Systems;  Systems,  Man,  and  Cybernetics),  the  AIAA  Journal  of  Aircraft,  the  Journal 
of  the  Human  Factors  Society,  and  Ergonoml/s ■  Particularly  complete  documentation  in  recent  years  has  come 
about  from  the  "Annual  Manual,"  officially  the  Annua?  Conference  on  Manual  Control,  sponsored  by  NASA  and 
the  USAF  (Refs.  These  international  annual  meetings  present  summaries  of  the  latest  word  in  manual 

control  technology  as  well  as  a  useful,  but  incomplete,  cross  section  of  applications.  Extensive  biblio¬ 
graphies  are  given  in  Refs.  This  literature  reflects  the  application  of  almost  every  conceivable 

control  theory  approach  to  modeling  the  human  operator,  as  well  as  the  application  o.  appropriate  models 
to  the  control  of  almost  any  imaginable  vehicle- 

As  befits  the  vapid  evolution  of  mathematical  models  for  the  human  operator  through  the  past  three 
decades,  a  series  of  summary  surveys  have  appeared  aperiodically  (Refs.  19,  32-38).  Thus,  not  only  the 
dynamic  characteristics  and  models  ®f  the  human  operator  but  some  synthesis,  as  depicted  in  summary  surveys, 
are  well  represented  in  a  readily  available  literature. 

C .  MODELING  PILOTED  SYSTEMS 

Engineering  models,  as  we  use  the  terms,  comprise  analytic  and  verbal  descriptions  of  the  phenomenon 
of  interest  ?o  expressed  that  the  key  variables  are  explicit;  the  ranges  and  rules  of  application  ar-»  well 
defined;  and  all  of  the  relevant  constructs  described  in  operational  terms  compatible  with  such  other 
models  as  may  be  pertinent .  When  this  has  been  done  ouch  models  are  (and  to  be  at  all  convincing  must 
be)  capable  ®f  being  validated,  «>gmented,  or  disproved  by  researchers  other  than  the  originators.  These 
models  are  will,  more  than  convenient  and  efficient  summaries  of  data;  they  are  predictive  tools  for 
rational  systems  synthesis.  In  fact,  the  most  compelling  justification  of  any  model  is  its  capacity  to 
subsume  past  experimental  results  and  to  predict' the  outcome  of  future  experiments  especially  contrived  to 
exercis®  the  motel  beyond  its  original  scope.  Once  these  engineering  models  have  been  validated,  they 
serve  several  purpo’oe.  rev  instance,  they: 

1.  Provide  "laws"  which  can  be  used  t>.  estimate  or  predict  pilot  dynamic  behavior. 

.  Characterise  past  observations  with  simple  analytical  forms,  thereby  achieving 
economies  in  description  of  experimental  data. 

7  •  'Jiant  itativelv  describe  the  connections  between  model  parameters  and  situation 
variables  ..hi  -h  affect  the  pilot. 

.  Establish  the  basis  for  a  rationale  .diich  -an  be  used  to  better  "understand"  the 
ways  in  which  the  pilot  behaves  as  a  control  system  component. 

lie'- a cse  the  purpose  of  engineering  models  is  application,  a  host  of  variables  which  appear  to  be 
interesting  intuitively  can  easily  tarn  out  to  te  irrelevant  to  the  particular  application  and,  thereby, 
t  the  m  del.  However,  -bother  they  be  relevant  or  irrelevant  is  not  an  intuitive  judgment.  The  decision 
tv  It fro'  a  careful,  alcost  tedious,  statement  of  potentially  relevant  'variables  and  the  establishment 
of  their  influence  ®it,i®r  v  pirieally  or  theoretically,  so  that  these  can  be  made  analytically  and/or 
verbally  explicit  in  the  model. 

Vf®  have  denoted  the  relevant  variables  as  key  variables  and  have  divided  them  into  four  categories: 
tack,  environmental,  pilot- centered,  and  procedural.  They  will  be  described  in  detail  further  along,  but 
the  point  for  now  is  the  need  to  recognise  and  take  into  explicit  account  the  principal  factors  which 
underlie  h-ran  operator  variability  —  thus  educing  to  a  minimum  the  inevitable  component  of  variability 
v,hi®h  remains  urns  signed.  The  engineering  model  then  becomes  more  generalizable  in  that  its  parameters 
are  affected  in  an  empirically  discoverable  way  by  the  identified  key  variables.  Further,  our  understand¬ 
ing  of  vaguely  defined  .  ey  variable  elements,  e.g.,  "motivation"  or  "fatigue"  as  typical  pilot- -entered 
'variables,  is  fjrther  enhanced  by  associating  quantitative  changes  in  measurable  pilot  parameters  with 
'hangeo  introduced  in  these  variables.  In  summary,  an  effective  model  of  human  dynamic  operation  must 
recogni  e  that  inter-  and  intra-subject  variability  exists  and  that  much  of  this  variability  can  be 
assigned  as  the  effects  of  key  variables  on  the  model’s  parameters. 
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D.  PILOT  ZKKAMIC  BEHAVIOR  IK  COKTRDL  TERMS 

A  pilot's  actions  are  many  and  their  range  is  broad.  We  are  .interested  in  those  actions  which  exert 
control  on  the  aircraft,  and  we  will  outline  the  model  structure  implications  of  these  activities.  Subse¬ 
quently,  we  will  indicate  the  variety  of  factors  which  can  influence  the  resulting  performance;  for  by 
limiting  the  conceptual  context  and  defining  these  factors  it  becomes  possible  to  convert  an  unmanageably 
complicated  general  problem  to  a  manageably  complicated  set  of  specific  problems. 

The  conceptual  context  derives  from  the  following  three  fundamental  concepts  which  constitute  the 
essential  substance  of  pilot-in-the-loop  systems  analysis: 

•  The  establishment  of  appropriate  control  loops  among  the  display,  pilot,  and  vehicle 
system  elements. 

•  The  optimization  of  the  loop  properties  by  suitable  pilot  adaptation. 

•  The  penalties  for  this  adaptive  control. 

The  initial  concept  is  that  to  accomplish  guidance  and  control  functions  (as  in  flying  approaches,  inter¬ 
cepts,  formation,  regulating  against  atmosphere  disturbances,  etc.)  the  human  pilot  est._  fishes  a  variety 
of  closed  loops  around  the  airplane  which,  by  itself,  could  not  accomplish  such  tasks.  The  greatest  amount 
of  information  and  experience  exists  for  the  single-loop  case  for  which  the  pilot  is  applying  his  full 
attention  to  one  control  task.  Because  many  piloting  problems  involve  one  axis  dominantly,  and  because  a 
wide  variety  of  displays,  configurations,  vehicle  dynamics,  and  manual  control  mechanisms  can  be  reduced 
to  i  single-loop  block  diagram,  concentration  on  this  simplified  structure  has  been  highly  productive . 

When  multiple  inputs  and  irredueiole  multiple  loops  exist,  more  complicated  analysis  cannot  be  avoided.  It 
is  often  the  case  that  the  input  data  upon  which  the  pilot  depends  and  the  loops  he  elects  to  close  are  not 
immediately  obvious  and  must  be  deduced  from  his  measured  dynamics  and  performance.  After  recognizing  such 
data  sources  and  loops  it  is  possible  for  the  system  engineer  to  enhance  their  utility  for  the  pilot  by 
appropriate  design  refinements .  In  this  fashion  the  peculiarly  human  contrioutions  to  the  system  may  be 
strengthened  and  the  pilot  will  be  treated  with  more  design  "respect." 

The  second  major  concept  is  that  to  be  satisfactory  these  various  closed-loop  systems  have  to  behave 
in  a  ^iiituoie  fauuori,  although  animate  and  inanimate  components  are  interacting,  the  overall  system  must 
share  certain  of  the  qualitative  features  of  all  "good"  closed-loop  control  systems  (stability,  accuracy, 
simplicity,  etc.).  As  the  adaptive  means  to  accomplish  this  end,  the  pilot  must  make  up  for  auj  dynamic 
deficiencies  of  the  effective  controlled  element  a.  a  whole  by  appropriate  adjustments  of  his  own  dynamic 
properties.  The  view  of  the  pilot  is  that  of  an  active,  sophisticated,  and  knowledgeable  system  component, 
not  an  afterthought  which,  in  some  quasi-mechanisti:  manner,  must  fit  into  the  control  system.  The  pilot 
expresses  his  operational  judgments  on  criteria  for  system  performance  and  on  satisfactory  system  stabili¬ 
zation  by  his  actions.  In  every  sense  he  is  in  control  and  serves  to  enrich  the  system's  task  capability 
repertory. 

The  third  major  concept  states  that  there  is  a  cost  to  the  pilot  for  this  adjustment:  in  workload- 
induced  stress;  in  concentration  of  his  faculties;  and  in  a  reduced  potential  for  coping  with  the  unexpected 
The  measures  of  the  cost  to  the  pilot  are  pilot  commentary  and  pilot  rating,  possible  task-associated  change 
in  the  pilot's  dynamics,  and  other  workload- sensitive  physical  and  psychophysiological  measures  such  as 
heart  rate.  The  cost  to  the  system  is  expressible  terms  of  diminished  operational  reliability  and  effec¬ 
tiveness.  fince  these  costs  can  be  high  indeed,  efforts  to  predict  their  onset  and  magnitude  by  measures 
based  on  the  pilot's  behavior  are  of  great  potential  benefit. 


CHAPTER  II 


ELEMENTARY  CONCEPTS  FOR  THE  QUANTITATIVE  DESCRIPTION  OF 
HUMAN  OPERATOR  DYNAMICS  IN  CONTROL  SYSTEMS 


When  n  human  is  successfully  engaged  as  an  active  element  in  a  feedback  control  system  the  combined 
behavior  of  man  and  machine  is  such  as  to  satisfy  overall  man/machine  system  purposes.  Machines  and  their 
manipulative  controls  differ;  the  external  environment  and  surround  from  which  the  operator  derives  sensory 
input  may  be  highly  diverse;  and  control  purposes  may  bear  little  resemblance  to  one  another.  Yet  the  human 
controller  counters  all  this  diversity  by  modifying  his  characteristics  appropriately  to  match  the  many 
possible  control  situations.  The  fact  of  these  adjustments  and  adaptation  is  characteristic  of  the  hunan; 
the  consequence  is  a  bewildering  range  of  behavior  which  makes  quantitative  description  of  the  human  opera¬ 
tor  enormously  complex  when  viewed  in  the  large.  Because  the  entirety  is  incomprehensible,  in  c"der  to 
gain  quantitative  understanding  it  is  necessary  to  partition  the  whole  into  parts.  If  we  are  lucky,  these 
parts  can  be  selected  so  as  to  have  behavioral  properties  which  efficiently  capture  the  dominant  effects 
and  are  susceptible  to  simple  quantitative  description.  In  this  way,  there  is  some  hope  of  obtaining  an 
ensemble  of  reasonably  simple  quantitative  behavioral  models,  each  of  which  has  adequate  generality  in  its 
restricted  sphere.  Then,  when  efficient  models,  each  suitable  for  one  of  several  related  facets,  are 
available,  further  inductive  generalization  may  be  feasible.  This  is  the  approach  we  will  pursue  here. 


The  situations  which  are  most  susceptible  to  quantitative  description  involve  control  of  inanimate 
elements  which  have  constant  characteristics,  and  an  operator  who  has  learned  by  extensive  practice  to  con¬ 
trol  the  system  effectively  with  little  fluctuation  in  the  man/machine  system  performance  measures.  When 
the  stimuli  impacting  on  the  system  are  random  or  random-appearing  time-stationary  processes,  the  skilled 
operator  will  ultimately  adjust  so  that  his  control  actions  (as  time  signals)  have  similar  random  and  time- 
stationary  properties,  at  least  over  all  but  a  very  short  or  very  long  period.  The  simple >t  manual  control 
system,  rig.  1,  which  fits  these  characteristics  is  the  so-called  single-loop  compensatory  system  with 

visual  stimulus,  which  is  the  paradigm 
for  many  much  more  complicated  appear¬ 
ing  configurations.  The  system  forcing 
function,  i(t),  is  a  random  or  random¬ 
appearing  time  function  which  has  sta¬ 
tionary  or  quasi-stationary  properties. 
Tlie  compensatory  display  presents  the 
operator  with  a  visual  stimulus  which 
shows  only  the  difference  between  the 
system  forcing  function  and  the  system 
output.  Typically,  the  operator's  task 
is  to  minimize  this  presented  error 
signal,  e(t),  by  trying  to  keep  it 
superimposed  on  a  stationary  point  or 
line  on  the  display.  This  is  accom¬ 
plished  by  manipulative  control  action, 
c(t),  which  affects  the  controlled  ele¬ 
ment  and  gives  rise  to  the  system  out¬ 
put,  m(t),  being  controlled.  The  usual 
purpose  o*'  a  system  of  this  nature  is 
to  make  the  system  output  closely 
resemble  the  system  forcing  function 
or,  in  other  words,  to  make  the  output 
follow  the  input.  The  quality  of  this 
following  is  indicated  by  the  system 
error  which  is  directly  proportional 

Figure  1.  Display  and  Functional  Block  Diagram  of  to  the  operator's  visual  stimulus. 
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In  order  to  obtain  some  piysioal 
appreciation  for  human  operator  behavior 
in  such  systems,  we  shall  so  tsider  a 

few  simple  cases.  Several  kinds  of  human  operator  dynamic  behavior  can  be  evoked  by  changing  the  controlled 
element  dynamics  while  retaining  a  common  control  purpose,  i.e.,  minimization  of  the  error  for  sil  the  con¬ 
trol  systems  studied.  For  each  control  situation,  a  simplified  characterization  of  the  human  •'perator  will 
be  determined.  Then  the  human  operator  models  for  the  several  specialized  control  situation*'  will  be  con¬ 
sidered  as  representative  elements  of  all  single-loop  compensatory  manual  control  systems  subjected  to  a 
stationary  ranaom-appearing  forcing  function  in  an  attempt  to  generalize.  The  generalization  achieved  will 
be  a  simplified  overall  model  of  manual  control  system  dynamics  which  connects  human  operator  and  controlled 
element  dynamic  forms  by  a  simple  equation.  This  equation  also  reveals  how  the  operator's  dynamics  are 
affected  by  controlled  element  form. 


The  simple  mathematical  models  for  manual  control  system  dynamics  developed  on  this  ad  hoc  basis  define 
a  particular  form  of  feedback  control  system.  The  mathematical  laws  for  this  type  of  system  will  then  be 
exercised  to  illustrate  their  use  in  analysis  and  prediction. 


A.  APPROXIMATE  DESCRIPTION  OF  HUMAN  OPERATOR  DYNAMICS 
FOR  ELEMENTARY  COMPENSATORY  SYSTEMS 


In  this  section  we  shall  consider  several  pieces  of  experimental  data  taken  with  a  single  highly  trained 
operator  for  a  variety  of  controlled  element  dynamics  composed  of  combinations  listed  and  described  in 
Table  1 .  To  make  matters  as  physical  and  specific  as  possible,  the  experimental  data  and  discussion  will 
emphasize  time  histories  and  time  domain  operations.  Nonetheless,  the  controlled  element  properties  given 
in  Taole  1  indicate  both  time  and  frequency  domain  characteristics  and  illustrate  how  the  element's  defining 
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TABLE  1 .  ELEMENTARY  CONTROLLED  ELEMENT  DYNAMICS 


controlled  Clement  dynamics 

TIME  DOMAIN  CHARACTERISTICS 
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parameters  are  reflected  in  each  domain.  This  serves 
as  a  con/enienn  introduction  to  the  frequency  domain 
which  will  be  used  in  subsequent  chapters. 

Although  the  controlled  element  dynamics  in 
Table  1  are  very  simple,  they  have  many  parallels  in 
practical  control  systems.  For  example,  the  rate 
(velocity)  controlled  element  is  an  idealization  for 
pitch  attitude  control  with  a  heavily  augmented  air¬ 
craft,  heading  con  rol  in  automobile  driving,  bank 
attitude  control  in  an  aircraft  which  has  a  rapid 
roll  subsidence  mode,  etc.  But  the  most  important 
feature  these  elementary  controlled  elements  share 
with  their  more  complex  cousins  which  represent  real- 
world  control  situations  more  adequately  is  that, 
with  the  exception  of  the  proportional  type,  none  of 
them  will  stay  put  by  themselves.  That  is,  they 
respond  to  an  operator  input  in  a  dynamic  fashion  so 
that  the  controlled  element  output,  m(t),  continues 
to  change  or  remains  non-zero  after  the  operator's 
output,  c(t),  has  stopped.  Thus,  to  cause  the  system 
output  to  behave  as  desired,  the  operator  must  exert 
control  activities  much  of  the  time.  In  the  compen¬ 
satory  situations  to  be  considered  here,  this  control 
action  is  inherently  of  a  feedback  system  nature 
because  the  operator's  sole  stimulus  is  the  error 
presented  on  the  display.  In  describing  the  opera¬ 
tor's  behavior  then,  the  key  question  is  how  the 
operator,  as  a  dynamic  element  in  the  closed-loop 
system,  equalizes,  scales,  and  adjusts  quantities 
derivable  from  the  error  signal  to  form  his  control 
output,  c(t).  In  other  words,  we  wish  to  investigate 
the  dynamic  functional  relationships  between  stimulus 
and  response  in  the  closed- loop  system  context. 
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Figure  2.  Man/Machine  System  Response 
Time  Histories;  Yc  =  Kc/s 


As  a  starting  point,  consider  the  signals  present  at  various  points  in  the  closed  loop  when  the  con¬ 
trolled  element  dynamics  are  Yc  =  Kc/s.  A  sample  set  of  such  signals  is  shown  in  Fig.  2.  These  were 
generated  using  as  a  subject  a  professional  instructor  pilot  who  is  also  highly  trained  in  the  controlled 
tasks  to  be  presented  here.  The  control  situation  is  compensatory  (Fig.  l)  and  involves  a  random¬ 
appearing  forcing  function  which  is  to  be  followed  by  the  controlled  element  output.  The  system  error  is 
presented  on  a  cathode  ray  tube  and  the  operator  manipulates  a  stiff  stick  grasped  between  thumb  and  fore 
finger  to  exert  control  over  the  controlled  element  dynamics  generated  in  an  analog  computing  setup. 
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The  first  thing  to  notice  about  the  typical  time  histories  is  that  the  system  output,  m,  follows  the 
system  forcing  function,  i,  very  closely.  Only  a  slight  time  lag  Keeps  the  output  from  being  a  nearly 
identical  duplicate  of  the  forcing  function.  The  actual  difference  between  these  two  quantities  is  shown 
by  the  error,  which  is  quite  small  and  fluctuates  about  zero  in  a  random-like  fashion.  It  is  in  response 
to  this  stimulus  that  the  operator  develops  his  output,  c.  Here  w>,  find  a  remarkable  thing.  The  operator 
output  tends  to  look  very  much  like  a  scaled  and  delayed  version  of  the  error.  In  fact,  if  we  simply 
advance  the  output  trace  by  a  time,  =  O.lt  sec,  the  two  signals  appear  to  be  directly  proportional,  with 
a  proportionality  constant,  Kp2,  except  for  a  few  tiny  "random  wiggles"  here  and  there  on  the  output.  If 
we  now  grace  this  observation  with  a  quantitative  description,  we  can  write  that  the  operator  output  is 
approximately  related  to  the  system  error  or  operator  stimulus  by: 

c(t  +  r?)  i  Kp^e(t)  0) 

or: 

c(t)  i  Kp?e(t  -  tj>)  (&) 

If  the  relationship  is  considered  exact  rather  than  approximate,  a  transfer  function  equivalent  can  be 
derived  by  Laplace  transforming  to  obtain: 


YjV<(s) 


(3) 


These  equations  all  indicate  that  quantitatively  the  operator  behaves  as  a  proportional  controller  (pure 
gain)  with  pure  time  delay. 

Although  it  is  easy  to  imagine  from  Fig.  Z  that  the  delayed  proportional  relationship  between  stimulus 
(cause)  and  response  (effect)  obtains,  the  operator's  output  is  not  precisely  given  by  the  operation  of  the 
transfer  function  on  the  stimulus.  Instead,  there  is  a  slight  difference  between  the  delayed  and  amplified 
error  signal  and  the  actual  output.  This  is  a  kind  of  residual,  or  "remnant,"  and  appears  to  have  a  random 
fluctuating  nature  akin  to  a  random  noise.  Since  it  is  generated  by  the  operator,  we  can  consider,  for  the 
moment  at  least,  that  the  remnant  is  operator- induced  noise.  Then,  the  above  equations  in  the  uime  domain 
can  be  rewritten  with  the  remnant  added  and  the  approximate  equality  removed. 

Besides  the  major  system  signals,  i(t),  e(t),  o(t),  and  m(t),  two  other  time  histories  appear  in  Fig.  £. 
These  are  the  differential  EMG,  obtained  by  differencing  muscle  action  potentials  from  the  major  agonist 
and  antagonist  muscles  involved;  and  the  average  EMG,  found  by  summing  the  two  muscle  action  potentials 
(Fig.  1 ) .  The  SEMG  time  trace  thus ^represents  the  operator's  neuromuscular  system  force  output  which  tends 
to  move  the  manipulator,  while  the  T,EMG  indicates  the  average  tension  or  tension  operating  point  within  the 
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Figure  J.  Muscle  Action  Potential  (EMG)  Processing  to  Obtain 
Effective  Muscle  Actuation  Signal  and  Average  Tension 


actuator  system.  In  Fig.  Z  the  AEMG  is  similar  in  general  trend  to  the  operator's  output,  with  a  fairly 
large  nearly-perlodic  wiggle  of  about  10-11  Hz  superimposed.  This  oscillatory  component  is  a  high-frequency 
mode  in  the  neuromuscular  system.  The  ZEMG  is  fairly  constant  but  not  zero,  indicating  that  some  steady- 
state  tension  is  present  and  that  the  agonist  and  antagonist  muscle  groups  are  loaded  against  each  other 
even  when  no  net  force  output  is  acting  on  the  manipulator.  The  neuromuscular  system  will  be  discussed 
later. 

If  the  controlled  element  free  s  is  moved  slightly  into  the  right  half  plane,  a  first-order  divergence 
is  obtained  (fourth  type  in  Table  1).  This  situation  is  unstable  and  requires  constant  attention  to  main¬ 
tain  control,  although  for  moderate  values  of  divergence  rate  the  control  required  is  similar  to  that  for 
Yc  =  Kc/s.  A  typical  example  is  shown  in  Fig.  U,  where  Yc  =  Kc/(s -  Z) .  A  careful  comparison  of  expanded- 
scale  original  records  reproduced  in  reduced  form  as  Figs.  Z'  and  K  revealed  that  the  operator  output  is 
approximately  proportional  to  the  system  errors  for  both,  but  that  the  delay  is  much  less  for  the  diver¬ 
gent  condition,  i.e.,  ‘  0.07  sec  as  opposed  to  tj  -  O.lt  sec.  There  is  also  some  higher-frequency  con¬ 

tent  to  the  system  output  trace,  indicating  that  the  system  output  response  reflects  the  presence  of  a 
closed-loop  v/stem  mode  in  addition  to  the  system  response  to  the  forcing  function. 
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Figure  I*.  Man/Machine  System  Response  Time 
Histories;  Yc  =  Kc/(s-2) 


Let  us  turn  now  to  typical  data  for  a  proportional 
controlled  element,  Kc.  These  are  shown  in  Fig. 

(These  data  ere  for  a  different  subject  than  all  the 
rest,  a  relatively  untrained  non-pilot.)  The  system 
output  follows  the  forcing  function  in  this  case  very 
well,  but  rhen  compared  with  the  Yc  =  K0/s  data  the 
high-frequency  fluctuation  is  more  apparent.  This  is 
because  the  operator-induced  noise  is  passed  directly 
through  the  controlled  element  without  any  filtering 
due  to  the  controlled  element  dynamics.  This  is  also 
evident  from  the  direct  correspondence  between  the 
operator’s  output,  c(t),  and  the  system  output,  m(t). 

Because  the  operator  does  so  well  in  this  kind  of 
system,  the  error  is  very  small.  Consequently,  the 
e(t)  trace  shown  has  been  amplified  to  show  its  nature 
on  a  larger  sc  3e .  If  we  now  examine  the  relation¬ 
ship  between  stimulus,  e(t),  and  response,  c(t),  we 
find  that  unliK.a  the  rate  case,  there  is  little  appa¬ 
rent  resemblance.  So,  the  human  operator  is  clearly 
not  acting  as  a  proportional  controller  operating  on 
the  error.  If,  however,  the  integral  of  the  erbor  is 
compared  with  the  operator  output,  we  see  a  great 
similarity.  Again,  the  output  trace  has  more  higher- 
frequency  "wiggles,"  but  the  main  trend  is  a  near 
duplicate  to  the  integral  of  the  error,  delayed  by 
x,  =  0.)8  sec.  Thus,  we  can  again  consider  a  cause- 
effect  proportional  relationship,  but  this  time 
between  the  integral  of  the  stimulus  and  the  response. 
In  equation  form: 


c(t  +  xi)  =  Kp,  J  e(t)  dt  (5) 


The  equivalent  transfer  function  is: 


Figure  5.  Man/Machine  System  Response 
Time  Histories,  Yc  =  Kq 
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(6) 


This  implies  that  the  operator  acts  as  an  integration 
or  smoothing  filter  on  the  stimulus,  which  is  quite 
different  indeed  from  the  proportional  action  exhi¬ 
bited  in  the  two  previous  situations. 
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For  the  fourth  example,  the  controlled  element 
dynamics  are  Y0  =  K^/s.  For  this  case,  as  shown  in 
Fig.  6,  the  system  output  follows  the  system  forcing 


Figure  6-  Man/Machine  System  Response  Time 
Histories;  Yc  =  Kc/s£> 
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function  well  enough,  in  general  trend,  although  with  some  displacement  to  the  right  lay 

But  here  the  point-to-point  correspondence  is  nowhere  near  as  good  a,  for  the  ^ve  exa^les.  ^nd, 
there  seems  to  be  an  almost  periodic,  function  in  the  system  output,  m(s),  or  e(t)  and  c(t)  ^r  that  matter, 
which  is  not  present  in  the  input,  as  if  the  closed-loop  system  possessed  an  under-damped  relatively  high 
frequency  mode  of  its  own.  Except  for  this,  and  the  Yc  =  Kc/(s-2)  examples,  the  output  has  been  a  faith¬ 
ful  enough  representative  of  the  input  to  suspect  that  operator-lndu .ed  noise  rather  than  any  closed-loop 
system  dynamics  was  affecting  the  match. 

Comparison  of  the  error  and  output  reveals  little  connection  between  the  two  except  a  general  oscilla¬ 
tory  quality.  On  the  other  hand,  if  the  integral  of  the  operator's  output  is  compared  with  I the  error, 
these  time  traces  are  similar.  A  roughly  proportional  relationship  appears  to  exist  ' when  ■ the  error  ii s 
delayed  by  t,  =  O.hJ  sec  and  then  compared  with  the  integral  of  the  operator's  output.  So,  again,  we  can 
demonstrate  1  roughly  proportional  relationship  between  stimulus  and  some  function  of  the  response,  i.e.: 


fc(t)  dt  i  Kpje(t  -  tj) 


(7) 


An  a  transfer  function,  this  will  be: 


its 


=  KpjSe 


(8) 


Thus  the  operator  behaves  as  if  generating  a  first-order  lead  (with  a  break  at  a  very  low  frequency)  on 
the  stimulus  or,  alternatively,  as  if  operating  as  a  proportional  control  on  stimulus  velocity. 

An  intermediate  example  between  the  rate  (Kc/s)  and  acceleration  (Kc/s2)  controlled  elements  is  the 
second-order  system: 


Ke 

s(Ts  +  1 ) 


(8) 


Time  history  examples  for  this  controlled  element  form  with  T  =  l/~  sec  are  given  in  Fig.  7-  The  system 
output  follows  the  forcing  function  quite  well,  with  the  higher  input  frequencies  being  somewhat  accentuated 
in  the  output.  The  operator  output  does  not  correspond  at  all  well  with  the  system  error  as  it  would  if 
the  controlled  element  were  Just  a  rate  control;  and,  in  fact,  the  integral  of  the  show n)  also 

does  not  show  the  type  of  correspondence 
expected  for  the  acceleration  control .  For 
this  intermediate  controlled  element,  the 
operator  output  lagged  by: 
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is  proportional  to  the  error  signal  delayed 
by  about  0.1 6  sec.  The  operator  character¬ 
istics  as  a  transfer  function  would  then  be: 


Kpv(s+3)e”T5,i>  0°) 


|S 


This  result  implies  that  the  operator  develops 
a  lead  which  is  approximately  equal  to  the 
first-order  lag  component  of  the  controlled 
element  dynamics. 

As  the  sixth  and  last  of  the  examples,  a 
more  complex  set  of  controlled  element  dyna¬ 
mics  will  be  considered,  i.e.: 


°o"pu°r 

c  _ 


AEMG  ' 


Average 

EMG 

ZEMG 


Figure  7-  Man/Machine  System  Response  Time 
Histories;  Yc  =  K0/[s(s/3  +  l)] 
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This  is  a  combination  of  the  first-order  lag 
and  first-order  divergence  elementary  con¬ 
trolled  element  dynamics  shown  in  Table  1 • 

The  unstable  root,  X,  creates  an  output  motion 
which  is  always  tending  to  diverge,  so  the 
system  requires  continuous  control  efforts  by 
the  operator.  This  need  for  continuous  opera¬ 
tor  attention  is  not  present  for  Y0  =  Kq, 
where  the  operator's  action  was  motivated 
primarily  by  a  desire  to  have  output  match 
input.  Any  lack  of  operator  action  in  the 
Y0  =  Kq  situation  will  not  incur  the  penalty 
of  a  system  output  divergence.  On  the  other 
hand,  to  the  extent  that  the  operator 
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motivation  in  both  systems  is  error  minimization,  then  his  actions  will  tend  to  be  similar  in  kind,  if  not 
in  detail.  The  actual  signals  in  the  control  system  for  this  controlled  element  are  depicted  in  Fig.  8 
for  c  =  5  and  X  =  1.J>.  Again,  the  output  follows  the  input  rather  well,  with  the  inevitable  delay,  but  the 
operator  output  is  not  in  any  recognizable  way  a  proportional  function  of  the  stimulus.  Yet,  if  the  con¬ 
troller's  output  signal  is  passed  through  a  first-order  lag  with  time  constant,  T  =  1/3  sec,  the  result 
appears  to  be  approximately  proportional  to  the  error  signal  delayed  by  tj,  =  0.20  sec.  Again,  a  propor¬ 
tional  relationship  is  reasonable,  und  the  appropriate  approximate  equations  describing  the  human  operator's 
action  are  given  by  the  transfer  function: 


%  =  Iff}  -  V3  +  5)e“T!iS  (]2) 


The  result  is  very  similar  to  that  obtained  in  Eq.  10  for  the  second-order  controlled  element  with  the 
free  s-1 . 

As  a  final  point  in  this  section,  we  should  emphasize  that  the  simplified  examples  discussed  are  by  no 
means  academic.  They  are  all  approximately  representative  of  various  flight  control  circumstances.  Some 
are  also  applicable  to  other  vehicular  control  situations,  such  as  automobile  driving.  Table  2  presents 
the  exemplary  controlled  element  forms  considered  above  with  a  listing  of  control  situations  which  these 
forms  idealize. 


The  several  examples  described  above  make  plausible  the  first  concept  tnat  some  linear  function  of  the 
operator's  output  can  be  related  approximately  to  another  linear  function  of  the  operator's  stimulus.  The 
second  concept  introduced  is  that  these  relationships  are  different  from  system  to  system;  that  in  fact 
they  depend  explicitly  on  the  controlled  element  dynamics.  Tne  third  concept  illustrated  is  that  the 
actual  relationship  between  stimulus  and  response  contains  a  random  component,  f.t  least  when  the  operator's 
output  is  considered  to  be  predominately  linearly  connected  with  the  input.  These  Kinds  of  data  form  an 
empirical  basis  f^r  the  description  of  human  operators  as  quasi-linear  systems  comprising  describing  func¬ 
tion  and  remnant  components  to  characterize  human  behavior.  This  will  be  formally  pursued  in  the  next 
chapter.  The  examples  taken  in  concert  can  be  explored  more  fully  in  an  attempt  to  develop  some  approx - 
mate  laws  of  human  behavior.  This  will  be  done  below. 

B.  AH  APPROXIMATE  CROSSOVER  MODEL  FOR 
HUMAN  DYNAMIC  OPERATIONS  IN  3IN0LE- 
LOOP  COMPENSATORY  SYSTEMS 
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Figure  8.  Man/Machine  System  Response  Time 
Histories;  Y0  =  Kc/[(s +  3)(s- 1 .5)] 


The  human  operator  transfer  functions 
developed  for  the  controlled  element  set 
investigated  above  are  compiled  in  Table  3. 

As  indicated  there,  the  human’s  transfer 
function  is  different  for  each  controlled 
element,  but  the  open-loop  system  function, 
i.e.,  YpYc,  is  essentiallj  the  same  in  form 
for  '‘our  of  the  six  cases.  In  fact,  in  fre¬ 
quency  response  terms,  all  the  cases  fall 
into  the  same  category  for  values  of  fre¬ 
quency  greater  than  X. 

Although  the  form  is  common,  the  compo¬ 
nents  are  quite  different.  As  will  be  seen 
in  subsequent  chapters,  the  "crossover 
frequency, "  o^,,  which  always  contains  the 
operator  gain,  Kp,  as  a  factor,  is  different 
for  each  of  the  controlled  elements  consid¬ 
ered.  We  have  not  thus  far  applied  the 
appropriate  analytical  tools  to  determine 
Kp  from  the  data  presented,  so  the  ramifica¬ 
tions  of  this  point  will  await  more  complete 
development  subsequently.  For  the  examples 
we  have  described  above,  the  effective  time 
delay  has  been  approximately  determined  from 
the  time  traces  by  superimposing  the  two 
nearly  proportional  time  history  quantities 
and  shifting  them  about  until  a  maximum 
degree  of  overlay  occurs.  This  is,  of 
course,  an  eyeball  rather  than  mathemati¬ 
cally  based  determination.  Nonetheless,  the 
effective  time  delays  determined  show  an 
interesting  ordering  across  the  five  con¬ 
trolled  elements  for  which  we  have  data  for 
the  same  subject.  Proportional  control 
action,  which  is  exhibited  when  the  con¬ 
trolled  element  transfer  function  itself  is 
an  integrator,  can  be  taken,  temporarily, 
as  a  baseline.  When  the  operator  has  to 
generate  low-frequency  lead,  a  larger  delay 
is  incurred.  Even  when  only  moderate  amounts 
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TABLE  2 

EXAMPLE  APPLICATIONS  OP  IDEALIZED  CONTROLLED  ELEMENT  FORMS 


CONTROLLED  ELEMENT 
POEM 

AEROSPACE  CONTROL 

AUTOMOBILE  CONTROL 

Attitude  (pitch  or  roll)  control  with 

K0 

an  attitude  corrmand,  attitude  hold 
coOT.and  augmentation  system  (CAS) 

Speed  control 

Kc 

Attitude  (pitch  or  roll)  control  with 

Heading  control  at  low  to 

s 

a  rate  command  CAS  or  with  a  damper 
stability  augmentation  system  (SAS) 

moderate  speeds 

K„ 

Roll  attitude  control  of  conventional 

aircraft  with  roll  subsidence  time 

s(Ts  +  1 ) 

constant,  T 

Attitude  (pitch  or  roll)  control  of  a 

Longitudinal  position  con- 

Kc 

spacecraft  with  damper  off 

trol 

Limiting  case  of  roil  attitude  control 

Lateral  path  control  with 

S2 

of  conventional  aircraft  as  roll 

only  positional  cues 

subsidence  time  constant  becomes 

(e.g.,  local  lane  markers) 

very  large 

and  little  preview 

Kc 

Heading  control,  with  rudder,  of  con- 

S  -  X 

ventional  aircraft  spiral  mode 

Kc 

Pitch  attitude  control  of  unstable 

Heading  control  of  oversteer- 

(s  +  a)(s  -  x) 

short  period 

ing  car  above  critical  speed 

TABLE  3 

SUMMARY  OP  HUMAN  OPERATOR  APPROXIMATE  CHARACTERISTICS 
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of  operator  lead  are  present,  as  in  the  Yc  =  K<./s(s  +  a)  and  Y0  =  Kc/(s  +  a){s  —  X)  cases,  the  delay  is 
greater  than  with  proportional  control.  Qualitatively,  at  least,  it  thus  appears  that  low- frequency  lead 
equalization  of  the  system  which  is  accomplished  by  the  operator  incurs  a  coot  which  can  be  measured  in 
terms  of  effective  time  delay  introduced.  In  the  other  direction,  the  introduction  of  a  mild  divergence 
in  the  controlled  element  resulted  in  a  smaller  delay  for  the  operator. 

In  subsequent  chapters,  more  elabors  te  mathematical  and  analysis  procedures  will  he  applied  to  many 
more  complex  situations.  As  a  consequence,  the  operator's  dynamic  behavior  will  be  quantified  in  much 
greater  depth  and  detail  than  that  presented  above.  Yet,  when  all  is  said  and  done,  a  common  thread 
through  almost  all  situations  which  we  shall  c  insider,  and  the  only  thread  which  has  any  pretensions  to 
general  applicability,  is  the  form  exposed  in  Table  5  above.  This  is  the  "crossovir  model,"  which  relates 
the  operator  and  controlled  clement  transfer  c  laracterlstics  by  the  equation: 


Yp(jo)Y0(jaO 


rn-e 


-jwe 


03) 


The  reason  :\r  the  "crossover"  appellation  is  connected  with  the  model's  frequency  range  of  validity, 
which  was  introduced  above  in  connection  with  the  divergent  controlled  elements,  and  which  will  be  elabo¬ 
rated  on  extensively  throughout  the  report. 


C.  MAIHEMATIC3  OF  THE  CROSSOVER  MODEL  —  INTRODUCTION 
TO  MAH  ''MACHINE  SYSTEMS  AHALYSXS  AHD  PREDICTION  OF 
HUMAN  OPERATIONS  IN  COMPENSATORY  SYSTEMS 


One  of  the  results  of  the  first  two  sections  ir.  this  chapter  is  the  crossover  model  as  an  approximate 
quantitative  description  of  man/machine  system  dynamics  for  single-loop  compensatory  systems.  The  model 
does  net  distinguish  explicitly  betw-eer.  controlled  element  and  operator  characteristics,  although  it  is 
plain  from  what  has  been  presented  that  the  two  parameters,  a0  and  -ce,  depend  on  controlled  element  dyna¬ 
mics.  It  will  later  turn  out  that  these  quantities  are  dependent  also  on  other  control  system  tnsu  vari¬ 
ables  and  on  experimental  situation  variables  which  have  an  impact  on  the  operator.  Discussion  of  these 
possible  dependencies  will  be  deferred  while  the  system  dynamics  and  performance  for  the  crossover  model 
as  a  control  system  is  derived  and  discussed.  To  accomplish  this,  c*  nventional  methods  of  feedback  systems 
analysis  will  be  used. 

The  purpose  of  this  exercise  is  twofold.  First,  it  serves  a  tutorial  function  in  presenting,  by  way 
of  an  elementary  example,  some  sit.  pie  methods  which  are  useful  in  the  analytical  treatment  of  manual  con¬ 
trol  systems.  Second,  all  of  the  generalized  data  on  the  dynamics  and  performance  of  the  crossover  model 
developed  can  be  used  to  make  estimates  or  interpretations  of  manual  control  system  properties  by  specializ¬ 
ing  the  results  with  par.,  ular  assignments  of  t  and  u^.  The  answers  obtained  will  be  valid  for  those  many 
circumstances  in  which  the  human  controller  behaves  in  a  fashion  similar  to  that  depicted  by  the  crossover 
model.* 

When  the  loop  is  closed  about  the  open-loop  crossover  model,  as  depicted  by  the  Fig.  9a  block  diagram, 
the  result  is  the  system  dynamics  summarized  in  the  system  survey  of  Fig.  9b  and  c.  The  survey  shows  a 
number  of  related  plots.  The  most  familiar  are  probably  the  G( jo.)  Bode  plot  (or  open-loop  ju>  Bede)  and 
the  conventional  root  locus.  As  shown  in  Fig.  the  ju.  Bode  diagram  gives  the  open-loop  amplitude  ratio 
and  phase  in  dB  and  degrees,  respectively,  versus  the  normalized  frequency,  tu>,  on  a  logaritliraic  scale. 

The  amplitude  ratio  plot  is  a  straight  line  with  a  slope  of  -20  dB/deeade.  The  phase  lag,  while  increas¬ 
ing  at  an  ever-increasing  rate  in  these  coordinates,  is  linear  with  frequency.  This  diagram  represents 
the  open-loop  frequency  response  characteristics,  and  can  be  manipulated  (e.g.,  using  a  Nichols  chart. 

Ref.  3)  and  interpreted  to  determine  the  closed-loop  frequency  response  properties.  The  conventional  root 
locus  shown  in  Fig.  9b  indicates  the  closed-loop  roots  of  the  system.  Each  closed-loop  root  corresponds 
to  a  different  open- loop  gain,  which  is  the  basic  parameter  along  the  loci.  Damping  ratio  is  directly 
indicated  for  a  particular  root  vector  by  its  angle  made  with  the  negative  real  axis;  this  angle  is  cos-  £. 
On  Fig.  9b  these  ^  values  are  called  out  rather  than  the  gains  as  parameters  along  tne  complex  locus. 

The  root  locus  and  .ja  Bode  are  directly  correlated  at  that  value  of  gain  for  which  the  system  is 
neutrally  stable.  At  this  point,  the  phase  angle  on  the  Bode  plot  is  — 180  deg,  and  the  normalized  gain, 
to;c,  is  rt/2.  The  reference  0  dB  line  shown  in  Fig.  9c  is  constructed  for  this  gain  value  as  it  corre¬ 
sponds  to  neutral  stability.  The  underdampod  quadratic  locus  branches  (only  the  upper  one  is  shown)  of 
the  root  locus  cross  the  jru.  axis  at  this  value  of  gain  on  their  way  from  the  left  to  the  right  half  plane. 

The  closed-loop  system  dynamics  relating  system  output  and  error  to  the  forcing  function  are  given  by: 


Gmi(s) 


M(s)  _  0(s) 

TCsJ  “  1  +  G(s) 


-TS 

<«ce 


s  +  u^e 


■TS 


—  1  when  |G(s)|  » 


(1*0 


*In  some  cases  the  operator-induced  noise  or  remnant  is  large  enough  to  be  of  major  consequence.  The/ 
answers  provided  by  the  crossover  model  may  still  be  valid,  but  only  to  provide  the  linearly  correlated  , 
(with  forcing  functions  or  disturbances)  components  of  signals  within  the  control  loop.  The  remnant  intro' 
duces  another  source  of  system  excitation  which  must  be  taken  into  account. 
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Perhaps  the  most  significant  property  of  feedback  systems  is  that  which  obtains  when  the  open-loop  transfer 
characteristic  is  much  larger  than  1,  for  then  the  system  output  is  almost  exactly  equal  to  the  system  input 
and  the  system  error  is  very  small.  From  examination  of  the  G(ju>)  Bode  plot  it  can  be  seen  that  this  occurs 
at  low  frequencies  and  is,  of  course,  the  reason  that  the  system  output  and  system  forcing  function  were 
similar  in  a'  the  examples  given  thus  far.  At  high  frequencies  |G( Ja>) |  «  1,  so  the  closed-loop  relation¬ 
ship  between  M(j<o)  and  I(Jui)  is  substantially  the  same  as  the  open-loop,  i.e.,  the  feedback  loop  is  effec¬ 
tively  inoperative.  For  the  crossover  model,  the  frequency  which  divides  these  two  regimes  of  near-ideal 
following  of  the  forcing  function  and  little  or  no  feedback  action  is  the  crossover  frequency,  cue.  In 
Fig.  9c  this  is  the  intersection  of  the  0  dB  line  established  for  a  particular  gain  with  the  G(jcu)  Bode 
plot.  The  name  "crossover  frequency"  comes  from  this  crossover  intersection  of  the  0  dB  line  by  the  open- 
loop  frequency  response  characteristic.  For  stable  operation  of  the  system  the  normalized  crossover  fre¬ 
quency,  tcoc,  can  range  free.  0  to  ,t/2. 

When  Tu'c  is  relatively  small  compared  to  it/2,  then  tcdc  is  also  tantamount,  to  the  closed-loop  system 
"bandwidth"  (the  frequency  at  which  the  output  amplitude  is  3  dB  less  than  the  amplitude  of  an  input  sinu¬ 
soid)  For  higher  crossover  frequencies  this  direct  equivalence  between  tio,.  and  system  bandwidth  uegrades 
because  of  a  peaking  in  the  closed-loop  frequency  response  near  the  crossover  frequency,  cut  even  for  this 
kind  of  a  system  rcoc  is  always  equal  to  or  less  than  the  bandwidth  and  thus  provides  a  lower  bound. 

While  the  conventional  G(ju,)  Bode  plot  is  useful  in  frequency  domain  descriptions  to  illustrate  the 
fundamental  properties  of  feedback  as  it  affects  output/input  relationships,  the  root  locus  best  emphasizes 
the  closed-lcop  system  roots.  The  root  locus  plot  shown  in  Fig.  9b  indicates  that  the  root  which  starts  at 
the  origin  for  open-loop  (zero  gain)  conditions  progresses  further  into  the  left  half  plane  as  the  gain  tu'c 
is  increased.  At  the  point  to  =  ~i  (for  which  tcuc  =  l/e),  this  branch  of  the  locus  meets  with  the  first  of 
an  infinity  of  branches  present  because  of  the  e“TS  term  (the  other  branches  are  not  shown).  The  quadratic 
formed  by  the  two  branches  then  increases  in  undamped  natural  frequency  and  decreases  in  damping  as  gain 
is  increased  until  neutral  stability  occurs  at  the  gain  tu^.  =  it/2. 


a)  Block  Diagram 
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Figure  9-  System  Survey  for  Crossover  Model 
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The  connections  between  the  dominant  closed-loop 
second-order  mode  and  open-loop  parameters  are 
readily  developed,  '’’he  damping  ratio  can  be  related 
directly  to  the  phase  margin  from  the  6(ju»)  Bode 
plots.  For  a  particular  open-loop  gain,  wu,;,  the 
phase  margin  is  given  by: 

~  tjoc  (itc) 


The  relationship  between  closed-loop  damping  ratio, 
5C,  and  phase  margin,  <j^,  is  shown  in  Fig.  10  (with 
open-loop  gain  as  an  alternate  ordinate).  This  could 
also  be  expressed  in  terms  of  the  gain  margin, 

|GMl<iB  >  which  is’ 

[i  ~(2%/«)r’L 


(rod)  (deg) 


Closed -Loop  Dorr.piog  Ratio  ,  £c 


It  is  apparent  from  the  Bode  root  locus  (Fig.  9a) 
that  for  the  values  of  normalized  gain,  tu)0,  which 
result  in  an  underdamped  closed-loop  system  dominant 
mode,  the  logarithm  of  the  closed-loop  undamped 
natural  frequency  of  the  quadratic  modes  is  linearly 
related  to  the  logarithm  of  the  crossover  frequency 
(or  gain  in  dB).  The  actual  relationship  can  be 
shown  to  be: 


Thus,  the  relationship  between  closed-loop  undamped  Figure  11.  Closed-Loop  Undamped  Natural 

natural  frequency  and  open-loop  gain  is  a  straight  Frequency  as  a  Function  of 

line  in  logarithmic  coordinates.  This  is  illustrated  Crossover  Frequency 

in  Fig.  11  for  the  range  of  tu^  from  1  (at  the  ren¬ 
dezvous  point  on  the  root  locus)  to  n/£  at  neutral 
stability. 

The  same  information  available  on  the  conventional  root  locus  is  also  present  in  the  so-called  Bode 
root  iocus.  This  comprises  two  elements:  the  "siggy  plot"  of  I G( -a ) ] versus  xoj  and  the  complex  branch 
plotted  as  a  dotted  lino  versus  t|s|.  The  direct  correspondence  between  the  two  root  locus  plots  is  indi¬ 
cated  on  Figs.  9b  and  9c  by  the  labeling  of  the  branches  1,  2,  and  J.  On  the  Bode  root  locus,  £c  is  a 
parameter  and  gain  is  the  ordinate j  whereas  on  the  conventional  root  locus  gain  is  a  parameter.  The  Bode 
root  locus  places  on  one  common  figure  both  the  frequ :ncy  response  and  closed-loop  pole-zero  relationships. 

In  system  dynamic  operations,  the  system  output  response  will  contain  two  components.  The  first  is  the 
forced  response,  i.e.,  an  output  which  derives  from  the  system  operating  upon  the  forcing  function.  The 
second  output  component  is  the  natural  modes  of  the  system  as  transients  excited  by  the  forcing  function. 
For  the  crossover  model  these  modes  are  an  overdarnped  second  order  (i.e.,  two  first  orders)  for  tup  less 
than  l/e,  and  on  oscillatory  second  order  from  that  gain  to  the  stability  limit,  xuc  =  n/2.  When  the 
system  dynamics  are  relatively  well  damped  and  the  forcing  function  bandwidth  is  low  compared  with  the 
crossover  frequency,  then  the  system  output  will  very  closely  resemble  the  system  forcing  function  alone. 
When  the  closed-loop  system  dynamics  are  more  oscillatory,  the  presence  of  t’  i  natural  modes  will  also  be 
seen  in  the  system  output  because  this  component  then  does  noo  rapidly  damp  out  and  is  constantly  being 
excited  by  the  forcing  function. 

If  we  now  view  some  of  the  example  cases  studied  in  the  first  section  we  will  recall  two  forms  of 
following  responses.  For  several  of  the  systems  the  system  output  follows  the  forcing  function  quite  pre¬ 
cisely  (except  for  a  very-high-frequency  component  in  the  output  which  is  not  explained  by  the  crossover 
model).  For  some  of  the  other  systems,  the  output  followed  the  forcing  function  fairly  well,  but  another 
mode  also  appeared  to  be  present.  We  can  identify  this  with  the  basic  closed-loop  mode  when  it  is  oscilla¬ 
tory  in  nature.  Thus,  for  those  cases  which  show  this  quality  in  the  responses,  the  closed-loop  damping 
ratio  was  on  the  low  side,  corresponding  to  a  high  gain  and  a  low  phase  margin. 

When  the  forcing  function  is  a  stationary  random  or  a  harmonic  process,  it  can  be  characterized  by  a 
power  spectral  density,  Oji-  Then  the  system’s  steady-sta‘e  performance  in  response  to  such  a  forcing 
function  will  be: 


where  Gei  Is  the  transfer  characteristic  between  the  forcing  function  and  the  error  (Eq.  ) .  For  a  rec¬ 
tangular  power  spectral  density  input  with  bandwidth  u>i  and  variance  of  this  becomes: 


1  pi  uj*  dto _ 

“i  Jq  a"— 2nac  sin  at 


09) 


This  integral  is  not  readily  evaluated  analytically.  However,  if  sin  wx  at,  which  is  an  approximation 
good  when  u^t  «  i ,  then: 
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(20) 


If  the  tanh 


is  expanded  and  carried  only  to  the  second  term,  the  extremely  simple  "1/3  Law"  results, 


i  .e. 


(21) 


This  gives  remarkably  good  results  for  a  large  number  of  practical  control  systems.  When  more  precise 
answers  are  needed  the  exact  results,  obtained  by  integrating  Eq.  19,  shown  in  Fig.  12  can  be  used. 


Figure  12.  Crosssover  Model  Steady-State  Performance 
for  Inputs  with  Rectangular  Spectral  Densities 


16 


CHAPTER  III 

X SY  VARIABLES  AMD  MEASUREMENTS 


A.  KEY  VARIABLES  WHICH  AFFECT  THE  PILOT’S  BEHAVIOR 

As  our  discussion  of  the  origins  of  human  dynamics  theory  has  indicated,  engineering  needs  provided 
the  impetus  and  engineering  techniques  the  methods  for  much  of  the  early  research.  The  nature  and  char¬ 
acteristics  of  this  approach  is  exemplified  qualitatively  by  the  presentation  of  Chapter  II.  In  fact,  the 
use  there  of  a  single  subject,  with  the  tacit  assumption  that  single-run  segments  are  representative 
responses,  exposes  the  obvious  need  for  detailed  consideration  of  intrasuoject  and  ’ntersubject  differ¬ 
ences.  Here,  psychological  thinking  has  had  a  major  impact  in  calling  attention  to  the  many  variables 
which  might  feet  the  intra-  and  inter-subject  variability  of  human  dynamic  performance  and,  by  so  doing, 
influencing  the  design  of  experiments.  A  strong  traditional  vie  ./point  in  experimental  psychology  has 
been  oriented  to  the  description  of  behavior  in  terms  of  discrete  events  rather  than  to  the  interacting 
continuous  description  demanded  by  a  closed-loop  analysis.  As  a  consequence,  much  effort  in  psychomotor 
research  in  psychological  laboratories  was  devoted  to  specifying  stimulus  and  response  pairs  in  motor 
behavior  under  cons-rainea  conditions.  Response  latencies  and  simple  performance  measures  were  favored 
measurements,  and  tne  description  and  parameters  of  learning  were  common  objects  of  study.  It  required 
thv  introduction  of  closed-loop  analysis,  which  even  in  simplified  form  provides  a  greatly  enriched  char¬ 
acterization  of  human  psychomotor  properties,  to  enable  the  wealth  of  sophisticated  behavior  which  the 
human  can  generate  tJ  be  revealed,  understood,  and  applied  to  the  rational  design  of  control  tasks.  Two 
examples  are  the  analysis  of  what  were  previously  discrete  reaction  time  measures  into  components  arising 
from  the  plant  dynamics,  the  neuromusculature,  and  the  processing  time  requirements  of  the  human  controller 
and  the  description  of  the  exploitation  of  the  human  ability  to  discern  temporal  patterns  of  signals  in 
improving  the  performance  of  a  variety  of  skilled  tasks. 

Although  human  controller  variabilities  can  be  magnified  or  minimized  depending  on  what  measurements 
are  taken  and  where  the  measurement  points  lie,  there  is  a  need  to  specify  appropriate  ranges  of  reliable 
application  for  human  dynamic  findings.  To  fill  this  need  we  begin  by  classifying  the  plethora  of  physical 
psychological,  physiological,  and  experimental  effects  which  can  affect  the  human's  relevant  outputs.  This 
organization  is  illustrated  in  Fig.  13  where  the  variables  —  in  some  cases  really  parameters  —  are  sub¬ 
sumed  under  four  categories  as  follows. 


Task  Variables 


Environmental  Variables  .  c  ir.  Procedural  Variables 


Figure  O*  Variables  Affecting  the  Filot/Vehiele  System 


1  ■  T»ik  Variables 

Task  variables  comprise  all  the  system  inputs  and  those  control  system  elements 
external  to  the  pilot  which  entr.r  directly  and  explicitly  into  the  pilot's  control 
task.  Stability  of  the  closed-loop  system  is  always  a  necessary,  though  not  suffi¬ 
cient,  control  strategy.  Consequently,  the  pilot's  dynamics  are  profoundly  affected 
by  the  display  and  controlled  element  dynamics,  because  his  properties  must  be  adapted 
to  provide  the  necessary  loop  stability.  The  general  nature  of  these  adaptations  has 
been  illustrated  in  Chapter  II.  The  characteristics  of  the  other  task  variables, 
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i  .e .,  disturbance  inputs  and  command  inputs  related  to  the  mission  and  control 
strategy,  also  exert  direct  influences  on  the  pilot  dynamics,  although  their 
effects  are  more  in  the  nature  of  adjustment  and  emphasis  than  of  change  in 
fundamental  dynamic  form. 

These  variables  constitute  an  enormous  range  of  possible  conditions.  As  a 
practical  convention  interest  is  often  focused  on  a  limited  number  of  input  sig¬ 
nals  having  well-defined  amplitude  distribution,  shape,  and  bandwidth,  on  a 
selection  of  controlled  element  dynamics  which  represent  or  idealize  practical 
vehicles,  and  on  a  restricted  range  of  both  functional  and  physical  displays.  The 
abstraction  of  these  task  variables  was  accomplished  by  a  combination  of  engineer¬ 
ing  insights  (the  variables  were  similar  to  gust  disturbance,  the  dynamics  were 
similar  to  aircraft  dynamics,  etc.)  with  the  need  to  make  accurate  measurements 
over  a  usefully  broad  range  of  frequencies  (the  use  of  subdued  high-frequency 
signals,  the  selection  of  appropriate  run  lengths,  etc.). 

2.  Environmental  Variables 

The  state  of  the  environment  external  to  the  pilot  is  shown  in  Fig.  1 J  as  the 
vector  e.  Included  as  components  of  this  vector  are  such  factors  as  ambient 
illumination,  vibration,  temperature,  acceleration  (to  the  extent  that  this  is 
superimposed  on,  rather  than  controlled  by,  the  pilot),  noise,  ambient  atmosphere, 
etc. 

3 •  Procedural  Variables 

The  procedural  variables,  denoted  by  the  vector  jt,  include  such  aspects  of 
experimental  procedure  as  instructions  and  background  indoctrination,  training 
schedule,  order  of  presentation  of  trials,  and  so  forth.  For  those  experiments 
for  which  the  subjects  are  experienced  engineering  test  pilots,  a  carefully 
planned  indoctrination  briefing  can  serve  to  heighten  their  motivation  as  well 
as  to  enable  them  to  generate  an  appropriate  mode  of  response.  In  many  cases, 
the  experimental  control  device  is  structurally  a  highly  abstracted  version  of 
the  actual  aircraft  being  simulated  and  a  careful  indoctrination  can  help  avoid 
a  situation  in  which  the  pilot's  control  responses  are  more  appropriate  to  a  pin¬ 
ball  machine  than  to  the  flying  of  an  aircraft.  There  is,  however,  a  danger  that 
in  motivating  a  subject  he  will  become  so  involved  in  the  experiment  that  despite 
his  intentions  to  the  contrary  he  will  influence  the  results  to  conform  with  what 
he  perceives  are  the  experimenter's  hopes.  Instructions  and  indoctrination  serve 
either  explicitly  or  by  means  of  the  mission  to  establish  the  performance  cri¬ 
terion  which  the  subject  will  use. 

The  experimental  design  and  tne  statistical  analyses  used  can  serve  to  obscure 
actual  effects  or  "reveal"  behavioral  phenomena  which  are  in  truth  artifacts  of 
improper  technique.  Meticulous  attention  must  be  paid  to  counterbalancing  time- 
dependent  effects  on  subjects  and  to  making  assumptions  explicit.  In  view  of  the 
generally  small  number  of  subjects  and  runs  generated  by  human  dynamics  experi¬ 
ments,  it  is  necessary  to  develop  techniques  for  controlling  extraneous  variables 
rather  than  depending  on  randomization  over  a  large  number  of  subjects  and  condi¬ 
tions  in  the  experimental  design. 

4.  Pilot-Centered  Variables 

The  operator-centered  variables,  denoted  by  the  vector  a,  include  the  char¬ 
acteristics  the  pilot  brings  to  the  control  task:  training7  motivation,  "set"  to 
respond,  physical  condition,  etc.  Many  of  these  factors  are  difficult  to  quantify 
in  terms  meaningful  to  a  given  experiment.  They  can,  however,  at  least  be  quali¬ 
tatively  graded  by  pretest,  interview,  etc.,  or  controlled  or  modified  by  proce¬ 
dures  (therefore  there  is  some  interaction  between  n  and  a).  "Set"  to  respond, 
for  example,  as  established  by  a  particular  set  of  experimental  experiences  can 
be  compensated  for  by  counterbalancing  the  order  in  which  subjects  are  tested. 

The  subject's  performance  criteria  are  another  pilot-centered  variable  which  can 
be  modified  by  procedures. 

B.  DESCRIPTION  0?  HUMAN  PILOT  BEHAVIOR 

The  most  obvious  aspect  of  human  dynamic  behavior  in  a  control  tass  is  the  pilot's  control  actions 
within  that  task.  There  exist  precedents  from  the  analysis  of  inanimate  systems  for  specifying  which  con¬ 
trol  actions  should  be  measured  for  a  description  of  the  controller's  behavior.  Such  measurements  do  not 
tell  the  complete  story  for  a  human  controller,  since  associated  with  the  control  actions  are  physiological 
and  psychophysiological  outputs,  the  vectors  <£  and  £.  These  include  status  indicators  of  the  human's 
internal  environmental  control  systems,  such  as  respiratory  rate  and  volume,  heart  rate  and  blood  pressure, 
rate  of  sweating  and  body  temperature,  etc.,  as  well  as  such  highly  structured  but  nonetheless  subjective 
indications  of  workload  and  pilot  behavior  as  Cooper-Harper  pilot  ratings. 

Both  the  measures  derived  from  inanimate  control  device  description  and  the  peculiarly  human  outputs 
listed  above  can  be  used  to  provide  operational  definitions  for  a  variety  of  verbal  concepts  commonly 
associated  with  human  behavior.  Skill,  for  example,  is  a  concept  which  has  been  described  in  such  intui¬ 
tive  terms  as  "sequence  of  deftly  timed  responses"  and  "the  outstanding  character  of  rapid  adaptation" 

(Ref.  39) .  The  availability  of  dynamic  descriptions  of  human  control  actions  enable  us  to  quantify  "deftly 
timed"  in  a  fashion  not  otherwise  possible-  Similarly,  the  human  ability  to  adapt  can  be  reduced  to 
readily  quantifiable  changes  in  the  mathematical  form  of  the  description  of  the  control  actions. 
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To  make  these  general  statements  more  concrete  presume,  for  the  moment,  that  the  crossover  model 
developed  in  Chapter  II  was  a  complete  description  of  operations  in  a  single-loop  manual  control  system. 
Then,  the  equations  of  motion  for  a  specific  system  would  be: 


e(t)  =  i(t)  -  m(t) 
m(t)  =  o)ce(t  —  t) 


(22) 


The  first  of  these  equations  app] ies  for  all  compensatory  situations;  whereas  the  second  requires  that  the 
crossover  model  bo  valid.  The  corresponding  open-loop  transfer  function  is: 


The  human  pilot's  adaptation  to  controlled  element  dynamics  is  implicit  in  this  relationship,  i.e.,  for  a 
particular  set  of  controlled  element  dynamics  defined  by  Yc(s),  the  human  will  adopt  a  transfer  function: 


Vs)  = 


mee~TS 

SYq(s) 


The  general  form  of  the  hvjnan's  response  would  thus  be  determined  by  the  specifics  of  Yc,  and  a  change  in 
this  task  variable  evokes  a  change  in  Yp  such  that  the  crossover  model  open-loop  transfer  function  is  con¬ 
served.  The  effects  of  changes  in  other  task  variables,  or  in  the  environmental  and  operator-centered 
variables,  will  be  to  modify  a>c  or  t.  These  and  t  modifications  are  themselves  the  quantification  of 
changes  or  differences  in  key  variables  on  the  pilot's  control  action.  In  measuring  the  effects  of  train¬ 
ing,  for  instance,  ioq  increases  with  trials  until  stable  conditions  are  obtained  for  that  particular  sub¬ 
ject  and  set  of  constant  task  and  environmental  variables.  For  circumstances  in  which  the  l/3  taw  applies 
(Eq.  21 ),  an  increase  in  <uc  with  trials  will  be  reflected  by  a  hyperbolic  decrease  in  the  system  rms  error. 
Thus,  in  the  context  of  this  hypothetical  example  in  which  the  crossover  model  is  a  total  representation 
of  human  control  actions,  the  human's  control  action  will  be  quantified  in  terms  of  the  crossover  model 
gain,  <cc,  and/or  effective  time  delay,  t. 

The  human  pilot's  actions  are  unfortunately  not  as  readily  described,  in  general,  as  this  hypothetical 
example.  For  one  thing,  the  effects  of  environmental  stresses  and  operator-centered  variables  may  mani¬ 
fest  themselves  in  the  £  and  £  vectors,  as  well  as  in  changes  in  the  control  actions.  Further,  the  human 
pilot  is  a  multi-input, "“multi-output  device  of  enormous  complexity  rather  than  a  single-channel  control 
mechanism.  Nonetheless,  it  should  be  plain  from  both  the  hypothetical  and  real-data  examples  discussed 
thus  far  that  a  great  amount  of  verbally  expressed  descriptive  material  and  concepts  can  be  quantified  and 
reduced  to  operationally  usab.le  terms  by  virtue  of  a  control  engineering  approach. 

A  complete  and  detailed  descriptor  of  a  human  in  this  context  is  one  which  relates  control,  physiologi¬ 
cal,  and  psychophysiological  outputs  to  control  and  environmentally-derived  inputs.  For  a  large  number  of 
practical  cases,  however,  the  £  and  £  outputs  are  of  secondary  interest  and  the  e  inputs  are  parameters 
over  the  measurement  period.  When  attention  is  focused  on  control  actions  as  the  major  interest,  and 
as  shall  be  done  henceforth,  the  pilot's  control  activity  is  capable  of  being  described  as  a  short-time 
stationary  process  for  a  very  large  number  of  circumstances. 

When  the  key  variables  are  fixed  and  the  signals  in  the  control  loop  are  approximately  time  stationary 
over  an  interval  of  interest,  the  pilot/vehicle  system  can  be  modeled  as  a  quasi-linear  system.  Quasi¬ 
linearization  is  one  of  the  most  fruitful  approaches  to  the  description  of  those  nonlinear  time-varying 
systems  in  which  the  relationships  between  pertinent  measures  of  system  input  and  output  signals  have  some 
linear  correlation  in  spite  of  the  possible  existence  of  nonlinearities  and  short-term  (relative  to  the 
observation  interval)  time  variations.  The  description  is  valid  for  only  the  specific  situation;  when 
conditions  (e.g.,  inputs)  are  changed  a  different  quasi-linear  description  is  needed. 

In  a  quasi-linear  system  the  response  for  a  given  input  is  divided  into  two  parts  —  describing  function 
components  which  correspond  to  the  responses  of  equivalent  linear  elements  driven  by  that  input,  and  a 
"remnant"  component,  which  represents  the  difference  between  the  response  of  the  actual  system  and  an  equi¬ 
valent  system  based  on  the  linear  element  (Ref.  !»0).  Quasi-linear  models  consisting  of  describing  function 
plus  remnant  descriptions  for  random-appearing  inputs  and  disturbances  have  been  the  basis  for  the  vast 
majority  of  man/machine  systems  analyses  and  have  also  received  the  lion's  share  of  experimental  effort. 

The  most  important  class  of  situations  in  closed-loop  control  of  aircraft  are  compensatory  tasks  in 
which  the  pilot  acts  on  displayed  error  quantifies,  £,  between  desired  command  inputs,  i,  and  comparable 
vehicle  output  motions,  m,  to  produce  control  actions,  c.  This  class  is  illustrated  in~Fig.  1h.  In  this 
block  diagram  the  dynamics  of  the  equivalent  controlled~element  and  displays  are  described  by  a  matrix  of 
transfer  functions,  jY0(j<jo)}.  The  signals  in  this  general  block  diagram,  i,  e,  e,  and  m,  as  well  as  the 
remnant,  ne  (considered  as  a  quantity  injected  at  the  pilot's  input)  are  all,  in  general,  vector  quanti¬ 
ties.  Finally,  the  transfer  characteristics  of  the  pilot  are  represented  by  the  matrix  of  quasi-linear 
describing  functions,  {Yp}.  The  describing  functions  and  remnant  depend  explicitly  on  the  task  variables, 
as  noted  in  the  functional  notation.  (While  an  explicit  functional  dependence  is  not  shown,  the  remnant  and 
transfer  characteristics  are  also  functions  of  the  operator-centered,  procedural,  and  environmental  variables.) 
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Figure  l4.  Quasi-Linear  Paradigm  for  the  Human  Pilot 


C.  MEASUREMENT  FUNHAMEHIALS 

The  detailed  nature  of  typical  measurements  in  manual  control  systems  is  most  easily  introduced  by 
simplifying  the  general  block  diagram  of  Fig.  14  to  a  single -loop  situation.  Such  a  figure  appears  at  the 
top  of  Table  4.  Here,  the  pilot's  properties  are  divided  into  actuation-load  dynamics  and  sensory  and 
equalization  dynamics  blocks.  The  manipulator  task  variable  is  intrinsically  included  in  the  former.  The 
signals  throughout  the  loop  are  characterized  in  several  ways:  as  a  time  function,  x(t);  as  a  mean-squared 
value,  xTtT)  as  an  amplitude  distribution  probability  density,  px,  where  px  dx  =  Prob(x^  <  x  <  Xg);  and 
as  a  power-spectral  density,  «w(tu).  The  time  functions,  i,  e,  c,  and  m,  are  ordinarily  available;  whereas 
special  means  (e.g.,  see  Fig.  3)  are  needed  to  obtain  EMG  signals  proportional  to  the  pilot’s  force,  F(t), 
which  drives  the  actuation  and  load  dynamics. 

For  stationary  situations  the  forcing  function  signal  i(t)  may,  in  principle,  have  two  components, 
random  and  periodic.  As  power-spectral  densities  these  correspond  to  a  continuous  power-spectral  density 
function,  and  a  sum  of  line  spectra.  In  responding  to  this  excitation,  the  pilot's  output  power- 

spectral  density  will  also  have,  in  general,  random  and  line  spectra  elements,  but  not  all  of  the  output 
will  be  linearly  correlated  with  the  forcing  function.  Accordingly,  the  remnant  spectral  density, 
constitutes  the  difference  between  the  linearly-correlated  and  total  pilot  output  power-spectral  density. 

In  principle,  the  rernnam;  spectrum  may  also  have  random  and  line-spectra  components;  although  the  spectral 
lines  will  not  appear  at  the  discrete  line  spectra  frequencies,  utt,  present  in  the  forcing  function.  The 
spectral  formulas  in  Table  1  represent  these  statements  mathematically. 

The  pilot's  dynamics  as  a  transfer  element  are  given  by  the  open-loop  describing  function,  Yp,  which 
is  also  represented  along  with  the  controlled  element  transfer  function  in  the  closed-loop  describing 
function,  H.  These  are,  in  general,  obtained  using  cross- spectral  measurements  as  shown  in  the  table. 

Noted  there  is  the  sometimes  used  estimate  for  Yp  of  $>ec/i>ee.  As  can  be  appreciated  from  the  formula, 
this  measurement  can  give  good  results  at  frequencies  where  the  forcing  function  power  and  its  consequences 
are  much  larger  than  the  remnant  power  and  its  consequences.  On  the  other  hand,  when  the  forcing  function 
effects  are  not  dominant,  the  answer  obtained  approaches  —  l/Yc.  Techniques  which  help  circumvent  problems 
caused  with  low  forcing  function  power  levels  are  given  in  Refs.  41-43. 

The  "linear  correlation"  and  signal-to-noise  ratio  are  useful  as  indicators  of  the  nature  of  pilot 
operations  and  as  a  means  for  the  determination  of  remnant.  These  quantities  depend  on  the  effective 
bandwidth,  aq>,  of  the  measuring  apparatus  or  data  processing  procedure.  Clearly,  as  this  becomes  more  nar¬ 
row,  the  value  of  p  increases  regardless  of  the  magnitude  of  the  remnant,  approaching  1  as  uq,  approaches 
zero. 

Finally,  relative  remnant  measures  finish  off  the  Table  4  summary. 

The  principles  and  practices  of  measuring  pilot  characteristics  are  extremely  well  developed  and  have 
an  extensive  literature.  A  representative  cross  section,  which  includes  both  time  domain  and  frequency 
domain  considerations,  is  given  in  Refs.  40,  44-3,4.  There  are  many  pitfalls  for  the  unwary  in  practical 
application  of  these  methods,  with  the  — 1  /Yc  example  noted  above  being  perhaps  the  most  common.  The 
references  cited  offer  a  cross  section  of  others. 

A  measurement  detail  of  considerable  importance  is  the  character  of  the  forcing  function  or  distur¬ 
bance.  When  compensatory  behavior  is  to  be  measured  the  ideal  attributes  include: 

1 .  Random  appearance  so  that  the  operator  cannot  detect  any  internal  coherence  in  the 
forcing  function  and  thereby  adopt  a  higher  level  of  behavior  (unless  the  intent 
is  to  induce  that  higher  level). 

2.  Frequency  content  which  encompasses  the  bandwidth  of  the  manual  control  system 
under  test,  thus  insuring  that  all  significant  modes  of  the  system  are  excited. 

3.  Evoked  outputs  easily  recognized  and  distinctly  different  from  constant-coefficient 
linear,  nonlinear,  constant-rate  sampled,  and  other  idealized  system  fox-ms  which 
may  offer  a  mathematical  basis  for  human  pilot  modeling. 

4.  Precisely  known  statistics,  to  the  extent  possible,  thereby  eliminating  a  source 
of  uncertainty  in  experimental  results. 
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TABLE  t .  MEASUREMENT  SUMW  1Y 


Yp.O, 


p  >  ^nn 


Forcing  System 


Pilots 

Force 


Function  &  pi,n(  Output  Pil0' 

mi  -  m  p'l0's!Try  Fit)  “* 


^eet 


and 

Equalization 

Dynamics 


Manipulator 

and 

Pilot  Load 
Dynamics 


Pilots 

Motion 

Output 

c(t) 

— - - )•» 

$cc  > 


Controlled 

Element 


System 

Output 

m(t) 


SPECTRA 


Forcing  Function:  ■t^t'*-')  = 


Tii(c')  +  n  £)  <ri(<4n)f>(a>  -  <%) 

m=l 


Pilot  Output:  *ee(tt>)  = 


|H|2  «ii 


+  *„„(«>) 


Linearly  Correlated 


Remnant 


DESCRIBIN3  FUNCTIONS 


Closed  Loop: 


M  ' 

Tcc(“>)  +  *  L  «f(<%)&(«>  -  u^)  +  *nn(“>) 

m=1 


Controlled  Element:  Y0  = 


■!'ic(  jo) 
®ii(«>) 

®icLM 

$ie(  du>) 

$im(j«>) 
$ic(  ja<) 


1  +  YpYc 


„  ,  fectj") 
;  Note:  :■ 

'J'ee(a) 

^mrn^)  _  ^cm^j115) 
‘S’mct  ju»)  “  |t'cc(U)) 


Yp»ii(ot)  -  Ygitnn(^) 
$ii(«>)  +  |YC|2  <Jnn(<") 


LINEAR  CORRELATION  AND  SIGNAL/NOISE 


ruk+o>b/2  ~ 

I  lWuV  4o 

^aVn-a-b/2  _ 

«?§(<%)  +  f ‘%+U''b/2  <Pcc(a>)  da, 


“WinW 
«<P§  +  a'b<J>cc(utn) 


S/N  = 


oS(“fa)  .  «<J>1  +  a^ccKi) 

— —  -6  ■  .  »  — .1  

1  -  p2(<%)  <»b«Pnn(«k) 


RELATIVE  REMNANT 


o|c  =  '  -  (Output)  ;  pfg  =  1  -  — -  (Error) 
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5*  Repeatability  to  permit  reproduction  of  meaningful  tests  under  varying  conditions. 

6.  Representation  of  physical  signals  such  as  turbulence,  radar  noise,  evasive  targets, 
with  which  the  pilot  must  contend. 

7-  Pseudo-randomness  so  that  the  forcing  function  will  appear  crgodic  when  its  time 
average  is  taken  over  a  specified  interval. 

8.  A  Gaussian  amplitude  distribution  so  that  Gaussian  input  describing  function  theory 
may  be  used. 

A  forcing  function  comprising  a  properly  selected  sum  of  sinusoids  which  have  frequencies  that  are  integral 
multiples  of  the  run  length  is  perhaps  the  best  approach  to  achieving  these  desirable  features.  In  second 
place  would  be  a  rect.»’;d,  and  hence  repeatable,  noise  source  with  precisely  known  statistics.  Sometimes, 
special  functions  or  disturbances  cannot  be  used,  and  the  experimenter  must  then  make  do  with  whatever 
excitation  sources  are  available,  with  the  concomitant  reduction  in  quality. 

The  third  desirable  feature  listed  above  can  be  especially  valuable  in  providing  insights  as  to  the 
most  appropriate  mathematical  models  to  use  to  describe  the  pilot's  actions.  Table  'j  lists  the  features 
of  come  system  types  which  would  be  revealed  by  measurements  of  the  kind  summarised  in  Table  L  when  sub¬ 
jected  to  a  periodic  input  with  frequency  iq,  and  amplitude  o^.  Careful  examination  of  data  for  inputs 
with  sums  of  sinusoids  in  the  context  of  the  differences  indicated  in  this  table  can  go  a  long  way  in 
selecting  the  best  paradigm.  Many  cf  these  insights  are  lost  if  a  forcing  function  with  a  continuous 
power-spectral  density  is  used.  In  this  case,  the  remnant  for  all  but  the  linear  constant-coefficient 
system  will  also  be  a  continuous  power- spectral  density. 


TABLE  5 


PROPERTIES  OF  VARIOUS  SYSTEMS 


TYPE  OF  SYSTEM 

DESCRIBING 

FUNCTION 

REMNANT 

P2 

pfc  or  pae 

Linear,  Constant  Coefficient 

*P 

o 

1 

1 

Linear,  Sampled  Data 
(Constant  Frequency) 

YP 

Line  Spectra 

2 

kOj&(ub  ±  BU's) 

1 

<  1 

Linear,  Random  Time  Variations 

*P 

Continuous  Spectra 

kof<e(«>) 

<  1 

<  1 

Nonlinear  Constant  Coefficient 

Yp(0i) 

Line  Spectra 

kf(  o^BOncq,) 

1 

<  1 

System  with  Noise  Injection 

Line  or  Continuous 

«nn  *  f(°i) 

<  1 

<  1 

PRECEDING  PAGE  BLANK 
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TABLE  SUMMARY  OF  SINGLE-LOOP  COMPENSATORY  SYSTEM  DESCRIBING  FUNCTION  PLUS  R04NANT  DATA 


INVESTIGATOR 
( Reference ) 

NUMBER  OF 
RANDOM- 
APPEARING 
PORC.  FUNC. 
TYPES  INVEST. 

NUMBER  OF 
CONTROLLED 
ELEMENT 
FORMS 

INVESTIGATED 

MANIPULATOR 

CONTROLLED  ELEMENT  DYNAMICS; 

REMARKS 

Tustin 

(7) 

P 

Spade  grip, 
spring  restraint 

Simulated  tank  gun  turret  tracking;  both 

Yc  i  Kc's. 

Russell 

(PO) 

P 

r, 

Handwheel,  no 
restraint 

7c  =  Kc,  Kc  /s,  Ke( s + z )  /( s + p) . 

Goodyear 
(  ?  2,  lU 

p 

Aircraft  control 
stick 

Simulated  aircraft  pitch  attitude  control 
in  both  stationary  and  pitching  simulator. 

Krendel,  et  nl . 

}  o ) 

■» 

7 

Aircraft  control 
stick 

Simulated  control  of  aircraft  lateral  and 
longitudinal  axes  in  tail-chase,  with  and 
without  airframe  dynamics:  two- 
dimensional  input. 

Elkind 

(  r  '  ) 

OQ 

Pencil- like 
stylus ,  no 
restraint 

Yc  =  Kc;  single-dimensional  input;  some 
remnant  data. 

Seckel,  et  al. 

(!■’) 

1 

1 

Aircraft  control 
wheel 

Attitude  control  of  aircraft  lateral  and 
longitudinal  axes  in  both  flight  and 
fixed-base  simulator;  two-dimensional 
input . 

Hall 

f  n,  iA,  ) 

j 

c  . 

Aircraft  control 
wheel 

Simulated  aircraft  pitch  attitude  control 
while  also  controlling  a  fixed  set  of 
lateral  characteristics;  many  Yc  para¬ 
meter  variations:  two-dimensional  input. 

McRuer,  et  al. 

■ 

*  . 

5 

Lateral  side 
stick 

Control  of  a  wide  range  of  idealized  dyna¬ 
mics  contrived  to  e”oke  n  complete  range 
of  operator  transfer  characteristics; 
many  Yc  parameter  variations:  good 
remnant  data. 

Stapleford, 
et  al.  . 

(•  ~)  ■ 

.  .1 

1 

1 

Lateral  side 
stick 

Simulated  aircraft  lateral  control  in  a 
tail-chase;  two-dimensional  input. 

Magdalene 
and  McRuer 
( 1  ® ,  ‘  0 ) 

P% 

Longitudinal  side 
stick;  various 
restraints 

t  r\ 

Yc  =  Kc,  Ke  s.  Kc  s'  ;  extreme  ranges  of 
manipulator  restraints. 

Jex  and 
McDonnell 

f'O.  '-‘V) 

r- 

Rigid  side  stick 

First-  and  second- order  suberiticnl 
tasks.  Good  remnant  data. 

Smith,  et  al. 

tz)  - 

? 

Aircraft  control 
stick 

Lateral  bank  angle  control;  simultaneous 
longitudinal  stabilization;  both  flight 
and  fixed-base  simulator;  single¬ 
dimensional  input. 

Tie  well,  et  'al. 

{ <> .  o») 

c 

Aircraft  control 
stick 

Lateral  bank  angle  control;  simultaneous 
longitudinal  stabilization;  both  flight 
and  fixed-base  simulator;  single¬ 
dimensional  input. 

Levison  and 
Elkind 

.(*7) 

7 

7 

Spring  restrained 
side  stick 

Yc  =  Kc ,  Ka 's,  and  Kc  V ;  single-axis 
base  runs  for  two-axis  experiments;  good 
remnant  data. 

Shirley 

!  V) 

u 

Spring  restrained 
stick 

Simulated  lateral  aircraft  control;  many 

Yc  parameter  variations;  complete  remnant 
data . 

Allen  end  Jex 

^,7) 

2 

I.aternl  side 
stick 

Yc  -  Kc  ,s  and  Kc's^;  compensatory  base 
runs  for  pursuit 'compensatory  comparison. 

Gordon-Smith 

trp.) 

7 

- 

Free-moving  and 
pressure  center 
sticks 

Yc  =  Kc  >'s  dynamics;  good  high-frequency 
!  neuromuscular  actuation  system  data. 

Staple  ford , 
et  nl. 

(m9) 

i 

> 

Aircraft  control 
stick 

Simulated  aircraft  pitch  attitude  con¬ 
trol;  base  runs  for  multiloop  single-loop 
comparison. 

Magdnleno 
and  McRuer 

(7’) 

i 

. 

5 

Lateral  side 
stick:  rudder 
pedals 

First-,  second-,  and  third-order  sub- 
critical  tasks;  good  neuromuscular 
actuation  system  data. 

CENTRAL  ELEMENTS 


NEUROMUSCULAR  SYSTEM 


& 


Figure  15-  Model  for  Human  Pilot  Dynamics  in  Single-Loop  Compensatory 
System  with  Random-Appearing  Inputs 


neuromuscular  system  shown  in  Fig.  If)  has  a 
feedback  element  labeled  as  spindle  and  tendon 
ensembles,  although  spindle  characteristics 
may  well  be  predominant  for  the  small  motion 
and  relatively  light  forces  involved  in  most 
measurements  thus  far  accomplished.  The  effec¬ 
tive  dynamics  of  the  closed-loop  neuromuscular 
system,  from  the  alpha  motor  neuron  command 
signals  to  mani;  ulator  force,  can  be  approxi¬ 
mated  over  a  wide  frequency  range  by  the  third- 
order  transfer  function  shown.  This  form  is 
also  compatible  with  small  perturbation  dyna¬ 
mics  based  on  analytical  models  of  muscle  and 
manipuYator  characteristics  (Refs,  i'8,  71, 

77-7r  )  •  The  parameter  values  are  strongly 
dependent  upon  the  steady-state  neuromuscular 
tension,  y0,  due  to  the  gamma  mot^r  system. 

The  gamma  system  also  affects  the  dynamics  of 
the  spindle  ensembles.  This  is  ceplcted  by 
the  arrows  indicating  variation  n  the  .’.Sp 
and  Psp  factors  in  the  neuromuscular  system 
feedback  block. 


Figure  1 6.  Functional  Diagram  of  Elemental 
Agonist/Antagonist  Neuromuscular  System 
Elements  Involved  in  Tracking 


This  degree  of  model  complexity  is  a  minimum  to  have  even  marginal  value  in  physiological  descriptions. 
It  is  also  often  needed  for  the  study  of  limb/manipulator  system  dynamics  in  aircraft  control,  as  in 
determining  the  effects  of  bobweights  and  primary  control  system  hysteresis  on  pilot-induced  oscillations. 
For  other  aspects  of  pilot/aircraft  analysis  the  neuromuscular  dynamics  are  so  high  in  frequency  as  to  be 
relatively  unimportant  in  their  details.  For  these  cases  a  pare  time  delay,  x jjj.j,  or  a  first-order  lag  can 
be  used  as  a  low-frequency  approximation  (at  frequencies  lower  than  any  or  the  neuromuscular  system  break¬ 
points).  This  delay  will  be  given  by: 


•  1  1 
=  tNi  +  —  +  7^-7— 
N  Fsp  *sp 


(~  ) 


Turn  now  to  the  central  and  input  elements  at  ft  of  Fig.  1;, .  As  shown  there,  the  pilot  can 

develop  a  neuromuscular  system  input  command  which  is  the  summation  of  a  lag,  proportional,  lead,  or 
double-lead  function  of  the  system  error.  The  integral  and  proportional  channels  have  a  basic  time  delay, 
tc,  associated  with  them.  The  higher  derivative  channels  have  additional  incremental  delays.  These  incre¬ 
mental  time  delays  constitute  the  dynamic  cost  of  pilot  lead  generation.  They  are  about  1  /;.  sec  for  rate 
(tp)  and  greater  than  l/?  sec  for  the  acceleration  channel  (t^). 

It  is  primarily  because  of  the  latency  differences  for  proportional,  rate,  and  acceleration  low-frequency 
lead  equalization  that  these  are  shown  as  separate  parallel  channels.  The  independence  schematised  is  over¬ 
simplified,  for  common  neurological  apparatus  is  undoubtedly  present  for  each  function.  These  are  modeled 
here  by  the  common  sensory  pathway  block  following  the  system  error  and  the  central  processing  and  computa¬ 
tional  block  following  the  four  parallel  channels.  A  slight  difference,  not  shown,  in  integral  and  pro¬ 
portional  channel  effective  latencies  is  indicated  by  existing  data  but  conceivably  these  channels  could 
be  combined  into  a  single -channel  adjustable  lag-lead  element. 

Besides  the  different  effective  time  delays,  the  other  evidence  for  parallel  channels  is  the  difference 
in  respons.  quality  as  a  function  of  low-frequency  equalization  supplied  by  the  pilot.  When  proportional, 
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integral,  or  a  lag-lead  combination  is  exhibited, 
tile  pilot's  output  amplitude  distributions  are 
Gaussian  when  the  input  amplitude  is  Gaussian. 

On  tne  other  hand,  when  very-low-frequency  leads 
are  present,  as  if  operations  were  through  the 
rate  or  acceleration  channels,  the  pilot's  output 
distirbutions  are  distinctly  non-Gaussian.  Typi¬ 
cally,  these  are  bimodal,  as  illustrated  in 
;.  17-  Time  traces  tend  to  appear  more  dis¬ 
crete  and  pulselike  as  well.  This  characteristic 
difference  is  also  reflected  in  changes  in  the 
remnant . 

The  association  of  incremental  latencies  with 
differences  in  the  neural  pathways  and  apparatus 
takes  us  back  a  century  or  more  to  the  very  earli¬ 
est  work  in  the  reaction  time  studies  of  Helmholtz, 
Wundt,  and  others  (Ref.  ^6) •  They  examined  a  large 
number  of  situations  which  exhibited  different 
reaction  times  and  tried  to  assign  these  differ¬ 
ences  to  various  component  functions.  In  their 
day,  their  attributions  often  turned  out  to  be 
wrong,  and  we  can  only  hope  that  those  indicated 
here  do  not  repeat  that  history! 

The  channel  gains  and  the  time  constant,  Tj, 
are  all  shown  os  variable  quantities.  To  a  first 
approximation,  they  are  adjusted  such  that  the 
crossover  model  applies.  To  a  higher  order  of 
approximation  the  adjustments  obey  the  Analytical/ 
Verbal  model.  Both  models  are  discussed  farther 
below. 


Figure  17 •  Example  of  Bimodal  Amplitude 
Distribution  for  Pilot  Output 
for  Yc  =  kc/s2 


C .  CROSSOVER  MODEL  TOR  SHKJtE-LOOP  SYSTEMS 


Consider  the  crossover  model  first.  Typical  data  for  velocity,  first-order  subcrit'.eal  divergence, 
and  acceleration  controlled  elements  are  shown  in  Fig.  l8  with  ±1o  bands.  These  data  rre  more  or  less 
typical  of  any  of  the  rate  and  acceleration  controlled  element  data  listed  in  Table  6.  Further,  the 
Icuc/ja'ds  amplitude  ratio  trend  in  the  region  of  crossover  is  characteristically  present  1n.  ,i>uch  of  the 
Table  6  data  with  more  complex  controlled,  elements.  In  fact,  for  those  duta  sources  from  the  last  decade, 
there  exists  a  remarkable  consistency.  In  most  cases,  investigators  have  taken  pains  to  show  tie-ins  with 
earlier  results  from  other  investigators.  For  instance,  typical  comparisons  made  in  Ref.  66  for  the  con¬ 
trolled  elements  Y0  =  Kc,  Kc/s,  and  Kq/s2,  which  cover  three  widely  varying  modes  of  control  behavior  from 
the  pilot,  are  reproduced  in  Fig.  I9a-c.  Other  comparison  data  are  shown,  for  example,  in  Refs.  56,  60, 
and  62-69.  For  a  given  controlled  element  such  differences  as  do  exist  between  data  sets  are  ordinarily 
attributable  to  differences  in  forcing  function  or  manipulator  task  variables.  In  other  cases,  the 
experimental  conditions  arc  such  as  to  permit  the  differences  to  be  ascribed  directly  to  the  effects 
environmental  or  operator-centered  variables. 


jve  *  475  rod/sec ,  re  1 018  sec ,  a  =  Ollrod/sec]  uc  =  4  9  rod /sec ,  r,  =  013  sec  ,  a  >0!9rad/sec]  uc  *325rcd/$ec  ,  r,*033sec  033 rod/sec] 


(a)  (b)  (c) 

Figure  18.  Data,  Simple  Crossover  Models,  and  a  Crossover  Models  for  Elementary 
Controlled  Elements;  =  2.5  rad/sec  (Ref.  56) 
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O  Ref.56  ,  o  Rtf.  21  ,  |  Ref. 66  (mean  25  runs  icr ) 


(a)  Yc  ■  I 


(b)  Yc  =  I/s 


(c)  Yc  1  I/s* 


Figure  19-  Comparisons  of  Typical  Data  from  Different  Investigators  (from  Ref.  66) 


Visual  examination  of  the  data  shown  in  Figs.  18  and  19  indicates  that,  to  a  first  order  at  least,  the 
simple  crossover  model  developed  in  Chapter  II  from  "eyeball"  fitting  of  time  traces  is  also  applicable 
when  the  more  sophisticated  spectral  analysis  procedures  ar''  ised.  The  model  is  apparently  a  better 
description  of  amplitude  ratio  than  of  phase  characteristics.  Also,  the  amplitude  ratio  for  Yc  =  Kc/( ja - 2), 
in  Fig.  18b,  tends  to  be  somewhat  flatter  than  -20  dB/decade. 

For  the  Table  6  data  in  general  the  crossover  model  is  neither  appropriate  nor  accurate  at  frequencies 
much  less  than  or  much  greater  than  the  crossover  frequencies.  Nonetheless,  it  usually  leads  to  essentially 
correct  closed-loop  dynamic  characteristics  because  the  actual  shape  of  the  open-loop  describing  function 
far  from  the  crossover  region  usually  has  little  effect  on  the  dominant  closed-loop  dynamics  (see  Eqs.  it, 

15,  and  21). 


Although  the  basic  YpYc  approximate  form  is  the  same  for  the  Figs.  18  and  19  data,  and  for  most  of  the 
cases  listed  in  Table  6  as  well,  the  numerical  values  of  crossover  frequency  and  effective  time  delay  are 
not.  These  are  both  functions  of  the  task  variables.  For  low-pass  forcing  functions,  idealized  as  rectan¬ 
gular  power-spectral  densities  with  bandwidth  and  rms  value  o*,  the  crossover  frequency  and  effective 
time  delay  can  be  represented  (Ref.  56)  as: 


u-c  =  <i'c0(Yc)  +  Axc(a-i)  ,  Ai.c(0)  -  0 
Te  i  t0(Yc)  -  Are(ai)  ,  Are(0)  =  0 


(26) 


As  the  forcing  function  bandwidth  is  reduced  the  phase  margin  also  decreases,  becoming  approximately  zero 
for  <n  =  0.  Although  this  trend  towards  zero  phase  margin  can  only  be  demonstrated  directly  for  <i'i  /  0,  it 
is  a  common  observation  that  signals  circulate  throughout  the  manual  control  loop  without  any  forcing  func¬ 
tion  as  long  as  the  operator  is  in  active  control.  In  the  absence  of  other  inputs  or  disturbances,  the 
p.-  sence  of  these  signals  implies  an  on-the-average  condition  of  zero  phase  margin.  So,  both  the  direct 
ar.  I  indirect  evidence  indicates  that  phase  margin  is  zero  when  is  zero.  Then,  the  neutrally  stable 
crossover  frequency,  u<;0,  and  the  basic  time  delay,  r0,  are  related  by: 


.  it 

^0"-Cq  =  2 


(i7) 


To  a  first  approximation,  is  nearly  zero  (i.e.,  Ao^  =  0)  for  controlled  element  forms  encountered 

in  aircraft  (although  it  is  not  quite  true  for  Yc  =  Kc).  If  this  approximation  is  taken  to  be  exact  for 
simplicity,  then  the  principal  variations  in  the  crossover  model  are  those  of  t0  with  controlled  element 
characteristics  and  Are  with  forcing  function  bandwidth. 

In  general,  t0  consists  of  effective  transport  lags  due  to  both  the  pilot  and  the  controlled  element. 
The  controlled  element  contribution  will  be  approximately  the  difference  between  the  controlled  element 
high-frequency  (much  greater  than  crossover)  lags  and  leads.  With  this  assumed  negligible  or  otherwise 
accounted  for,  the  net  to  depends  on  the  amount  of  lead  generation  required  of  the  pilot  to  establish  the 
approximate  -20  dB/decade  slope  in  the  region  of  crossover.  When  the  lead  required  has  been  determined, 
the  latency  can  be  estimated  from  Fig.  20a.  The  curve  shown  there  represents  a  grand  average  from  a  very 
large  number  of  experiments,  essentially  all  of  those  involving  crossover  region  data  and  spring-restrained 
manipulators  listed  in  Table  6.  Zero  lead  units  is  a  pure-gain  proportional  control;  plus-one  units  is 
first-order  lead  or  rate  control;  plus-two  units  second-order  lead,  etc.  The  time  delays  shown  include 
the  effective  dynamics  of  the  neuromuscular  system  as  an  additive  component  to  the  central  latencies.  The 
total  effective  low-frequency  time  delay  ranges  from  about  0.3  sec  for  the  pure  integral  control  case 
(-1  lead  units)  through  0.8  sec  for  the  acceleration  control.  The  time  delays  are  plotted  as  inverse 
functions  because  in  this  form  they  are  normally  distributed  and  the  variability  of  the  inverse  time  delay 
is  a  constant. 

Figure  20b  ?hows  the  remaining  variation  of  Eq.  26,  that  of  the  incremental  latency.  Ate,  with  ap . 

This  result  is  based  primarily  on  data  taken  by  several  investigators  with  the  so-called  STI  forcing 
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Forcing  Function  Bandwidth,  uq  (rod/sec) 


Figure  20a.  Variation  of  Crossover  Model 
Dynamic  Stimulus-Response  Latency  with 
Degree  of  Pilot  Lead  Equalization 


Figure  £0b.  Variation  of  Crossover  Model 
Incremental  Latency  with  Forcing 
Fur:tion  Bandwidth 
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function  spectrum  made  up  of  10  sinusoids  (Ref-  iC) .  With  other  forcing  function  types  a  somewhat  more 
complicated  adjustment  is  made-  This  will  be  discussed  in  connection  i  ith  the  "precision  model"  presented 
in  the  next  section.  With  the  approximations  made  thus  far,  the  Atc  variation  in  Fig.  20b  is  proportional 
to  tne  phase  margin  variation  with  ^  since  fyj  1  cq^ATe*  This  reduction  in  effective  time  delay  as  for¬ 
cing  function  bandwidth  is  increased  is  due  primarily  to  an  increase  in  steady-state  neuromuscular  tension, 
70,  (shown  entering  the  neuromuscular  block  in  Fig.  1‘).  Thus,  the  decreased  delay  is  predominantly  a 
change  in  the  neuromuscular  dynamics  which  would  be  reflected  in  the  low-frequency  approximation,  tjj.j,  to 
those  dynamics. 

The  foregoing  has  related  tne  parameters  of  the  crossover  model  to  the  task  variables,  primarily 
and  Yc .  The  numerical  values  were  obtained  under  test  conditions  —  specified  in  terms  of  operator- 
centered,  0,  procedural,  it,  and  environmental,  e,  variables  as  discussed  in  Chapter  }  —  so  as  to  insure 
that  Yp  determinations  would  be  useful  for  system  analysis  and  design.  The  operators,  pilots  or  well- 
coordinated  simulator  pilots,  were  highly  motivated,  well-trained  males  paying  full  attention  to  the  con¬ 
trol  tasks.  The  procedural  variables  were  selected  to  enhance  generalizabiiity.  Hostile  or  disruptive 
environmental  vaslables  were  avoided.  Consequently,  the  foregoing  data  are  representative  of  the  best 
high-dynamic  performance,  low- /ariability  parameter  values  which  pilot  can  be  expected  to  generate.  For 
usual,  flight-encountered  situations  the  values  of  uc  are  high,  of  Te  low. 

Ample  evidence  exists  for  the  dependence  of  crossover  model  parameters  on  0  variables.  An  increase  in 
gain  and  decrease  in  ie  occurs  with  training  and  has  been  noted  by  several  investigators  —  for  example, 
Ref.  86.  The  pilot's  performance  criteria  affect  measurements  of  Yp  by  determining  where  he  can  afford  to 
relax  control  and  where  he  must  conform  strictly  to  the  system's  demands.  Thus,  for  example  Ref.  y6, 
large  run- to-run  variability  is  exhibited  by  a  pilot  in  Yp  in  the  low-frequency  range  for  control  of  Kq/s 
jr  Kq/ss;  however,  this  variability  becomes  confined  to  a  very  tight  band  in  and  about  c.q  where  good  per¬ 
formance  demands  become  contraining.  Similarly,  the  variability  expressed  by  the  same  pilot  is  practically 
negligible  when  Yc  is  an  unforgiving  divergence.  Interestingly,  pilots  exhibit  individualized  consistent 
departures  from  the  crossover  model  at  low  frequencies  where  performance  and  stability  are  scarcely 
affected.  Such  behaviir  is  in  the  nature  of  personal  "style,"  and  in  effect  it  represents  decisions  on 
where  to  apply  trim  control. 

When  it  is  desired  to  measure  the  effects  of  a  or  c  variables  upon  Yp,  the  pilot  must  be  constrained 
so  that  his  dynamics  are  tightly  related  to  the  selected  variables.  A  highly  effective  manner  for  accom¬ 
plishing  this  is  by  employing  controlled  elements  of  the  form  Yc  =  l/sMs-X),  k  =  0,  1,  2.  Since  this 
divergent  controlled  element  tightly  constrains  the  allowable  pilot  equalization  near  the  region  of  gain 
crossover,  these  critical  tasks  leave  the  pilot’s  effective  time  delay,  xe,  as  the  sole  determinant  of 
system  stability.  When  the  divergence  is  gradually  increased  until  control  is  lost,  this  "critical" 
divergence  time  constant  is  a  measure  of  re  (Ref.  6l ). 

These  tasks  provide  a  sensitive  method  for  reflecting  u  or  e  variables  onto  a  restricted  measure  of 
the  pilot's  dynamic  performance.  This  technique  is  being  exploited  in  several  laboratories  (a  representa¬ 
tive  group  can  be  found  by  scanning  Refs.  2,  26-28)  to  study  a  variety  of  such  c  variables  as  heat  and 
noise,  and  0  variables  such  as  induced  by  amphetamine,  alcohol,  and  other  drugs;  task-induced  stresses, 
fatigue,  enforced  bedrest,  etc.;  and  even  certain  task  variables  related  to  displays.  The  application  of 
the  critical  task  to  measure  pilot  workload  and  to  aid  in  understanding  pilot  opinion  ratings  will  be 
presented  later  in  this  chapter. 

The  simplified  crossover  model  discussed  above  is  deficient  in  two  major  respects.  First,  for  the 
controlled  elements  with  non-zero  poles,  the  open-loop  describing  function  in  the  region  of  crossover  tends 
to  be  somewhat  less  than  -20  dB/decade.  For  these  cases  the  data  are  better  fitted  with  open-loop  describ¬ 
ing  functions  which  contain  the  controlled  element  dynamics  explicitly,  as  illustrated,  for  example,  in 
Fig.  18b.  Second,  the  phase  at  low  frequencies  is  not  described  too  well  by  the  simple  delay.  These  plase 
lags  at  low  frequencies  usually  do  not  substantially  affect  the  closed-loop  characteristics  because  magni¬ 
tude  |YpY0|  is  much  greater  than  unity  at  the  frequencies  where  the  lags  are  present.  However,  the 
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ubiquitous  nature  of  the  low-frequency  phase  lags,  and  their  occasional  importance  on  closed-loop  dynamics, 
demands  attention. 

This  can  be  accomplished  simply  by  adding  a  catch-all  increment  to  the  low-frequency  phase  which  takes 
into  account  dynamics  having  amplitude  ratio  breakpoints  below,  or  in  the  lowest-frequency  portion  of,  the 
measurement  bandwidth.  The  approximation  derives  from  considering  that  the  low-frequency  phase  lag  is  due 
to  an  equal  number  of  logs  and  leads  having  breakpoints  ordinarily  below  the  lower  frequencies  at  which 
the  crossover  model  applies.  Thus,  for  M  leads  and  lags  occurring  at  l/Tpea^  and  l/Tiagi,  the  phase  will 
be: 

M  -l  M 

A<*low  =  S  tan~  Tleadi“  -  £  tan~  Tlagio>  (28) 


When  all  of  the  breakpoints  occur  below  the  measurement  bandwidth  or  otherwise  outside  the  bandwidth  of 
crossover  model  validity,  the  pi.  se  angle  within  the  measurement  or  validity  bandwidth  will  be  approximately: 


A?low  = 


where 


a 


The  effective  time  constant,  l/a,  describes  the  effect,  within  the  crossover  region,  of  leads  and  lags  below 
that  frequency  band.  In  this  ’sense  it  is  analogous  to  Te,  which  lumps  high-frequency  phenomena  into  a 
simple  low-frequency  approximation  suitable  within  the  measurement  bandwidth.  In  the  describing  function 
the  low-frequency  term  is  represented  as  e'i0- A». 

The  open-loop  describing  functions  of  Fig.  18  are  curve-fitted  using  both  the  simple  and  the  a  cross¬ 
over  models,  which  differ  only  in  whether  or  not  they  include  the  low-frequency  phase  correction.  The  range 
of  validity  for  the  a  component  is  obviously  less  than  the  lover  frequency  at  which  the  phase  margin  will 
be  zero. 

The  relationship  of  the  crossover  model  and  Fig.  15  for  simple  controlled  elements  is  summarized  in 
Table  7.  For  more  complex  controlled  elements  the  channels  are  appropriately  combined  to  achieve  the 
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TAB  IS  7 

CONNECTION  BETWEEN  COMPLETE  AND  CROSSOVER  PILOT  MODELS 


CONTROLLED  ELEMENT 
CHAHACXciftioX  iC  b 

PRIMARY 
PILOT  CHANNEL 
ACTIVATED 

CROSSOVER  REGION 
PILOT  AMPLITUDE 
RATIO 

EFFECTIVE  PILOT 
TUGS  DELAY 

Kc 

Integral 

Kl 

Tidw 

m 

Ko 

s 

Proportional 

Kp 

Tc  +  TJJM 

Kc 

S2 

Rate 

tR  +  Tc  +  tNM 

Kc 

S2 

Acceleration 

KA(d«>)2 
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crossover  model  properties.  Thus,  when  a  lag-lead  equalization  is  needed,  the  integral  and  proportional 
channels  can  both  be  used  to  give  a  net  equalization  given  by: 


Kj 

t?TT7  +  kp 


1 

(Kp  +  Kj) 

— ^2 TrS  +  1 

Ko  +  Kj  1 

( 

l’js  +  1) 
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Similarly,  when  a  lead  characteristic  is  required,  the  proportional  and  rate  channels  can  be  used  together. 
Remarkably,  there  is  very  little  evidence  •  hat  more  than  two  channels  are  used  simultaneously  to  create 
complex  equalization.  This  observation  ma'  well  be  an  artifact  of  the  particular  experiments  which  have 
thus  far  been  accomplished,  for  in  almost  a.l  of  these  the  crossover  model  characteristic  can  be  achieved 
with  no  more  than  two  channels. 

1).  PRECISION  MODEL  TOR  SINGLE-LOOP  SYSTEMS 

In  order  to  cover  a  broader  frequency  range  than  permitted  by  the  crossover  model,  a  so-called  precision 
model  is  used.  The  pilot's  operation  of  his  acceleration  channel  is  extremely  difficult,  cannot  long  be 
maintained,  and  is  associated  with  extremely  unfavorable  subjective  ratings  and  excessive  workload.  Because 
of  these  features,  this  channel  is  probably  never  useful  on  aircraft-like  controlled  elements.  Consequently, 
the  precision  model  need  not  take  it  into  account.  With  this  proviso,  the  general  describing  function  form 
which  describes  the  transfer  characteristics  of  the  human  pilot  in  single-loop  compensatory  situations  for 
all  of  the  Table  6  data  base  is  given  by: 


VERY-LOW- 

PURfi  TIME  SERIES  FREQUENCY  NEUROMUSCULAR 

GAIN  DELAY  EQUALIZATION  LAG-LEAD  ACTUATION  SYSTEM 


where: 


Tk  *k 


(32) 


and: 


Tn  i 


tN, 


25w 


The  neuromuscular  actuation  portion  presumes  the  neglect  of  the  very-high-frequency  lead-lag  of  the  spindle/ 
tendon  organ  ensembles  implied  in  fig.  15.  It  can,  of  course,  be  approximated  at  low  frequencies  by  the 
first-order  neuromuscular  lag  or,  at  even  lower  frequencies  still,  by  the  pure  time  delay.  The  very-low- 
frequency  lag-lead  is  occasionally  indicated  in  very-low-frequency  data  by  at  least  the  low-frequency  lead 
breakpoint,  l/TK-  It  can  also  be  approximated  by  the  single  parameter,  a.  There  is  evidence  that  a  and 
*N  arc  covarying  quantities  (Ref.  77)  with  forcing  function  bandwidth  and/or  neuromuscular  average  tension. 
Therefore,  the  effect  represented  by  the  very-low-frequency  lag-lead  may  be  neuromuscular  in  origin;  although 
no  physiological  connections  have  yet  been  established. 


The  major  mid-frequency  action  elements  of  the  precision  model  are  the  gain,  Kp,  latency,  t,  and  the 
adjustable  lead-lag  or  lag-lead  represented  by  the  ratio  (T^ja'+  1  )/(Tjjo>+  1 ).  These,  of  course,  derive 
from  various  approximate  summations  of  the  "central  elements"  channels  of  Fig.  15. 

Commonly  used  simplifications  for  the  precision  model  take  advantage  of  t-he  .very-low-  and  high-frequency 
approximations  indicated.  For  conditionally  stable  systems  the  low-frequency  phase  can  be  an  important 
feature  of  the  manual  control  system,  and  an  appropriate  simplified  version  of  Eq.  52  is  then: 


y  c.  K_  1  J  \  e~ J[w(t  +TM)  +0/0)} 

P  “  \Tpjo)  +  1/  e 


(33) 


This  form  is  also  adequate  for  most  other  systems  as  well.  For  systems  wherein  low-frequency  performance  is 
essentially  unaffected  by  the  low-frequency  phase  lag  term,  e~Ja/U),  Eq.  55  can  be  simplified  to: 
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In  either  of  Eqs.  33  or  the  e"  is  interchangeable  with  (Tfjjuj  +  1 )~' .  Furthermore,  if  oil  <  1. 
e“J“®  i  -[(T/2)ju>  —  1  }  /  [(T/2)ju<  +  1  ]. 

Essentially  all  the  components  of  the  precision  model  and  its  simplified  versions  are  adjusted  as 
functions  of  the  controlled  element  dynamics  and  forcing  function  spectra,  the  latter  can  have  many  shapes 
and  sizes,  yet  we  have  thus  far  emphasized  spectra  which  can  adequately  be  defined  by  a  bandwidth,  u.^,  and 
a  forcing  function  rms  amplitude,  .  Those  are  unambiguous  quantities  only  when  the  forcing  function  has 
a  rectangular  spectrum,  achieved  by  using  many  sinusoids  as  in  Ref.  21,  or  an  essentially  rectangular  spec¬ 
trum  plus  a  very  limited  number  of  extremely  low-amplitude  sinusoids  at  higher  frequencies,  as  with  tne  STI 
forcing  function.  For  other  spectral  shapes  the  bandwidth  can  be  defined  in  several  ways.  The  best  metric 
yet  found  (Ref.  35)  to  put  different  shapes  on  a  comparable  effective  bandwidth  basis  is: 
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This  definition  for  effective  forcing  function  bandwidth  works  very  well  as  a  means  tc  consolidate  data  for 
the  20  forcing  function  spectral  shapes  considered  in  Ref.  33 .  It  also  reduces  to  u.-^  for  rectangular 
specti'a . 

To  specialize  the  precision  model  to  a  particular  controlled  element  and  forcing  function  combination, 
the  adjustment  rules  are  applied.  Those  given  below  are  similar  to  the  Ref.  72  set,  modified  by  conclusions 
based  on  data  of  Refs.  68  and  78  and  the  analysis  of  Ref.  79. 

Equalization  selection  and  adjustment.  A  particular  equalization  is  selected  from  the 
general  form  K(TjJa)+  I  )/(Tjj(u  +  1 )  such  that  the  following  properties  obtain: 

(a)  The  system  can  be  stabilized  by  proper  selection  of  gain,  preferably  over 
a  very  broad  region. 

(b)  Over  a  considerable  frequency  range  in  the  crossover  region  (that  frequency 
band  centered  on  the  crossover  frequency,  ay.)  |TpYe!dB  nas  approximately  a 
—20  dB/decade  slope. 

(c)  | YpY0|  »  1  at  low  frequencies  to  provide  good  low-frequency  closed-loop 
response  to  system  forcing  functions  (commands). 

Examples  of  form  selection  and  basic  adjustment  are  provided  in  Table  8. 


TABLE  8 

TYPICAL  PILOT  EQUALIZATION  CHARACTERISTICS 


CONTROLLED  ELEMENT 
APPROXIMATE  TRANSFER 
FUNCTION  IN  CROSSOVER 
REGION 

PILOT 

EQUALIZER  FORM 

EQUALIZER  ADJUSTMENTS  j 

LOW-FREQUENCY 
(a)  «  0>c) 

MID-FREQUENCY 
(too  Region) 

HIGH-FREQUENCY 
(-  *  a^) 

Lag-Lead 

1 

Tl 

— 

Tj,  to  partially 
offset  t  +  Tjj 

Juj 

High-Frequency 

Lead 

— 

— 

T|,  to  partially 
offset  t  ■*  Tfj 
(Tjp>c  <  l) 

Kc 

(jo)2 

Low-Frequency 

Lead 

1 

tl 

— 

Tl  not  available 
to  offset  t  +% 

Kc 

Mid-Frequency 
Lead  (T  >  1 ) 

— 

TL  =  T 

— 

ju>(Tju>  +  1 ) 

High-Frequency 
Lead  (T  <  t ) 

— 

— 

Tl  to  partially 
offset  ~  +T(j  +T 

Kc 

Low-Frequency 
Lead  con  «  “ 

1 

tl 

— 

— 

(S)  +  i ju>  + 1 

Lag-Lead 

1 

Ti 

— 

Tl  to  partially 
offset  t  +Tjj 
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Effective  time  delay.  After  the  appropriate  equalization  form  has  been  adopted,  the  net 
effect  in  the  region  of  crossover  of  high-frequency  (relative  to  crossover)  leads  and 
lags  can  be  approximated  by  replacing  these  terms  in  Eqs.  32-31*  with  a  pure  time  delay 
term,  e_^urTe.  The  effective  time  delay,  te,  is  the  sum  of  all  the  human  pilot's  pure 
time  delays  and  high-frequency  lags  less  the  high-frequency  leads,  i.e.. 


Te 


t+TNi  +^-TLhi 

•t  •*  TfJ  -  1Lhi 


(for  the  crossover  model) 
( fo.B  Eqs.  33  and  34) 


The  notation  Ti,j,i  means  that  only  those  T^'s  used  to  partially  compensate  for  high- 
frequency  phase  lags  (e.g.,  see  Table  8)  are  involved;  otherwise,  T^hi  =  0.  In  general, 
xe  depends  on  both  the  controlled  element  dynamics  and  the  forcing  function  bandwidth. 
These  dependencies  are  approximately  serial,  viz., 


Te(Yc,oi)  =  t0(Yc)  -ATe(u>ie)  where  Aie(0)  =  0 


(a)  Estimation  of  To*  t0  can  be  estimated  from  the  effective  order  of  Yc  in 
the  crossover  region  using  Fig.  20a. 

(b)  Incremental  Te  due  to  forcing  function.  Te  is  essentially  equal  to  t0 
when  the  forcing  function  bandwidth,  u is  zero  or  very  small.  As  u;^  is 
increased  the  neuromuscular  lag,  Tjj,  and/or  the  equalizer  lead,  TLhi,  are 
adjusted  to  reduce  the  net  value  of  te.  A  first-order  approximation  for 
this  effect,  good  for  all  controlled  elements,  is  given  in  Fig.  20b  or  by: 


Are  i  O.O0u>ie 


(36) 


where  Are  is  in  seconds  and  <uie  is  in  radians/second. 

Crossover  frequency,  tug. 

(a)  Rectangular  and  quasi-rectangular  forcing  function  spectra  (discrete 
power-spectral  densities  which  are  essentially  rectangular  and  low-pass 
continuous  spectra  with  a  high-frequency  cutoff  equivalent  to  a  third  or 
higher  order  lag  filter). 

(1)  Easic  crossover  frequency,  mc0*  The  basic  crossover  frequency 
for  quasi-rectangular  forcing  function  spectra  is  found  by 
adding  the  phase  angle,  -uxr0,  due  to  the  base  effective  time 
delay,  to  the  phase  angles  of  the  contr  lied  element  and  the 
previously  estimated  Yp  equalizer  characteristics.  Estimates 
for  o>c0  and  the  associated  pilot  gain  are  then  made  from  the 
conditions  for  neutral  stability. 

(2)  Phase  margin.  The  phase  margin  for  this  forcing  function 
category  corresponds  to  the  incremental  time  delay,  Ate,  of 
Eq.  36: 

=  (0.08-ieK.o  (37) 

(b)  low-pass  with  a  roll-off  of  less  than  third  order  and  augmented  (chslf- 
type)  continuous  input  spectra. 

(l)  Nominal  crossover  frequency,  cue*  With  equalization  and 
effective  time  delay,  Te,  selected  as  a Dove,  the  nominal 
crossover  frequency,  and  associated  pilot  gain  is  esti¬ 
mated  from  the  condition  to  provide  minimum  mean-squared 
error. 

(c)  cuq  regression-  When  mie  nears  or  becomes  greater  than  0.8  u^0  for  the 
quasi-rectangular  forcing  fui.  ttion  case  or  when  <uie/u\.  is  greater  than  1 
for  the  low-pass  and  augmented  low-pass  spectra,  then  the  crossover  fre¬ 
quency  regresses  to  values  much  lower  than  and  u\j,  respectively. 

(d)  «o  invariance  properties. 

(l)  Mo  ■*—  10.  independence.  After  initial  adjustment,  changes  in 
controlled  element  gain,  K0,  are  offset  by  changes  in  pilot 
gain,  Kp;  i.e.,  system  crossover  frequency,  tc,  is  invariant 
with  Kc. 
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(2)  at  — u>i  indepsndonce  •  System  crossover  frequency  increases 
only  slightly  with  forcing  function  bandwidth  until  cross¬ 
over  frequency  regression  occurs. 

(e)  Threshold  properties.  With  very  low  stimulus  amplitudes,  a  threshold  char¬ 
acteristic  should  be  included  in  series  with  the  pilot's  describing  function. 

Also,  when  full-attention,  nearly  continuous  control  actions  are  not  required, 
an  indifference  threshold  is  likely  to  be  present.  Both  of  these  lower  o.c 
from  what  would  be  estimated  using  the  above  adjustment  rules. 

The  coc  regression  phenomenon  mentioned  in  the  adjustment  rules  refers  to  a  reduction  of  pilot  gain  and, 
hence,  of  crossover  frequency  when  the  forcing  function  bandwidth  becomes  too  large.  The  physical  reason 
underlying  this  phenomenon  is  best  described  by  referring  to  Fig.  12.  Here,  it  is  seen  that  for  a  normal¬ 
ized  forcing  function  bandwidth,  xu^  less  than  about  0.08toc,  an  increase  in  gain  results  in  a  decrease 
in  normalized  mean-squared  error.  Wien  this  approximate  inequality  is  reversed,  the  normalized  mean- 
squared  error  becomes  greater  than  1  as  gain  is  increased.  In  fact,  ideally,  for  the  assumed  rectangular 
forcing  function  spectra  used  to  compute  Fig.  12,  a  gain  of  zero  would  be  indicated  to  minimize  the  error 
for  these  cases.  This  result  is  academic  for  pilot/vehicle  control  systems,  because  some  pilot  gain  is 
needed  to  maintain  control;  but  the  trend,  nonetheless,  for  high  forcing  function  bandwidths  is  to  reduce 
gain.  The  same  regression  effect  can  occur  for  other  than  rectangular  forcing  function  spectra  (Refs.  68, 
?8,  79).  This  regression  effect  has  practical  consequences  whenever  the  pilot  is  required  to  track  broad¬ 
band  signals. 


The  adjustment  rules  given  above  are  generally  adequate  for  the  pilot's  lower-frequency  dynamics  in 
tasks  with  spring-restrained  manipulators.  The  higher-frequency  properties  due  primarily  to  the  neuro¬ 
muscular  actuation  system  are  included  only  to  the  extent  that  is  a  component  of  Te-  The  complete 
model  of  Fig.  15  provides  for  a  much  more  elaborate  characterization  of  the  neuromuscular  system,  subject 
only  to  the  proviso  that  these  dynamics  reduce  to  the  when  viewed  from  the  low-frequency  end. 


Recent  data  (Refs.  68  and  71)  and  modeling  efforts  (e.g.,  Refs.  71  and  75-75)  can  be  used  to  estimate 
more  completely  the  neuromuscular  actuation  system  dynamics  for  a  given  set  of  manipulator  characteristics. 
An  example  of  what  can  be  involved  for  the  muscle-manipulator  dynamics  (the  forward  loop  element  in  the 
neuromuscular  system  of  Fig.  15)  is  indicated  in  the  simplified  schematic  diagram  of  the  muscle  and  mani¬ 
pulator  elements  shown  in  Fig.  21.  There  the  effective  driving  force,  CpAfa,  is  proportional  to  the  change 
in  average  firing  rate,  Afa,  of  the  alpha  motor  neuron  ensembles  involved.  This  is  essentially  the  signal  x 
in  Fig.  15.  The  muscle  characteristics  are  shown  as  functions  of  P0,  the  steady-state  isometric  tension  of 


the  muscle  system  operating  point, 
the  gamma  motor  neuron  system  dis¬ 
charge,  7Q.  The  tension  changes 
result  in  a  modified  and  thus 
underlie  the  variation  of  xe  with, 
for  example,  u4e-  Typical  describ¬ 
ing  function  data  for  the  muscle/ 
manipulator  dynamics  are  shown  in 
Fig.  22a  and  b  for  rudder  pedals 
and  a  hand  manipulator,  respectively 
(Ref.  71).  These  data  indicate  that 
the  muscle /manipulator  dynamics  for 
rudder  pedals  and  hand  manipulators 
are  similar  in  form  and  numerically, 
in  spite  of  the  difference  in  limb 
size  and  function.  The  data  also 
provide  an  exemplary  indication  of 
the  numerical  values  involved  in 
this  neuromuscular  system  component. 


The  changes  in  this  average  tension  are  caused  primarily  by  changes  in 
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Figure  21 .  Schematic  of  Limb /Manipulator 


When  the  neuromuscular  system  feedback  elements  (the  spindle/tendon  organ  ensembles  of  Fig.  15  with 
specific  numerical  values  ZSp  =  11  rad/ sec,  PST>  =  kO  rad/sec,  and  x3p  -  J‘0  msec  as  given  in  Ref.  71 )  are 
combined  with  the  muscle/manipulator  dynamics  (Fig-  22a)  into  a  closed-loop  system,  the  resulting  neuro¬ 
muscular  system  dynamics  are  shown  as  the  curves  in  Fig.  2'j.  The  data  points  on  the  figure  are  those  for 
the  total  human  operator  describing  function,  measured  for  this  controlled  element.  The  closed-loop 
neuromuscular  system  curves  provide  excellent  curve  fits  to  these  data  at  the  higher  frequencies.  On 
Fig.  25  the  closed-loop  neuromuscular  system  dynamics  are  1/Tjj)  =  12  rad/sec,  u>jj  =  19  rad/sec,  and 
tjj  2  0.15.  These  are  representative  values  for  stiff  spring-restrained,  nearly-isometric  situations 
with  very  low-inertia  manipulators. 


The  neuromuscular  system  dynamics  will  change  markedly  as  the  manipulator  load  dynamics  are  modified. 

One  of  the  most  important  of  these  possible  mod'fications  is  reduction  in  stiffness  of  the  spring  restraints. 
This  is  a  common  feature  of  aileron  controls,  as  opposed  to  elevator  and  rudder  controls.  When  the  spring 
forces  are  light  the  manipulator  approaches  the  free-moving  (isotonic)  extreme .  In  these  cases,  the 
spindle/ tendon  organ  ensemble  feedback  elements  are  joined  by  joint  receptor  ensembles.  These  feed  back 
to  highe  centers  than  the  spindle  organs  before  they  influence  the  alpha  motor  neuron  commands.  They 
accordingly  introduce  into  the  neuromuscular  system  dynamics  additional  delays  which  are  not  present  with 
the  isometric  situation.  Available  data  from  Refs.  59/  68,  and  7)  indicate  that  the  effect  of  this  pro¬ 
prioceptive  feedback  required  of  the  pilot  when  the  manipulator  is  free-moving  is  to  increase  the  effec¬ 
tive  time  delay  by  approximately  0.1  sec.  This  can  be  added  directly  to  the  previously  discussed  time 
delay,  x0,  of  Fig.  20a.  It  amounts  to  an  additional  time  delay  cost  incurred  by  forcing  the  pilot  to 
close  a  positional  loop  about  the  manipulator. 


33 


2  5  10 

w(  rad  /sec) 


50  (00 


u(rad/sec) 


(a)  Rudder  Pedal  (b)  Hand  Manipulator 

Figure  22.  Muscle /Manipulator  Describing  Function,  Gm,  for  Yc  =  l/(s-i) 


E.  REMNANT  CHARACTERISTICS 


The  remnant  component  of  the  human  pilot’s  response 
is  generally  defined  as  that  portion  not  accounted  for 
by  his  describing  function.  Vfhen  all  of  the  system 
elements  other  than  the  pilot  can  be  described  mathe¬ 
matically  in  constant  coefficient  linear  and  time¬ 
stationary.  terms,  remnant  must  arise  from  the  pilot's 
action  alone.  Then,  in  principle,  remnant  could  result 
from  the  following  sources: 

1 .  Pilot  responses  to  inputs  other  than 
the  supposed  system  forcing  function. 

2.  Nonlinear  transfer  behavior. 

3.  Non-steady  behavior. 

h .  Injection  of  "noise"  into  the  loop. 

The  first  can  be  ignored  for  the  single-loop  situa¬ 
tions  with  a  solitary  forcing  function.  For  the 
others,  the  best  data  to  consider  are  those  taken 
with  a  limited  number  of  sinuoids.  Then  the  considera¬ 
tions  of  Table  3  apply.  Existing  remnant  and  describ¬ 
ing  function  data  obtained  with  this  type  of  forcing 
function  indicate  that: 


Figure  23 .  Closed- Loop  Neuromuscular  System 
Model  and  Total  Yp  Dataj  Rudder 
Pedal,  Yc  =  l/(s—  1 ) 


1 .  Remnant  is  a  continuous  and  reasonably  smooth  spectrum  indicating  no  spectral  lines 
which  might  be  associated  with  periodic  sampling  or  strongly  nonlinear  behavior 
(Ref.  56). 

2.  For  full-attention,  continuous -control  operations  there  is  very  little  nonlinear 
dependence  of  the  describing  function  on  forcing  function  amplitude  over  a  wide 
range,  i.e.,  Yp  /  Yp(ci)  (Refs.  21  and  56).  Threshold  phenomena  are  present 
with  low-velocity  stimuli  (Ref.  7h),  and  at  very  low  display  gains  (Refs.  136). 

An  "indifference"  threshold  (Refs,  ih  and  19)  has  also  been  observed.  These  can 
be  significant  in  special  circumstances. 

3.  The  remnant  data  for  a  wide  variety  of  controlled  elements  and  forcing  function 
amplitudes  coalesce  best  when  all  the  remnant  is  reflected  to  the  pilot's  input 
(Ref.  5 6). 
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k.  When  open-loop  input-referred  remnant  is  normalized  by  the  variance  of  the  system 
error,  the  remnant  data  tend  to  coalesce  further  (Refs.  80-82).  A  sometimes- 
refined  normalization  is  with  respect  to  the  variance  of  the  principal  pilot 
channel  being  exercised,  e.g.,  o|,  when  the  rate  channel  is  used  to  control 
Yc  =  K/s2,  has  also  been  observed  (Refs.  80-8l). 

5-  The  major  sources  of  remnant  appear  to  be  located  in  the  human  operator  prior  to 
the  neuromuscular  system  (Refs.  71  and  83).  This  is  indicated  in  Fig.  22a  by  the 
interleaving  of  open-  and  closed-loop  cross-spectral  measurements  based  on  x  and 
the  forcing  function  as  references.  These  measurements  are  analogous  to  the  two 
expressions  for  Yp  presented  in  Table  k . 

6.  Some  evidence  for  pulsing  behavior  when  very-low-frequency  lead  generation  is 

required  is  present  from  output  amplitude  distributions  and  time  traces  (Ref.  56). 

When  these  findings  are  compared  with  the  entries  in  Table  3,  the  conclusion  is  obvious  that  same  variety 
of  random  time-varying  behavior  is  present  in  the  central  elements  of  the  pilot.  That  this  is  the  most 
likely  source  of  remnant  is  indirectly  corroborated  by  other  studies  (e.g.,  Refs.  8k-87).  The  time  varia¬ 
tion  could  comprise  fluctuations  in  gains  and/or  in  the  time  delays.  Figure  13  indicates  the  latter,  as 
shown  by  the  Tfj(t).  This  tentative  assignment  of  the  fluctuations  to  the  effective  delay  is  consistent 
with  trial-to-trial  variations  exhibited  in  reaction  time  measurements  but  is  otherwise  somewhat  arbitrary. 
By  considering  it  a  random  change  in  the  delay,  the  remnant  cause  can  be  interpreted  as  a  random  change  in 
phase,  akin  to  a  random  frequency  modulation,  or  to  variations  of  sampling  rate  in  a  sampled  data  inter¬ 
pretation  of  the  pilot  (Ref.  88). 

When  interest  is  centered  primarily  on  power- spectral  densities  or  their  time  domain  covariance  equiva¬ 
lents  and  the  quantities  derivable  therefrom,  the  effect  of  any  remnant  source  associated  with  the  transfer 
characteristic  can  be  modeled  as  an  injected  noise,  just  as  shown  in  Fig.  it.  Remnant  power  spectra  are, 
therefore,  most  conveniently  inserted  at  this  point  and  are  often  referred  to  accordingly  as  "observation 
noise."  While  this  use  of  the  remnant  as  an  inserted  noise  source  is  adequate  for  calculations  of  mean- 
squared  error  if  finer-grained  detail  is  needed,  such  as  the  probability  distributions  of  signal  amplitudes 
throughout  the  pilot/vehicle  system,  more  attention  has  to  be  paid  to  the  actual  remnant-generating  process 
as  a  random  time-varying  delay  or  gain. 

For  single-loop  systems  an  approximation  to  the  remnant  form,  $nne,  when  reflected  to  the  pilot's  input, 
is  given  as  (Ref.  89): 


^nrce  r  (0.1  to  0-3)  when  integral  and  proportional 

02  ~  (^2  +  ^2)  channels  are  used 

(58) 

*nn<;  .  (0.1  to  0.3)  . 

— __  =  e - - - when  the  rate  channel  is  used 

of  (a.2  +  1 ) 


As  demonstrated  in  Fig.  2k  these  analytical  forms  bound  most  of  the  available  data.  Note  that  the  remnant 
is  smaller,  although  somewhat  more  broadband,  when  low-frequency  lead  is  not  required  of  the  pilot.  Thus, 
another  penalty  for  low-frequency  lead  generation  is  seen  to  be  an  increased  remnant.  The  high-frequency 
asymptotes  for  both  no-lead  and  high-lead  situations  are  common.  Consequently,  the  remnant  increase  when 
lead  equalization  is  necessary  is  subsumed  here  primarily  by  the  reduction  in  break  frequency. 

As  already  explained,  because  the  kind  of  remnant  described  scales  with  pilot-stimulus  variance,  it 
clearly  derives  from  some  signal  conditioning  operations  within  the  pilot  such  as  the  time-varying  time 
delay.  Accordingly,  this  component  of  remnant  is  a  "processing  noise."  In  Chapter  V  another  remnant  com¬ 
ponent  due  to  scanning  will  be  introduced,  which  adds  to  the  processing  noise  in  multiloop  situations. 

If  remnant  were  totally  of  a  processing  noise  variety,  it  would  disappear  when  no  forcing  function  or 
disturbance  is  present.  There  is  a  great  deal  of  evidence,  however  (e.g.,  Ref.  61),  that  some  remnant 
remains.  This  is  wideband  and  independent  of  the  stimulus  signal  variance.  This  "residual  remnant"  is 
the  "motor"  which  keeps  the  signals  throughout  the  loop  fluctuating  in  the  absence  of  any  external  driving 
source . 

Another  occasionally  observed  remnant  component  is  dither,  that  is,  a  sinusoidal -like  oscillation 
inserted  by  the  pilot  (Refs.  1“  and  90).  This  consciously-applied  input  is  often  present  as  an  attempt 
on  the  part  of  the  pilot  to  effectively  linearize  control  system  nonlinearities. 

Although  the  primary  source  of  remnant  is  the  effective  stochastic  variation  of  pilot  characteristics 
throughout  the  measurement  run,  threshold  nonlinearities  are  occasionally  important  as  well.  One  type 
(Ref.  136)  can  be  significant  as  a  stimulus  resolution  feature  when  display  gains  are  very  low.  Another, 
more  severe,  variety  is  a  task-dependent  "indifference  threshold."  For  many  flight  control  situations  the 
pilot's  control  task  is  intermittent  adjustments  of  controls  when  the  aircraft-alone  has  departed  too  far 
from  some  preconceived,  mission-oriented  level.  The  main  effect  here  is  a  substantial  reduction  in  gain 
if  the  operation  is  viewed  as  a  continuous  process.  Oftentimes  it  may  be  more  pertinent  to  treat  these 
adjustments  as  discrete,  using  one  of  the  models  presented  in  Chapter  VI. 
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Figure  £«.  Normalized  Remnant  Spectra  for  Various  Experiments 
and  First-Order  Model  Bounds 


T.  CONNECTIONS  BETWEEN  PILOT  RATING,  WORKLOAD,  AND 
PHOT  DYNAMICS  FOR  SINGLE-LOOP  SYSTEMS 


The  previous  sections  have  provided  a  comprehensive  review  of  the  dynamic  properties  of  single-loop 
pilot  models.  These  are  sufficient,  in  themselves,  for  a  mathematical  analysis  of  pilot/aircraft  systems. 
But  pilots  are  vocal  as  well  as  dynamic  elements  —  they  toil,  and  spin,  and  talk!  And  for  many  systems 
the  talk  is  what  counts,  rather  than  the  dynamic  details  of  the  toiling  and  spinning.  Accordingly,  we 
will  present  an  Introduction  to  pilot  rating  and  workload  connections  with  pilot  dynamics.  The  discussion 
will  be  short  to  keep  matters  as  concise  as  possible)  artistic  because  ratings  are  fundamentally  ordinal 
scales  subjectively  applied  and  hence  difficult  to  quantify)  and  incomplete  because  that  portion  of  the 
data  base  which  contains  pilot  dynamics  and  pilot  ratings  and/or  comments  is  very  small  indeed. 

To  develop  closed-loop  analysis  procedures  which  permit  the  assessment  of  flying  qualities,  system 
accuracy,  etc.,  as  quantities  to  be  traded  off  with  pilot  dynamics  and  workload  requires  some  generalized 
criteria.  That  is,  the  assessment  procedures  should,  as  a  minimum,  cover: 

•  Measures  of  mission/task  performance 

•  Pilot  workload 


•  Effects  of  aircraft  dynamics 

•  Effects  of  control  augmentation  systems 

These  are  the  kinds  of  factors  which  are  taken  into  account  by  a  skilled  test  pilot  in  providing  a  pilot 
commentary  and  an  associated  rating  using,  for  example,  the  Cooper-Harper  Scale  (Ref.  93),  part  of  which 
is  shown  in  Fig.  25-  It  is  apparent  from  Fig.  25  that  pilot  compensation  (equalization)  and  effort 
(workload)  are  key  factors  in  the  rating  scale.  The  scale  is  especially  useful  as  an  index  for  comparing 
competing  vehicles  on  a  workload,  pilot  compensation,  basis. 

Closed-loop  tasks  are  ordinarily  critical  from  the  standpoint  of  pilot  compensation  or  skill  required, 
and  are  often  critically  involved  in  high  workload  phases  of  flight.  Consequently,  we  can  expect  some 
connections  between  subjective  ratings  and  the  pilot  and  pilot/vehicle  system  dynamics  and  performance. 
These  connections  are  intrinsically  empirical,  and  correlations  are  made  somewhat  awkward  because  the 
rating  scale  is  ordinal.  The  latter  point  can  be  circumvented  because  the  Cooper-Harper  scale  can  be 
related  to  an  interval  scale  (Ref.  95),  thereby  permitting  the  use  of  parametric  statistics  when  using 
and  discussing  ratings.  There  are  many  other  features  of  pilot  ratings  that  are  important  to  the  analyst) 
these  are  well  covered  from  the  pilot's  viewpoint  in  Ref.  93,  and  from  a  theoretical  standpoint  in  Ref.  95. 

When  the  tasks  being  rated  can  be  treated  as  frozen  conditions  or  are  otherwise  more  or  less  time 
stationary  (if  orly  for  a  short  interval),  a  general  form  of  rating  functional  which  explicitly  contains 
some,  and  implicitly  contains  all,  of  the  desired  features  is  given  by: 
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Figure  25.  Cooper-Harper  Handling 
Qualities  Rating  Scale 


The  subscript  notation  used  here  implies  that  i  €  I 
motion  and  task  measures  are  controlled  by  k  £  K 
pilot  loops  actuating  j  €  J  control  points.  As  dis¬ 
played  in  Eq.  39  the  functional  form  is  general 
enough  to  include  the  existing  (e.g..  Refs.  9)1-100) 
approaches  to  quantitative  flying  qualities  rating 
criteria  functions.  The  key  closed- loop  system 
quantities  in  the  rating  functional  are  measures  of 
mission  and  task  performance.  These  are  conveniently 
described  by  a  set  of  dominant  weighted  (via  Ai)  air¬ 
craft  motion  deviations  and  total  task  accuracy  or 
error  indications. 

__  The  pilot  activity  component  of  pilot  effort, 

&?,  is  particularly  dependent  on  the  level  of  pilot 
gain.  For  a  given  gain,  tS  increases  directly  with 
gust  disturbance  spectrum  amplitude  and  remnant 
amplitude.  Accordingly,  both  and  the  qf  quanti¬ 
ties  will  reflect  turbulence  and  remnant  levels. 

The  pilot  equalization  component  of  pilot  work¬ 
load  is  represented  in  Eq.  39  by  the  slope  (in  dB 
per  octave  or  decade)  of  the  pilot's  amplitude  ratio 
evaluated  at  a  particular  frequency  (generally  near 
crossover).  This  is  by  no  means  the  only  measure 
available  to  describe  the  dynamic  quality  of  the 
pilot’s  effort.  References  96-100,  for  example,  use 
pilot  lead  time  constants  as  measures;  for  particular 
situations  with  a  sufficient  data  base,  this  may  be 
a  desirable  alternative.  Then,  the  rating  functional 
can  take  a  form  such  as  illustrated  in  Fig.  26. 
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Figure  26.  Typical  Pilot /Vehicle  Model  and  Form  of  Pilot  Rating  Functional 
(Specific  Values  for  Hover  Task,  Refs.  96-97) 

(Refs.  5-6) 
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Unfortunately,  there  do  not  now  exist  for  all  closed-loop  tasks  set.-  ■  f  pilot  rating,  system  perfor¬ 
mance,  and  pilot  equalization  data.  Consequently,  for  many  situations  we  must  rely  on  available  correla¬ 
tions  which  do  no*-  include  quantities  like  the  q?. 

The  pilot  adapts  to  the  vehicle  and  forcing  function  characteristics.  He  therefore  reflects  in  his 
adapted  describing  function  form  many,  if  not  all,  of  the  vehicle  dynamic  characteristics  and  closed-loop 
pilot; vehicle  system  properties.  Consequently,  as  a  first  approximation,  a  functional  relationship  can  be 
set  up  between  piiot  ratings  and  the  objective  system  factors  in  terms  of  the  pilot  dynamic  characteristics 
alone.  In  this  connection,  it  has  been  found  that  the  dominant  rating-sensitive  pilot  parameters  are  the 
low-frequency  lead  equalization  and  the  crossover  gain.  Typical  relations  of  this  nature  are  shown  in 
Fig.  c,7 •  These  can  be  used  directly  to  estimate  pilot  rating  onee  the  pilot  dynamics  estimates  are  made. 


Figure  27.  Pilot  Rating  Decrements  as  Functions  of  lead  Equalization  and  Gain 


The  descriptive  phrases  listed  under  "Demands  on  the  Pilot  in  Selected  Task  or  Required  Operation"  in 
the  Cooper-Harper  Scale  (Fig.  )  directly  parallel  the  correlations  with  pilot  lead  equalization.  There 
is  also  a  strong  connotation  of  increasing  pij.ot  effort  and  w-  r’.load  in  the  scale  phrases.  But,  workload 
is  difficult  to  define  and,  consequently,  to  quantify.  In  the  spirit  of  offering  a  general  definition 
which  can  be  measured  and  predicted,  it  lias  been  suggested  (Ref.  10 1 )  that  workload  margin  he  defined  as 
the  ability  (or  capacity)  to  accomplish  additional  (expected  or  unexpected)  tasks.  For  example,  the  pilot 
opinion  rating  scale  satisfies  this  definition  up  to  its  "uncontrollable"  limit  point.  Furthermore,  a 
number  of  auxiliary'-  tasks  also  satisfy-  this  definition  of  workload  in  that  the  decrements  in  auxiliary 
tasr.  scores  give  an  index  of  demand  on  the  primary  task.  One  particular  measure  offers,  at  the  moment, 
unusual  promise  in  integrating  many  of  the  measures  into  one  basic  context.  This  is  excess  control  capa¬ 
city,  a  major  connector  with  pilot  rating  and  main  task  effective  time  delay(s). 

The  notion  what  among  the  causal  factors  of  pilot  rating  is  the  pilot's  attention  or  effort  need“d  to 
maintain  performance  is  supported  by  an  experiment  which  measured  a  parameter  uniquely  related  to  excess 
control  capacity  (Ref.  9‘  ).  A  secondary  subcritical  tracking  task  was  used  to  "load"  the  pilot  so  that 
his  performance  on  the  primary  task  began  to  deteriorate.  A  block  diagram  of  these  tasks  is  shown  in 
Fig.  28.  The  difficulty  of  the  secondary  task  was  made 
proportional  to  primary  task  performance.  Thus,  when 
the  pilot  was  keeping  primary  task  error  performance 
less  than  a  criterion  value,  the  secondary  task  diffi¬ 
culty  was  automatically  increased  by  increasing  the 
rate  of  divergence  of  the  secondary  instability.  Con¬ 
versely,  when  the  pilot  was  so  busy  with  the  secondary 
task  that  ordinary  error  was  larger  than  the  criterion 
value,  the  secondary  task  difficulty  automatically 
decreased.  The  final  stationary  level  of  secondary 
difficulty  was  determined  by  the  sensitivity  of  the 
primary  task  performance  to  loading.  The  final"score" 
is  -s,  the  stationary  value  of  the  secondary  unstable 
pole  (■)  in  rad/sec.  The  scores  obtained  from  this 
cross-coupled  secondary  task  represent  its  degree  of 
difficulty;  consequently,  they  also  represent  the 
"degree  of  ease"  of  the  primary  task  or  the  excess 
control  capacity  available  with  respect  to  the  pri¬ 
mary  task. 

The  achievement  of  the  critical  limiting  score  for 
which  >c  -  ]Ae  the  cross-coupled  secondary  task 
indicates  a  condition  of  maximum  available  excess  con¬ 
trol  capacity.  We  speak  of  the  secondary  task  as  a 
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Figure  28.  Single-Loop  Primary  Task  with 
Secondary  Cross-Coupled  Loading  Task 


"critical"  task  in  this  Uniting  case  where  Te  is  the  sole  determinant  of  system  stability.  Thus,  any 
activity  by  the  pilot  which  demands  an  increase  on  re  on  the  whole  task  can  be  expected  to  prevent  him 
from  achieving  his  critical  limiting  score  on  the  cross-coupled  secondary  task. 


Secondary  scores  obtained  for  a  variety  of  pri¬ 
mary  controlled  elements  are  presented  in  Ref.  91>. 
Figure  29  shows  how  the  scores  for  the  best  gain 
configurations  of  each  controlled  element  compare 
with  the  Cooper  ratings.  The  agreement  is  extremely 
good.  Even  the  subcritical  task  itself  in  the  role 
of  the  primary  task,  which  has  been  a  notable  cul¬ 
prit  in  other  correlations,  seems  to  be  correlated 
linearly  with  the  other  data.  In  Fig.  29  a  score 
>s  =  0  corresponds  to  100  percent  of  the  pilot's 
attention  being  devoted  to  the  primary  task  or  no 
excess  control  capacity,  whereas  a  limiting  score 
(,\s  =  !r.f')  means  that  no  attention  is  required  to 
maintain  primary  task  performance  or  that  lOO  per¬ 
cent  excess  control  capacity  is  available. 

These  limited  experimental  data  offer  a  com¬ 
pelling  direct  connection  of  pilot  rating  to 
measures  of  workload,  thereby  quantifying  ordinal 
scale  pilot  ratings  in  these  workload  terms  for 
continuous  control  tasks.  The  particular  experi¬ 
ment  on  which  these  nice  correlations  are  based 
is,  unfortunately,  difficult  to  perform  and  general¬ 
ize.  A  simpler  form  wherein  the  side  task  is  not 
cross-coupled  but  instead  is  used  to  set  a  series 
of  subsidiary  workload  levels,  each  calibrated  by 
a  value  of  >3,  is  much  easier  to  apply. 
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Figure  29-  Subjective  Pilot  Rating  Versus 
First-Order  Cross-Coupled 
Instability  Score 
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CHAPTER  V 

QUASI-LIKEAR  MODELS  FCR  COMPENSATORY  MULTIIOOP  SYSTEMS 


The  multiloop  system  pilot  models  discussed  in  this  chapter  are  fundamentally  extensions  of  the  quasi- 
lincar  models  for  compensate ry  single-loop  systems  treated  in  Chapter  XV.  To  ease  the  transition  from  the 
single-loop  case  we  vill  first  describe  some  of  the  needs  for  multiloop  control  in  aircraft.  This  will  be 
followed  by  a  review  of  the  data  base  for  visual  inputs.  A  continuous-attention  multiloop  model  for  non¬ 
scanning  visual  situations  is  then  presented,  with  emphasis  on  the  adjustment  of  the  multiloop  structure. 

This  basic  multiloop  model  has  several  variants,  depending  on  the  nature  of  its  perceptual  stages  ir.  a 
particular  instance  or,  in  plain  words,  depending  on  how  much  scanning  is  involved.  The  first  modifica¬ 
tion  to  the  multiloop  model  ‘akes  scanning  into  account  in  a  first-order  fashion  via  modifications  to  loop 
gains  and  remnant. 

The  second  modification  to  the  basic  multiloop  model  differs  from  the  first  in  that  the  feedbacks 
derive  from  proprioceptive-vestibular  signals  associated  with  vehicle  motion.  This  is  essentially  a  parti¬ 
cular  "multimodality"  model  where  pilot  motion  inputs  Join  vision  in  providing  control  cues  on  which  the 
pilot  can  operate. 

Unlike  the  single-loop  models  which  are  based  on  three  decades  of  evolutionary  improvement,  some  features 
and  adjustment  rules  for  the  multiloop  models  are  quite  new  and  not  yet  firmly  based  on  extensive  data.  Con¬ 
sequently,  we  can  expect  much  that  is  written  here  to  change  in  the  next  few  years.  Perhaps  nowhere  is  this 
more  likely  to  occur  than  with  the  very  tentative  connections  of  multiloop  pilot  dynamics  with  pilot  rating 
described  in  the  section  of  the  chapter.  Nonetheless,  the  quasi-linear  pilot  models  for  multiloop  systems 
described  here  "fly"  tracking  stages  (via  multiloop  describing  functions  which  close  loops  about  the  air¬ 
plane),  introduce  pilot-induced  noise  (remnant)  in  the  process,  look  about  to  gather  information  and, 
incidentally,  add  more  noise  (scanning  traffic  and  scanning  remnant),  have  some  workload  margin  (excess 
control  capacity),  and  "speak"  rationally  (pilot  rating). 

A.  NEED  IDS  MULTILOOP  CONTROL 


At  the  outset  it  is  desirable  to  define  what  we  mean  by  multiloop.  This  can,  perhaps,  best  be  accom¬ 
plished  with  the  aid  of  rig.  '0  (Ref.  r7)-  The  first  block  diagram.  Fig.  JOa,  depicts  the  single-loop 
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Figure  JO.  Examples  of  Single  and  Multiloop  Manual  Control  Systems 
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configuration.  The  second  system  shown  consists  of  two  single -loup  compensatory  systems  with  no  coupling 
in  the  controlled  element.  A  typical  example  of  such  a  multiaxis  situation  is  longitudinal  and  lateral 
stabilization  of  an  aircraft  in  straight,  wings-level,  horizontal  flight  using  aileron  and  elevator.  A 
number  of  the  studies  contained  in  Table  C  were  for  controlled  elements  of  this  nature  (Refs.  >9,  lC-lS, 

2',  C2-0*,  and  69).  This  multiple  single-loop  structure  was  assumed  (for  data  analysis  purposes)  and 
demonstrated  to  be  appropriate  (as  based  on  data  results)  for  these  cases.  Multiloop  systems  differ  from 
multiple  single-loop  systems  in  that  interaxis  coupling  is  inherent  in  tne  controller  or  the  controlled 
element  dynamics,  as  indicated  in  Fig.  JOc-f.  These  Kinds  of  systems  represent  a  higher  level  of  pilot/ 
vehicle  complexity  but  are,  nonetheless,  required  for  many  piloting  tasks. 

Fundamentally,  there  are  three  inspirations  for  multiloop  control.  The  first  is  a  desire  to  exert 
control  over  more  than  one  variable}  in  general,  this  requires  one  point  of  control  application  per  con¬ 
trol  variable  as  well  as  meeting  controllability  and  observability  limitations.  Although  this  desire  could 
hypothetically  be  met  with  multiple  single-loop  systems,  a  true  multiloop  <:ontr..±  is  often  the  practical 
result. 

The  second  reason  is  the  use  of  auxiliary  quantities  (controlled  element  outputs)  as  feedbacks  in  lieu 
of  series  compensation  of  primary  quantities.  Common  examples  of  thiw  parallel  ( feedbag. i  equalization  are 
the  use  of  pitch  attitude  (instead  of  altitude  rate)  to  supply  path  damping  in  an  altitude  control  system, 
and  bank  angle  (instead  of  heading  rate)  to  provide  path  damping  for  a  heading  control  system.  Often,  the 
auxiliary  feedback  will  have  advantages  in  other  respects.  For  instance,  it  may  provide  a  more  stable 
feedback  for  other  than  the  primary  modes  considered,  may  be  easier  to  sense,  may  suffer  less  from  noise 
contamination,  or  may  be  itself  a  suitable  outer-loop  feedback  for  some  flight  control  system  mode.  For 
piloted  control  systems,  the  auxiliary  feedback  in  lieu  of  series  compensation  can  be  profoundly  important. 
This  is  perhaps  best  appreciated  by  recalling  that  pilot  generation  of  lead  equalization  incur.,  penalties 
in  incremental  time  delay,  increased  remnant,  decreased  system  performance,  increased  pilot  workload,  and 
poorer  pilot  ratings. 

The  third  reason  for  multi loop  control  Is  to  a-hieve  coupling  or  decoupling  purposes.  This  is  a  modi¬ 
fication  of  the  effective  controlled  element  transfer  function  numerators  by  auxiliary  .. ,ntrol  from  another 
control  point. 

Because  this  list  of  inspirations  for  nultiloop  control  is  particularly  favorable  for  manual  control, 
piloted  situations  may  more  often  be  multiloop  than  corresponding  automatic  conditions. 

The  key  to  multiloop  pilot  models  and  actions  is  also  the  first  fundamental  concept  of  pilot 'vehicle 
analysis:  that  the  pilot  "constructs"  feedback  loops  about  the  effective  -.-ntroiled  element.  The  feedback 
quantities  actually  selected  by  the  pilot  wo.ll  be  "those  necessary  to  satisfy  the  guidance  and  control  needs 

and  certain  pilot-centered  requirements.  The  guidance 
and  control  needs  are  situation-specific.  Satisfac¬ 
tion  of  these  needs  always  involves  a  task-  or  purpose- 
centered  outer  loop,  with  possible  subsidiary  inner- 
loop  and  other  axis  closures  as  needed  to  make  the 
principal  feedbacks  work.  Examples  of  appropriate 
outer  and  subsidiary  loo:  structures  for  typical  pre¬ 
cision  tracking  and  precise  flight  path  control  tasks 
are  given  in  Fig.  i 1 .  Illustrated  there  are  system 
configurations  pertinent  to  air-to-air  combat  and/or 
weapon  delivery  and  to  formation  flight.  The  task  and 
purpose-centered  outer  loops  are  shown  <■ ith  solid  lines 
in  this  figure,  whereas  the  subsidiary  inner  loops  are 
shown  with  dashed  lines.  For  the  tasks  noted,  the 
figures  may  not  contain  all  the  inner  loops  and  they 
do  not  indicate  the  crossfeeds  or  other-axis  closures 
which  might  be  desirable  from  the  standpoint  of  gui¬ 
dance  and  control  needs.  These  additional  features  are, 
more  often  than  not,  highly  vehicle-specific.  None¬ 
theless,  the  unlabeled  blocks,  each  of  which  represents 
pilot  actions  and  thus  a  describing  function,  are  suffi¬ 
ciently  large  in  number  to  imply  an  enormously  complex 
analysis  and  measurement  situation.  Fortunately,  the 
analytical  complexity  is  more  apparent  than  real,  since 
there  are  efficient  analytical  techniques  available 
specifically  designed  to  handle  these  types  of  systems 
(Ref.  102). 


In  multiloop  pilot/aircraft  control  systems,  the 
essential  features  of  the  system  structures  are  the 
feedbacks  themselves  and  their  equalization.  For  suc¬ 
cessful  systems,  i.e.,  systems  which  demonstrate  uni¬ 
form,  reliable,  high-quality  performance  in  a  given 
task,  etc.,  the  possible  feedback  structures  are  very 
limited.  They  derive  primarily  from  guidance,  control, 
and  regulation  demands,  and  secondarily  from  dynamic 
response  characteristics  desired  by  the  pilot.  From 
the  systems  view  it  is  the  satisfaction  of  these  require¬ 
ments  that  is  important  rather  than  the  means  employed. 

In  other  words,  the  feedback  loops  closed  is  the  cen¬ 
tral  issue,  whether  the  closures  are  accomplished  auto¬ 
matically  or  manually. 


CATAGORY  A  TASKS  ( Ropid  Moneuvering ,  Precision  Tracking , 
Precise  Plight  Path  Control  ) 


a)  Offensive  Air-to-Air  Combat  or  an  Air -to- Ground 
Weapon  Delivery 


b)  Defensive  Air ->o- Air  Combat 


c)  Formation  Flight 


Figure  31 •  Illustration  of  Outer  and  Subsidiary 
Loops  for  Typical  Precision  Tracking  and 
Precise  Flight  Path  Control  Tasks 
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Stated  verbally,  the  key  guidance  and  control  requirements  for  most  pilot/aircraft  precision  control 
tasks  may  include: 

•  Establish  and  maintain  the  aircraft  on  some  specified  spatial  pathway  or  beam  (e.g., 
localizer  and  glide  path). 

•  Reduce  flight  path  errors  to  zero  in  a  stable,  well-damped,  and  rapidly  responding 
manner. 

•  Establish  a  (perhaps  accelerating)  equilibrium  flight  condition. 

•  Limit  the  speed  or  angle  of  attack  excursions  from  this  established  equilibrium 
flight  condition. 

The  regulation  requirements  are  similar,  i.e.: 

•  Maintain  the  established  flight  path  in  the  presence  of  disturbances  such  as  gusts, 
crosswinds,  and  wind  shears. 

•  Provide  a  degree  of  short-time  attitude  stability  in  the  presence  of  disturbances. 

These  re  juirements  relate  primarily  to  the  relatively  low-frequency  path  modes  of  the  pilot/aircraft  system. 
In  essence,  they  define  outer  control  loops  involving  those  vehicle  motion  quantities  which  define  the 
desired  equilibrium  state  of  motion.  More  often  than  not,  such  outer  loops,  when  closed  about  unmodified 
aircraft  dynamics,  do  not  result  in  stable,  weld-damped,  rapidly  responding  systems.  Instead,  as  already 
remarked,  equalization  of  either  a  series  or  a  parallel  nature  is  needed  to  assist.  Parallel  equalization 
is  cost  common  and  is  achieved  by  the  use  of  inner  loops  which  feed  back  such  quantities  as  attitude, 
angular  velocity,  and  sometimes  linear  acceleration.  These  inner  loops  dominate  the  high-frequency  char¬ 
acteristics  of  the  aircraft  ^controller  system. 


3.  MULTILOOP  DATA  BASE 

The  data  base  available  for  the  construction  of  mathematical  pilot  models  for  multiiovp  compensatory 
situations  ..ith  visual  inputs  is  given  in  Tabic  j.  Because  of  the  instrumental  measurement  and  analytical 
difficulties  (Ref.  •'))  inherent  i:.  finding  unique  pilot  describing  functions  in  multiloop  systems,  only 
these  three  studies  have  been  accomplished.  However,  their  combined  scope  is  sufficient  to  cover  many 
cases  of  pilot 'aircraft  system,  interest,  both  laterally  and  longitulinally .  Several  other  experimental 
series  avoid  some  of  the  measurement  and  analysis  difficulties  by  assuming  adjustable  fixed-form  pilot 
characteristics  (e.g.,  Refs.  100,  I0  -lOy).  These  provide  useful,  and  in  some  cases  definitive,  data  even 
though  the  pilot  dynamic  characteristics  are  not  uniquely  defined. 

Another  aspect  of  tr.e  Table  0  data  which  significantly  enhances  their  utility  and  broadens  their  scope 
is  that  each  seriec  was  preceded  by  an  extensive  pre-experimental  analysis  activity.  Those  preliminary 
studies  estimated  the  pilot  dynamics  and  pilot/ vehii.j.0  system  dynamics  arid  performance  before  the  experi¬ 
ments  were  conducted  uing  the  then  best  available  version  of  the  multiloop  pilot  model.  The  estimates 
..ere  used  both  to  help  configure  critical  aspects  of  the  experiments  and  as  the  basis  for  the  setup  of 
analog  pilots  to  permit  the  assessment  of  .measurement  and  data  analysis  techniques,  after  the  experiments 
had  been  completed  and  tne  data  analysed,  it  was  possible  to  compare  the  pre-expenmertul  estimates  with 
the  actual  experimental  results.  Difference-,  bet. .ecru  expectation  and  realization  then  Lead  tv  rectifica¬ 
tion  in  tertv;  of  the  model,  as  well  as  to  a  more  direct  understanding  of  the  results,  ’his  application  of 
I ilot, 'vehicle  system  analysis  to  provide  a  pre-experimental  estimate  of  what  to  expect  and  for  post- 
experimental  rationalization  of  what  occured  is,  of  course,  one  of  the  principal  uses  of  pilot- vehicle 
analysis;  but  nowhere  is  this  more  advantageous  than  when  dealing  with  new  model  development. 

C.  MULTILOO?  VISUAL  COMPENSATORY  SYSTEM  PILOT  MODEL 

As  a  eou.W'i'iense  of  the  kinds  of  studies  referred  to  in  the  last  section,  a  series  of  adjustment,  rules 
similar  to  those  used  for  the  single-1.. op  model  can  be  stated.  These  are  listed  b.  ow . 

1.  The  multiloop  situation  with  full  visual  field  is  similar  to  that  shown  in  Fig.  1‘, 
if  the  signals  i,  e,  c,  and  m  are  interpreted  as  vector  quantities,  and  with  the 
understanding  that  the  perceptual  portions  of  the  central  elements  for  each  stimulus 
used  operate  in  parallel. 

2.  The  feedback  quantities  available  to  the  pilot  for  possible  use  consist  of  those: 

a.  Directly  sensed  within  the  general  visual  field. 

b.  Observable  via  visual  displays. 

c.  Directly  sensed  using  modalities  other  than  vision. 

Quantities  which  can  be  perceived  from  the  full  visual  field  or  other  modalities 
will  show  no  scanning  penalties  whereas  those  which  require  instrument  scan  or  modi¬ 
fication  of  the  eye  fixation  point  will  introduce  decrements  (described  later). 

; .  The  feedback  loops  preferred  are  those  which: 

a.  Can  be  closed  with  minimum  pilot  equalization. 

b.  Require  minimum  scanning  to  sense  the  feedback  quantity. 

c.  Permit  wide  latitude  in  the  pilot's  adapted  characteristics. 
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TABLE  9.  SUMMARY  OF  MULTILOOP,  VISUAL,  COMPENSATORY  SYSTEM  DESCRIBING;  FUNCTION  AND  REMNANT  DATA 


INVESTIGATOR 

CONTROLLED  ELEMENT  DYNAMICS 

AND  CONTROL  TASK(S) 

GENERAL  RESULTS  AND  REMARKS 

Stapleford, 
McRuer,  and 
Magdaleno 
(Ref.  '  7) 

Single -input,  two-control  point, 
integrated  display  aircraft 
lateral  control. 

Bank  angle  tracking  with  good 
laterul  (roll  subsidence  l/Tj>  A  9) 
dynamics  and  5  different  levels  of 
dutch  roll  damping  ratios: 

Cp  =  0.7?  (single-loop  control), 

PR  r  A;  =  -0-07 3  (marginal 

multiloop  required),  PR  A  6;  and 
tD  =  -0*35  (near  nultiloop  control 
limits),  PR  A  8. 

Effective  system  was  Fig.  JOf  with 
ij  =  Yp,rt  =  0. 

1 .  Single-loop  pilot  model  is  applicable  to 
multiloop  system  command  (outer)  loop  and, 
with  possible  re  ervation,  to  inner  loops 
as  well. 

2.  When  several  feedback  possibilities  are  pre¬ 
sent,  the  pilot  will  select  those  which  per¬ 
mit  best  dynamic  performance  with  least 
pilot  equalization  and  effort. 

3-  Control  crossfeeds  were  adopted  to  reduce 
inadvertent  oscillation  of  subsidiary  mode. 

i. .  Pilot  outer-loop  dynamics  for  all  config¬ 
urations  were  similar,  but  ratings  degraded 
as  required  level  of  multiloop  activity 
increased. 

Stapleford, 
Craig,  and 
Tennant 
(Ref.  C9) 

Two-input,  single-control  point, 
integrated  display,  aircraft 
longitudinal  control.  Altitude 
(approach)  control  in  the  presence 
of  attitude  and  gust  disturbances. 

Effective  system  was  Fig.  ?0d. 

) .  Attitude-nlor.e  (single-loop  task)  and  atti¬ 
tude  inner-loop  (multiloop  task)  closures 
were  very  similar. 

2.  Successive  closure  of  single-loop  models 
are  appropriate  for  multiloop  situations. 

3.  Series  (Fig.  ?0d)  closure  model  is  most 
appropriate  in  that  describing  function 
data  then  take  their  simplest  form. 

Weir  and 
McRuer 
(Ref.  10?) 

Single-loop,  single-control  point, 
integrated  display,  pitch  attitude 
control}  effective  system  Fig.  *0a. 

2  single-loop,  two-control  point, 
integrated  (flight  director)  plus 
full  panel  display  all-axis 
approach  task}  effective  system 

Fig.  ?0b  with  scanning  added. 

Multi -input,  single-control  point, 
full  panel  display,  longitudinal 
approach  control  in  the  presence 
of  attitude  and  gust  disturbances 
(lateral  axes  under  autopilot  con¬ 
trol);  effective  system  Fig.  '9d 
with  scanning  added. 

Multi-input,  two-eontrol  point,  full 
panel  display,  all-axis  approach 
task  with  attitude  and  gust  dis- 
turoanees;  effective  systems  for 
lateral  and  longitudinal  each 
correspond  to  Fig.  70d  with  scan¬ 
ning  added. 

Eye  movements  measured  as  well  as 
dynamic  properties. 

1.  Pitch  attl  ude-alone  data  are  consistent 
with  single-loop  model,  but  lower  in  gain 
than  the  Ref.  ($  results. 

2.  Data  for  pitch  attitude-alone,  flight  direc¬ 
tor  all-axis  control,  and  one  longitudinal 
control  eases  are  reasonably  well  fitted  by 
crossover  model  forms. 

?.  The  two  pilot  subjects  exhibited  different 
styles:  relatively  low-gain  inner  and  high- 
gain  outer  loops  as  contracted  to  high-gain 
inner  and  low-gain  outer  loops. 

«•  Results  were  consistent  with  satisfying 
the  guidance  and  control  requirements. 

*  •  Series  (Fig.  'Od)  structures  appear  per¬ 
tinent  for  both  lateral  and  longitudinal 
control  operations. 

( .  Both  inner  and  outer  loop  gains  for  this 
study  were  lower  than  the  comparable 

Ref.  69  results. 

7.  The  major  contributor  to  remnant  appears 
to  be  inner-loop  operation. 

<•.  Where  distinct  inner-  and  outer- loop  closures  can  be  defined  by  ordering  the  band- 
wjdths  (e.g.,  the  higher  the  bandwidth,  the  more  inner  the  loop),  a  series  aultiloop 
structure  applies. 

.  Pilot  equalization  for  the  outer  loop  of  aultiloop  systems  is  adjusted  per  the  cross¬ 
over  model,  with  the  proviso  that  the  effective  controlled  element  transfer  function 
include  the  effects  of  all  the  inner-loop  closures.  Tho  crossover  model  is  also 
directly  applicable  to  many  inner-loop  closures. 

Crossfeeds  are  commonly  adopted  by  the  pilot  to  directly  negate  the  excitation  of 
subsidiary  coupled  modes. 

7.  The  adjustment  of  the  variable  gains  in  each  of  the  loops  is,  in  general,  such  as  to 
achieve  basically  simple  (i.e.,  effectively  second-  or  third-order)  well-damped  domi¬ 
nant  modes,  and  nearly  uncoupled  sets  cf  aircraft  responses.  Outer-loop  gains,  in 
particular,  may  be  lower  than  maximum  1  andwidth  for  this  reason . 

3.  When  scanning  is  not  present,  the  remnant  is  primarily  associated  with  the  inner  loop 
and  is  essentially  '.he  same  as  that  for  a  single-loop  system  equivalent  to  the  inner 
loop  alone.  (Scanning  remnant  is  discussed  below.) 
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Although  scanning  is  avoided  where  possible,  multiple  fixations  are  sometimes  required  to  sense  the 
appropriate  controlled  element  output  quantities.  The  basic  process  of  scanning  during  multiloop  control 
tasks,  sampling  the  fixated  and  parafoveal  information,  and  reconstructing  the  scanned  signal  is  very 
complex  and  little  understood  in  detail.  However,  the  essence  of  past  work  in  this  area  (Refs.  91 >  10J, 
110-11*)  shows  that  in  the  process  of  extracting  the  feedback  information  from  the  displays: 

1 .  A  fairly  stationary  scanning  strategy  evolves  for  a  given  task  and  full  field/ 
instrument  array. 

2.  The  pilot’s  output  control  notions  are  much  more  continuous  than  a  discrete  sampling 
of  input  signals  coincident  with  eye  fixations  would  seem  to  imply  from  the  pure 
stimulus-response  sequence. 

*.  The  first-order  effects  of  scanning  are  to  reduce  the  pilot  gain  and  increase  remnant 
in  the  scanned  channels.  “  ”™  ~  “ 

The  development  of  multiloop  pilot  models  complete  with  scanning  effects  has  gone  through  several 
stages  of  theoreti:al  analysis  and  experimental  validation.  Many  of  the  phenomena  observed  empirically 
can  be  modeled  theoretically  with  two  different  multi-input,  multi-output  pilot  model  forms  called, 
respectively,  the  "switched  gain"  model  and  the  "reconstruction  hold"  model.  The  current  version  of  these 
models  is  described  in  Ref.  106.  We  shall  consider  here  the  gist  of  the  switched  gain  model  only,  because 
it  is  appropriate  to  most  critical  multiloop  piloting  situations  where  scanning  is  relatively  limited. 

This  model  is  also  conceptually  the  simpler  of  the  two. 

This  form  cf  scanning  model  is  termed  /.witched  gain  because  it  incorporates  a  quasi-randcai,  finite- 
ivell  sampling  or  switching  process  between  the  pilot's  foveal  gain  and  his  effective  parafoveal  gain  on 
each  of  the  several  displays  Involved.  Figure  J2  illustrates  the  model  with  a  block  diagram.  The  foveal 
path  is  closed  during  the  foveal  dwell  interval,  and  the  parafoveal  path  is  closed  during  the  foveal  inter¬ 
rupt  interval.  Each  of  these  paths  will,  in  general,  exhibit  different  gains,  equalization  and  effective 
time  ielaj  3  before  the  paths  are  combined  in  the  higher  neural  centers  to  send  a  signal  to  the  actuation 
describing  function . 


Figure  5?.  "Switched  Gain"  Multiaxis  Scanning  Model  for  Compensatory 
Multiaxis  Tracking  with  One  Among  Several  Displays 


The  conceptual  block  diagram  in  Fig.  12  can  be  remarkably  simplified  by  recalling  (Ref.  ilk)  that  any 
quasi-randomly  sampled  and  processed  signal  can  be  modeled  by:  1)  replacing  the  sampling  or  switching 
process  b j  a  continuous  transmission  path,  and  2)  adding  an  uncorrelated  wideband  noise  process  which  has 
a  power-spectral  density  proportional  to  the  variance  of  the  (displayed)  signal  before  sampling  Since 
the  quasi-random  scanning  process  has  a  finite  foveal  dwell  interval,  the  wideband  noise  process  will 
exhibit  a  low-pass  power  spectrum  with  a  first-order  break  frequency  which  is  inversely  proportional  to 
the  average  fove  ’  dwell  interval  (Ref.  112).  The  power- spectral  density  of  this  foveal/parafoveal 
switched  gain  scanning  remnant  is  given  in  simplified  terms  by: 
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where: 

0“  is  mean-squared  value  of  the  signal  scanned 
?s  is  mean  scanning  interval 
Tfe  is  effective  dwell  fraction  =  Td^/'fs 
Tde  is  effective  dwell  interval 

0Tg  is  standard  deviation  in  Ts 

(1-6)  is  approximately  Cji„/Ts,  the  scanning  variability  rg.tlo 


Measurements  of  this  switched  gain  remnant  in  Ref.  91  have  shown  that  it  is  so  predominant  compared  with 
the  other  sources  of  remnant  that  the  other  sources  cannot  even  be  identified.  This  makes  for  great  simpli¬ 
fication  of  the  remnant  In  the  equivalent  switched  gain  model. 

Representation  of  the  pilot's  describing  function  in  the  ..witched  gain  model  can  also  be  greatly  . impli- 
fied.  The  foveal  gain  exceeds  the  parafoveal  gain  in  all  measurements  which  have  been  made  (Refs.  91,  110- 
I1l).  This  is  probably  because  of  the  large  displacement  and  increased  rate  thresholds  in  parafoveal  per¬ 
ception  by  comparison  with  foveal  perception .  The  switched  gain  model  is  represented  simply  by  multiplying 
the  ratio  of  parafoveal  gain  to  foveal  gain  (f.)  by  the  interrupt  fraction  (l  -  t,)  and  adding  the  product  to 
the  dwell  fraction  (r,)  to  obtain  the  effective  dwell  fraction,  viz.: 


he  =  h+fi(l-*i)  (“0 

where  f  -  -Cp/‘-cf  =  rat.”  j  of  crossover  gains  for  continuous  parafovea^  tracking  relative  to  continuous  foveal 
tracking  (0  •  f  <  1 ) .  Ihe  effective  crossover  gain  for  the  equivalent  switched  gain  model  for  one  foveal/ 
parafoveal  channel  is  ’,eucf,  where  is  the  foveal  crossover  gain  in  continuous  single-axis  tracking  of 
the  same  display  an'  controlled  element  constrained  by  the  same  task  variables.  (For  a  step-by-step  appli¬ 
cation  of  the  sca-'-.ing  model  see  Ref.  137  •) 

There  are  no  apparent  phase  penalties  associated  with  switched  gain  scanning  as  long  as  _arafoveal  per¬ 
ception  is  not  completely  inhibited.  Inhibition  can  occur  either  by  requiring  a  multitude  of  different 
widely-separated  fixations  with  a  time  constraint  or  by  inducing  "tunnel  vision"  on  one  or  two  displays. 

Even  so,  measurements  reported  in  Ref.  9';  where  parafoveal  perception  was  inhibited  by  blanking  the 
parafoveally-viewed  display,  show  only  small  effective  time  delay  increments  (<1ts)  on  the  order  of  0.0>  to 
0.1J  sec  attributable  to  scanning  as  the  parafoveal-to-foveal  gain  ratio  (f.)  approached  zero. 

The  switched  gain  model  has  been  quite  successful  in  modeling  behavior  on  a  main  task  in  laboratory 
experiments  with  induce!  natural  scanning  between  a  primary  tracking  tasx  and  a  secondary  subcrltical  track¬ 
ing  task  (Ref.  91 )  and  on  foveal  and  parafoveally-viewed  displays  (Refs.  1 1 0- ill).  For  most  critical  multi¬ 
loop  aircraft  control  tasks  the  scanning  involved  ordinarily  needs  to  consider  only  two  or  three  major 
elements  (e.g.,  other  aircraft  in  tail  chase,  own  wing  man)  within  the  scan  pattern,  so  this  model  should 
be  very  useful. 

The  principal  effects  of  scanning  are  thus  seen  to  be  loop  gain  reductions  and  remnant  increases.  Both 
factors  separately  degrade  performance,  and  their  combined  effect  can  be  devastating.  Limiting  conditions 
of  complete  saturation  can  readily  be  approached  when  only  a  very  few  loop  closures  are  achieved  via 
scanning. 

D.  MJLTIMODALITY  PILOT  MODEL  CHARACTERISTICS 

The  description  of  pilot  dynamics  to  this  point  has  explicitly  assumed  fixed-base  conditions.  Many 
important  piloting  tasks  include  accelerating  and  rotating  flight.  Consequently,  the  effects  of  motion 
cues  on  the  pilot's  dynamic  response  can  be  very  important.  Generally  speaking,  these  can  be  divided  into 
five  categories: 

•  The  motion  may  provide  an  alerting  and  triggering  stimulus  which  activates  an  internal 
command  generator.  This  operates  in  conjunction  with  the  parallel  quasi-optimal  con¬ 
trol  path  of  the  dual-mode  model  (next  chapter).  It  is  perhaps  most  important  for 
recovery  maneuvers. 

•  Motions  indicative  of  status,  such  as  buffet  or  stick  shaking,  provide  alerting  and 
a  consequent  increase  in  neuromuscular  tension.  This  has  the  effect  of  reducing  the 
effective  time  delay  in  the  neuromuscular  system,  thereby  permitting  the  pilot 
to  operate  with  a  higher  gain. 

•  Essentially  steady-state,  moderate  g  levels  can  improve  the  pilot's  dynamic  response 
potential  via  a  neuromuscular  tension  increase  similar  in  effect  to  that  noted  above. 

At  higher  g  levels  the  pilot's  dynamic  capabilities  are  generally  degraded,  e.g., 
gains  are  decreased. 


•  Motion  effects  which  conflict  with  the  srisual  modality  can  cause  illusions  which 
distort  the  pilot's  perception  of  the  state  of  affairs.  These  can  be  so  severe 
as  to  affect  the  pilot's  control  capability  (Refs.  71*  and  115). 

•  The  vehicle  motions  sensed  by  the  pilot  are  used  as  the  basis  for  closed-loop 
control  (Ref.  70). 

The  first  three  of  the  effects  noted  above  are  taken  care  of  in  the  multiloop  model  described  above 
either  as  a  stimulus  or  as  an  adjustment  in  one  of  the  model  parameters.  The  last  two  effects  are  mors 
complicated  and  require  more  extended  discussion.  For  this  purpose,  we  will  draw  heavily  on  Ref.  10, 
which,  in  turn,  relies  extensively  on  data  from  Refs.  1l6  and  117.  To  start,  recognize  that  the  pilot 
contains  neurological  elements  capable  of  sensing  rotary  and  linear  accelerations.  The  primary  neurologi¬ 
cal  elements  involved  as  sensors  are  in  the  vestibular  apparatus,  although  other  sensors  and  pathways  may 
also  be  involved. 

The  rotary  motion  sensor,  usually  associated  with  the  semicircular  canals,  has  a  basic  second-order 
response  to  angular  velocity.  As  shown  in  Fig.  33&,  it  is  in  essence  a  highly  overdamped  angular  accelero¬ 
meter.  Over  the  frequency  range  from  0.2  to  10  rad/sec, 
the  output  signal  is  proportional  to  angular  rate,  so 
the  sensor  can  function  as  a  "rate  gyro."  For  pro¬ 
longed  turning  the  signal  washes  outj  thus,  spurious 
sensations  occur  in  steady  rotations  or  when  the  turn¬ 
ing  motion  stops.  The  threshold  element  shown  in 
Fig.  33a  is  somewhat  of  an  oversimplification  since 
the  threshold  valve,  iuj,  has  both  inter-axis  and 
mter-subject  variability  and  is  also  dependent  on 
the  magnitude  and  time  character  of  the  input  accel¬ 
eration.  It  is  of  the  order  of  1-3  deg  per  sec,  which 
is  sufficiently  low  to  assure  the  presence  in  the 
pilot  control  loops  of  the  rotary  motion  sensors  in 
large  ■\mplitude  recovery  maneuvers,  yet  perhaps  large 
enough  to  make  closed-loop  motion  effects  unimportant 
for  precision  control  near  the  stall. 

Because  the  rotary  motion-sensing  apparatus  gives 
rise  to  p  rate- tyro-like  cue,  the  necessity  for 
visually-generated  attitude  lead  is  reduced.  This 
is  particularly  important  (and  has  been  well  demon¬ 
strated  in  Ref.  70)  when  the  effective  controlled 
element  dynamics  are  K/s^-iike  or  worse.  It  is  also 
worth  mentioning  that  the  primary  effect  of  the  rotary 
motion  feedback  can  be  con/erted  into  an  equivalent 
visual-only  or  fixed-base  situation.  For  instance, 
a  fixed-base  crossover  model  can  be  applied  to  moving- 
base  tracking  by  modifying  the  effective  time  delay 
and  resulting  crossover  frequency.  Specifically,  the 
phase  lag  is  reduced  roughly  equivalent  to  a  time  delay 
reduction  of  0.1  to  0.2  sec,  and  the  magnitude  of  the 
pilot  describing  function,  Yp,  is  increased  to  provide  a  crossover  frequency  increase  of  0.3  to  1.J>  rad/seo. 
This  means  of  accounting  for  the  rotary  motion  cues  is  both  effective  and  extremely  simple  to  apply  in 
analysis . 
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Figure  37.  Angular  and  Linear  Motion 
Sensing  Dynamics  (Ref.  70) 


The  linear  acceleration  sensors  are  ordinarily  associated  with  the  utricles.  Less  is  known  about  the 
dynamics  of  this  sensory  apparatus,  although  the  effective  model  between  linear  acceleration  and  the  pilot's 
output  acting  on  that  linear  acceleration  is  indicated  in  Fig.  3?b.  The  thresholds  in  this  pathway  are 
very  small,  of  the  order  of  0.01  g  or  less,  and  therefore  have  a  negligible  effect  in  most  vehicular  con¬ 
trol  situations.  On  the  other  hand,  the  lag  t-uiie  constant  of  about  2/3  sec  is  more  often  than  not  too 
large  for  the  signal  to  be  useful  for  continuous  closed-loop  control  purposes.  Further,  the  utility  for 
control  purposes  of  acceleration  cues  is  highly  dependent  on  the  pilot's  location.  Consequently,  while 
experiments  have  demonstrated  the  significant  advantages  of  rotary  cues  for  control  purposes,  only  very 
limited  situations  have  done  the  same  for  linear  accelerations  (Ref.  118).  Nonetheless,  every  situation 
where  a  linear  acceleration  feedback  of  this  nature  would  be  advantageous  to  the  pilot  shows  some  evidence 
of  its  use.  The  point  is  that  there  are  few  such  situations. 

To  structure  a  mulsimodality  multiloop  pilot  model,  the  motion  cue  pathways  illustrated  in  Fig.  33  eve 
simply  added  to  the  vector  version  of  Fig.  13. 

In  th  tiscusdon  above,  the  linear  and  angular  motion  sen-^rs  are  treated  as  if  they  are  parallel 
pathways  to  she  visual  modality.  This  is  only  part  of  the-  story,  as  is  shown  in  Fig.  31* •  There,  the 
"nystagmus  crossfeeds"  from  the  canals  and  utricle  to  the  oculomotor  system  will  be  noted.  The  nystagmus 
crossfeeds  produce  involuntary  eye  motions  as  a  function  of  the  excitation  of  the  vestibular  apparatus. 

Such  motions  are  important  in  di sorientation  and  illusions  which  result  from  the  initiation  or  sudden 
cessation  of  large  amplitude  maneuvers.  There  are  also  other  flight  operations  which  have  no  ordinary 
earthbound  equivalent  and  which  give  rise  to  similar  illusory  phenomena.  Several  of  these  are  of  great 
importance  in  producing  conflicting  cues  for  flight  control.  Some  examples  are  given  in  Table  10.  Most 
of  these  rlienomena  have  only  qualitative  connections  with  even  the  most  sophisticated  versions  of  the 
multimodality  multiloop  pilot  model. 
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Figure  J1*-  Pilot /Vehicle  System  for  Head-Fixed,  Pressure-Manipulator,  Attitude  Control  'Tasks 


TABLE  10.  EXAMPLES  OF  POSSIBLE  MOTION  ARTIFACTS  IN  THE  PILOT/VEHICLE  SYSTEM 


SITUATION 

VESTIBULAR  SENSORS 
PRIMARILY  INVOLVED 

REMARKS 

UTRICLE 

(Gu) 

SEMICIRCULAR 
CANALS  (Gso) 

Steady  turn 

X 

X 

Illusion  of  straight  and  level 
flight;  G3C  and  Gu  terms  are 
washed  out. 

Straight  and  level  after 
steady  turn 

X 

Sensation  opposite  of  turning; 
visual  sensation  of  tilt.  Gsc 
and  Gvc  are  dominant  in  creat¬ 
ing  these  sensations. 

Steady  acceleration  in  hori¬ 
zontal  flight,  pushover  from 
steady  climb 

X 

Sensation  of  nose-up  change  in 
attitude. 

Deceleration 

X 

Sensation  of  pitch-down  change 
in  attitude. 

Straight  and  level  alter  high 
angular  path  (>  60  deg/sec) 
aerobatics 

X 

Sensation  of  turning. 

Straight  flight  after  long 
time  high  rate  rolls 

X 

Nystagmus,  blurred  vision, 
reversal  of  background. 

High-frequency  pitching 
rotations,  etc. 

X 

Blurred  vision 

E.  MUITILOOP  PILOT  RATING  CONSIDERATIONS 


A  number  of  schemes  to  obtain  multiloop  pilot  ratings  from  single-loop  considerations  have  been  proposed 
(Ref.  37);  at  present  the  most  promising  (Ref.  119)  is  based  on  excess  control  capacity  concepts.  Assume 
that  the  relationship  between  pilot  rating  and  excess  control  capacity,  \n  -  Xs/>c,  given  by  Fig.  29  is 
applicable  to  each  loop  of  a  multiloop  manual  control  system.  Then,  for  each  such  loop  a  pilot  rating,  R, 
can  be  estimated  using  the  single-loop  correlates  previously  discussed  (Fig.  27)  or,  when  available,  rating 
functionals  like  those  in  Fig.  26  and  an  excess  control  capacity,  ,'n,  assigned  accordingly. 
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Single-axis  capacity,  or  attention,  values  can  be  combined  to  yield  the  combined  axis  value  by  a  multi¬ 
plication  process  (Ref.  120),  i.e.,  the  multiaxis  excess  capacity,  Ann,,  given  by  the  product  of  the 
excess  capacities  for  the  individual  axes: 


m 

II  An^ 


(*»2) 


and  for  R  =■  A  +  BAn  as  a  linear  fit  of  the  Fig.  29  data: 


Hm  =  A  +  BAnjn 
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Combined  ratings  are  always  greater  than  (or  equal  to)  individual  ratings,  since  combined  An's  are  always 
less  than  any  individual  An.  Also,  the  maximum  value  of  R,n  never  exceeds  A,  i.e.,  for  large  Ri  <  A, 
U(R-A)— 0. 

The  lr„ical  value  for  A  is  10. 0;  B  is  determined,  using  the  empirical  data  to  be  equal  to  -8.J.  This 
results  in  a  good  overall  fit  to  all  the  available  multiloop  rating  data  (Ref.  119). 


CHAPTER  VI 


RECONFIGURING  MANUAL  CONTROL  SYSTEMS: 
BY  DESIGN  AND  BY  PILOT  ADAPTATION 


In  Chapter  I  the  stage  was  set  for  subsequent  elaborations  by  presenting  some  of  the  caveats  which 
historically  confronted  efforts  to  model  the  human  pilot  mathematically.  This  pessimism  was  based  on  the 
human  pilot's  manifest  ability  to  learn  new  modes  of  behavior  and  to  adapt  them  to  the  changing  demands  of 
the  vehicle  of  which  he  was  in  control  —  an  ability  which  is  strongly  nonlinear  when  compared  with  the 
responses  of  the  inanimate  components  of  the  control  system.  It  is  pilot  abilities  such  as  the  foregoing 
together  with  multimodal  sensory  perception  and  output  behavior  which  lead  to  a  number  of  possible  control 
loop  structures  for  piloted  vehicles  in  any  given  situation.  Through  higher-order  processes,  such  os  judg¬ 
ment  and  memory,  the  pilot  can  evolve  and  modify  his  performance  criteria,  select  relevant  inputs,  decide 
between  competing  control  loop  structures,  and  optimize  his  fine-grained  behavior  with  respect  to  several 
criteria.  Clearly,  these  are  marvelous  abilities,  and,  in  fact,  the  reasons  for  the  versatility  and 
effectiveness  of  manned  control  systems.  Understanding  this  behavior  can  serve  to  develop  rules  for  the 
pilot/vehicle  system  analyst  for  selecting  the  pilot  model  appropriate  to  a  given  situation  and  c>  1  facili¬ 
tate  the  development  of  system  configurations  which  elicir  and  exploit  this  behavior.  As  we  shalj.  see, 
this  development  is  an  activity  shared  between  the  designer  and  the  pilot.  In  what  follows  we  describe  a 
dual-channel  controller  and  illustrate  the  compatibility  of  such  a  configuration  with  pilot  behavior  for 
both  random-appearing  and  quasi-predictable  inputs.  A  theory  for  the  development  of  control  skills  is 
developed,  and  the  stages  in  this  development  related  to  known  capabilities  of  pilots.  Finally,  the 
behavior  of  the  pilot  in  responding  to  transient  inputs  is  placed  within  this  context  for  skilled  control. 

A.  SOME  PROPERTIES  OF  A  DUAL-CHANNEL  CONTROLLER 

The  well-developed  theory  for  pilot  single-loop  control  behavior  presented  hitherto  has  treated  the 
pilot  as  if  he  were  responding  to  information  arising  from  e(t)  alone.  Precise  input/output  measurements 
were  made  for  this  single  input  channel  mode  of 
behavior,  but  these  measurements  were  not  capable 
of  distinguishing  different  cognitive  operations 
which  the  pilot  might  elect  to  perform  upon  his 
inputs.  In  Chapters  III  and  IV  we  showed  how  a 
model  for  neuromuscular  control  and  actuation 
could  be  developed  from  physiological  data  and 
measurements  on  the  intact,  non-physiologically-  i 

instrumented,  pilot.  We  will  now  present  models  — 

for  the  internal  cognitive  organization  of  the 
pilot's  information  used  for  control. 

Consider  a  situation  in  which  both  the  for¬ 
cing  function,  i(t),  and  the  error,  e(t),  are 
available  for  use  in  control.  Figure  35  pre¬ 
sents  this  situation,  where  Yp^  operates  on  the  Figure  35*  A  Dual -Channel  Controller 

input  channel  and  Ype  on  the  error  channel. 


If  a  systems  analyst  were  presented  the  controller  in  Fig.  35>  as  an  unknown  to  be  identified,  he  would 
proceed  as  discussed  in  Chapter  III  and  compute  describing  functions,  one  of  which  we  will  call  Yp. 

•ini  yc(*Pi+Ype) 

^  =  "T-YcYPi  1  (uu) 


Were  the  controller  responding  to  eft)  alone,  i.e.,  Yp<  =  0,  tnen  Y(?/Yc  =  Ype  =  Yp. 
mode  the  closed-loop,  M/I,  and  error/input  describing  function,  E/t,  are  the  following: 


In  the  dual-channel 


_IL_ 

1  +  Yp 

Yc(YPi  +  *Pe) 
1  +  YcYPe 


1  ~  Vc YPi 
1  +  YcYPe 


Subject  to  the  constraints  imposed  by  stabilization  demands  we  would  like  Yp  to  be  very  large  so  that 
lM|  /  III  1  and  |E|  /  III  0.  This  could  j  accomplished  if  the  block  Yp^  were  "smart"  enough  to  adjust 
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Itself  so  that  Yp^Yg  i  1 .  Increasing  Yp  to  achieve  the  same  end  would  be  both  potentially  destabilizing 
and  an  inefficient  use  of  information.  For  YP^  to  be  calibrated  to  Yc,  functional  Knowledge  of  Yc  would 
have  to  be  in  YPj,'s  "memory."  This  knowledge  could  be  wired  in  or  it  could  evolve  with  experience. 

If  the  systems  analyst  were  to  suspect  that  he  was  measuring  a  controller  capable  of  acting  in  a  dual¬ 
channel  mode,  he  could  make  predictions  about  the  form  of  Y[-,.  He  would  expect  the  compensatory  Ype  to  be 
present  and  active  if  the  controlled  element  were  unstable  or  if  the  forcing  function  were  such  as  inher¬ 
ently  to  induce  a  remnant  (e.g.,  random-appearing).  In  these  events  the  Ype  loop  would  have  to  be  closed. 
On  the  other  hand,  if  the  controller  had  a  functional  knowledge  of  Yc  so  as  to  compensate  for  system  lags, 
and  knew  i(t)  sufficiently  well  to  enable  ''prediction"  of  its  future  course,  Ype  could  be  negligible  for 
varying  intervals  of  time  and  YpAeil.  A  feedforward  of  this  sort  could  exhibit  prodigious  performance 
by  acting  on  estimates  of  the  future  input  to  compensate  for  lags  in  itself  and  the  plant  under  control. 

The  hypothetical  behavior  described  has  parallels  in  human  pilot  operations.  There  are  various  means 
by  which  this  pilot  dual-channel  multimode  behavior  can  be  facilitated.  In  the  following  sections  we  will 
review  these  means. 

B.  THE  PURSUIT  DISPLAY 


The  most  direct  means  for  helping  the  pilot  become 
a  dual-channel  controller  —  bearing  in  mind  that  the 
dual  channels  are  cognitive  organizations  of  input 
information  which  can  be  facilitated  and  supported  but 
not  compelled  by  engineering  design  —  is  to  present 
the  necessary  inputs  to  him.  The  pursuit  display  does 
this  by  presenting  the  pilot  a  moving  target  which  he 
pursues  in  an  effort  to  capture  it  w.th  a  cursor  he 
commands  through  the  vehicle  dynamics.  Figure  *6  is 
a  schematic  of  this  display.  Three  possible  inputs 
are  illustrated  in  Fig.  ' 6 ,  but  only  two  are  indepen¬ 
dent;  so  that  the  dual-channel  model  does  not  result  in 
any  loss  of  generality.  Although  the  body  of  research 
on  pursuit  display  tracking  for  which  dynamic  measure¬ 
ments  exist  is  limited,  there  are  measurements  from 
which  we  can  determine  whether  the  pilot  is  capable  of 
generating  operations  which  enhance  the  performance  of 
a  dual-channel  controller. 

1 .  Random-Appearing  Inputs 

In  order  to  measure  the  elements  in  Fig.  A;  one  must  contrive  situations  which  a  priori  make  one  or  the 
other  block  dominant  (Refs.  21  and  121 ),  make  assumptions  about  the  form  of  Ypi  or  Ype  (Ref.  122),  or  intro¬ 
duce  an  additional  uncorrelated  input  so  as  to  be  able  to  solve  for  two  unknown  blocks  (Refs.  C]  and  70). 

The  first  approach  has  been  used  for  data  with  a  pure-gain  controlled  element  for  ^  •-  rad/sec  inputs; 
with  the  assumption  that  since  the  task  was  quite  easy  and  the  error  small,  signals  fed  back  through  Ype 
were  negligible.  Under  these  circu..istances  measurements  of  11, 1,  the  closed- loop  transfer  function,  are 
identical  with  measurements  of  the  feedforward,  Ypj, ,  and  the  pure-gain  plant.  On  Fig.  77a,  the  fidelity 
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Figure  J7a.  Ypi,  Inferred,  for  Yc  -  Kc 


Figure  77b.  Indirectly  Measured  YPi 
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with  which  Ypi  approximates  Y0  =  K  (0  dB,  0  deg  phase)  over  a  large  range  is  remarkable.  The  data  were 
generated  in  different  laboratories,  but  were  .malyzed  to  the  same  end  —  to  determine  whether  Yp^  could 
be  fit  by  an  optimal  linear  predictor  in  the  minimum  mean-square  sense  which  operated  only  upon  the  posi¬ 
tion  and  velocity  of  the  input  signal.  It  was  demonstrated  (Refs.  21  and  121)  that: 


Ypj.  =  (a  +  Wu,)e-*P*B 


(“7) 


where  the  prediction  interval,  Tp>  was  0.2  sec.  When  a/(-  was  less  than  about  5,  Yp*  1  Yc-1 .  As  the  task 
increased  in  difficulty,  the  approximations  which  led  M/I  =  no  longer  held,  since  the  error  loop  became 
significant. 


On  the  assumption  that  the  feedback,  Ype,  in  a  dual  channel  does  not  differ  significantly  from  the  Yp 
of  a  single-channel  compensatory  mode  for  the  same  subject  with  similar  Yc  and  .ip,  the  feedforward  for 
the  pursuit  case  can  be  determined  as: 


‘P- 


Vll£p 

Yc(1  +  Yt,) 
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Here,  Yp  is  measured  in  the  pursuit  situation  and  Yp  for  the  corresponding  compensatory  case. 

This  calculation  has  been  made  for  several  controlled  elements.  On  Fig.  37b,  which  was  generated  using 
an  unstable  second-order  system,  we  note  a  close  correpondence  between  predicted  and  computed  Ypp  magni¬ 
tudes  (Ref.  122).  Certain  data  have  indicated  that  when  the  pilot-generated  feedforward  cannot  effect  an 
adequate  prediction  of  i(t)  so  as  to  allow  a  lag-free  inversion  of  Yc,  a  delay  is  present,  i.e.: 


YpiYc 


e-d<p(a>) 


(1*9) 


This  effect  accounts  for  the  relatively  poor  rhas  •  match  on  Fig.  ?7b. 

The  substitution  of  Yj  for  Yp@  in  the  calculations  of  Fig.  '7b  was  justified  by  measurements  using  an 
additional  uncorrelated  iruut  —  the  third  technique  (Ref.  07).  That  is,  a  disturbance,  d(t),  added  to 
c(t)  in  Fig.  T5  served  to  specify: 
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and  an  open-loop  describing  function  for  the  dual-channel  controller,  which  can  be  designated  Yjjp,  was 
specified  by  i(t): 
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Equations  50  and  5'  were  then  solved  for  Y:.,„  and  Y;,^.  The  pilot's  pursuit  behavior  was  scarcely  affected 
by  the  additional  input  (Ref.  67).  .?Pe  exhibited  a  somewhat  lower  o<;  and  greater  pnase  margin  than  Yp. 
This  could  account  for  some  of  the  discrepancies  in  Fig.  37b. 

2.  Predictable  Inputs 

The  most  effective  operation  of  the  feedforward  requires  that  i(t)  be  predictable.  As  predictability 
deteriorates,  or  involves  more  internal  data  processing,  the  lag  in  Eq.  1*9  becomes  irreducible.  Thus, 
although  Fig.  37a  indicates  only  minor  phase  lags  at  frequencies  less  than  1  rad/sec,  it  is  clear  from 
Eq.  47  that  the  prediction  interval,  Tp  =  0.2  sec,  must  generate  phase  lags  of  increasing  impact  at  higher 
frequencies.  The  data  at  higher  frequencies  bear  this  out.  The  lags  inherent  in  generating  a  prediction 
of  the  input  signal  are  analogous  to  the  dynamic  costs  of  low-frequency  lead  generation  illustrated  in 

Fig.  15. 


There  exists  in  practice  a  large  number  of  examples  of  quasi-predictable  inputs  which  the  pilot 
might  be  able  to  reproduce  with  sufficient  fidelity  so  as  to  minimize  or  perhaps  eliminate  these  time 
costs  for  predictions.  Such  inputs  arise  in  a  variety  of  piloting  tasks.  Examples  are:  following  the 
optical  landing  beam  of  an  aircraft  carrier  plunging  through  deep  ocean  swells;  compensating  for  low- 
frequency,  lightly-damped  vehicle  modes;  fighting  pilot-induced  oscillations;  terrain-following  flight 
over  rolling  countryside,  etc.  An  example  of  such  quasi-predictable  inputs  are  the  actual  ship  motion 
time  histories  shown  in  Fig.  38*  A  carrier  pilot  would  like  to  match  perfectly  the  corresponding  deck 
motions  with  his  aircraft  landing  gear  if  this  could  be  accomplished  through  adequate  vehicle  response 
bandwidth  and  suitable  displays  and  controls. 

There  are  many  examples  of  how  well  the  pilot  can  match  similar  quasi-predictable  motion.  Pure  sine- 
waves  or  simple  combinations  of  a  few  sinewaves  are  subjectively  and  actually  clearly  predictable  and,  as 
such,  thus  provide  an  upper  limit  for  the  pilot's  feedforward  operation.  A  striking  example  of  how  well 
the  pilot  can  generate  a  prediction  of  i(t)  is  presented  in  Fig.  39,  where  the  human  continued  to  produce 
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J —  "tights  Out” 


Figure  39*  Demonstration  of  Pattern  Generation 
(Adopted  from  Ref.  12h) 


a  highly  similar,  though  phase-shifted,  version  of 
the  original  harmonic  signal  when  he  could  no  longer 
see  either  i(t),  e(t),  or  m(t)  (Ref.  12h).  Further 
evidence  of  such  internal  signal  generation  is  found 
in  Reft  123,  1<5,  and  126). 

Data  from  the  various  previous  studies  of  pure 
sinewave  tracking  (Refs.  37,  75,  and  127)  can  be 
coalesced  and  presented  in  describing  function  form 
to  illustrate  the  pilot's  ability  in  this  extreme 
case  of  signal  predictability  (see  Fig.  to).  Limit¬ 
ing  conditions  exist  at  the  low-frequency  end  where 
the  signal's  rate  of  change  is  so  low  as  to  require 
the  pilot  to  use  the  YPe  loop  for  continual  correc¬ 
tion,  and  above  ~20  rad/sec  where  neuromuscular 
limitations  come  into  play.  This  range  of  high 
fidelity  in  frequency  generation  (note  the  essen¬ 
tially  zero  phase  lags)  is  corroborated  by  a  variety 
of  pursuit  lapping  output  studies  in  response  to 
either  auditory  or  visual  inputs  (Ref.  125).  It  was 
found  that  the  standard  deviation  of  the  time  inter¬ 
vals  between  i(t)  and  m(t),  analogous  to  phase  varia¬ 
tions,  ranged  from  It  to  6  percent  of  the  command 
tapping  rate  up  to  a  commanded  period  of  about 
300  ms. 

The  foregoing  can  be  summarized  in  the  following 
rales  for  pilot  tracking  of  sinewave  inputs. 

•  In  the  region  below  about  0.3  Hz,  it  is  advantageous  to  use  a  Ypj,  =  1/Ye  feedforward 
operation,  with  Ype  present  for  vernier  corrections  and  remnant  suppression.  In 
this  region  rhythms  are  hard  to  reproduce  accurately.  At  very  low  frequencies 
(below  0.1  Hz)  the  operator  operates  to  reduce  the  error,  i.e.,  uses  the  Ype  com¬ 
pensatory  block  solely. 

•  From  0.5  to  about  1.0  Hz,  rhythm  detection  permits  activation  of  an  internal  "pattern 
generator"  block,  probably  using  proprioceptively  perceived  patterns.  The  feedforward 
block,  Ypp,  aids  in  producing  the  "pattern  generator"  response.  The  compensatory 
loop,  Ype,  can  still  be  closed  tightly  to  phase  lock  the  output  and  suppress  remnant. 

•  From  about  1  to  just  under  2  Hz,  the  pattern  generator  loop  is  active,  but  as  the 
frequency  exceeds  the  operator's  compensatory  crossover  limit  he  must  open  the  com¬ 
pensatory  loop  to  avoid  exciting  undue  overshoot  errors.  In  this  region  control 
approaches  ultimate  predictive  behavior  with  Ype  acting  intermittently  to  prevent 
frequency  drift,  but  unable  to  prevent  random  phase  and  amplitude  errors. 

•  From  about  2  Hz  to  the  neuromuscular  response  limit  of  9  to  10  Hz,  the  operator  only 
uses  an  internal  pattern  generator  loop  to  roughly  approximate  the  displayed  ampli¬ 
tude  and  frequency,  usually  undershooting  both. 

Inputs  can  be  scaled  for  suojective  predictability  along  the  dimension  of  signal  shape  or  waveform 
features  and  their  "coherence"  or  waveform  time  variations;  such  scaling,  however,  is  not  highly  developed. 
This  scaling  of  expected  inputs  would  be  helpful  to  the  analyst  and  to  the  designer  who  might  be  able  to 
either  specify  or  induce  YPj_  forms  and  control  configurations  appropriate  to  classes  of  signals.  Because 
of  their  apparent  similarity  to  such  practical  inputs  as  are  illustrated  on  Fig.  38,  bandpassed  random 
noise  provides  an  attractive  idealization  of  an  input  of  varying  predictability.  Examining  its  properties 
could  assist  the  analyst  in  predicting  the  performance  of  control  systems.  It  has  been  found  (Ref.  123) 
that  an  effective  rating  for  subjective  predictability  can  be  computed  from  the  bandpass  (half-power  width), 


Figure  20.  Frequency  Response  Data  for  a  Single 
Sinewave  Input  at  Various  Frequencies 
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A,  and  center  frequency,  f0,  of  the  narrowband  filter  through  which  the  randan  process  is  passed*  The 
properties  of  the  ratio,  v  =  A/fo,  determine  subjective  predictability  as  follows: 


Predictable 

v  <  0.3 

Non-predietable 

v  >  0.8 

For  Inputs  with  v  <  0.3,  there  is  evidence  that  the  pilot  generates  a  feedforward  acting  like  a 
Kalman  filter  predictor  for  Y0  =  K0: 


Ypi  = 


sin  itv 
kv 


e-Mo 


(52) 


where  l/T0  is  the  center  frequency  of  the  narrowband  filter.  The  prediction  interval,  To,  is  the  cost 
imposed  on  the  pilot  by  conditions  under  which  he  cannot  internally  generate  the  appropriate  prediction 
but  must  operate  upon  the  signal  instead. 

3.  Are  Pursuit  Displays  Desirable? 

Whether  a  designer  should  attempt  to  facilitate  the  generation  by  the  pilot  of  a  feedforward  loop  by 
presenting  him  with  a  pursuit  display  is  not  an  easily  answered  question.  Ambiguities  arise  because  the 
physical  display  does  not  guarantee  the  extent  to  which  the  pilot  will  utilise  tne  added  information 
channels .  Table  1 1  lists  the  available  data  comparing  pursuit  and  compensatory  displays  for  which  describ¬ 
ing  functions  have  been  measured.  Such  comparison  data  which  are  in  the  form  of  performance  measures  alone 
have  not  been  listed,  because  there  is  no  way  by  which  the  selected  pilot  control  structure  can  be  adduced. 
For  example,  in  Ref.  78  it  was  found  that  the  influence  of  the  Yp^  loop  was  negligible  in  the  control  of 
a  Kc/s,  and  that  Ype  predominated.  In  other  words,  there  was  very  little  cognitive  difference  between  a 
physical  compensatory  and  physical  pursuit  display  under  the  measurement  conditions.  The  performance  data 
confirmed  this  indifference,  but  without  the  dynamic  data  these  performance  measures  would  have  been 
uninterpretable.  Comparisons  based  on  relative  mean-square  error  must  be  interpreted  with  care,  because 
the  pilot  may  respond  to  additional  control  criteria.  With  this  in  mind  we  can  review  comparisons  based 
on  e2/i2  from  Refs.  21,  67,  78,  and  12?.  The  most  definitive  findings  are  for  Yc  =  Kc  where  pursuit  dis¬ 
plays  are  superior  to  compensatory  displays  in  the  range  from  1  to  10  rad/sec.  Outside  this  range,  where 
compensatory  tracking  is  at  all  feasible  the  compensatory  display  is  to  be  preferred.  For  Kc/s  over  the 
range  of  measurements  there  is  not  much  difference  based  on  performance  measures  alone.  What  little  data 
exists  for  Kq/s2  indicate  a  slight  preference  for  the  pursuit  display.  There  is,  however,  evidence  that 
the  pilot  requires  less  control  movement  for  Kc/s2  in  the  pursuit  than  compensatory  configuration,  and  he 
appears  to  like  this  condition. 

Other  bases  for  comparison  between  the  two  displays  ere  presented  in  Table  12.  As  it  can  be  seen, 
there  are  many  desiderata  in  selecting  a  control  display  configuration  in  addition  to  performance  measures. 
Perhaps  the  most  important  for  aircraft  applications  is  the  indication  of  status  information  as  well  as 
command  and  error  present  in  the  pursuit  display. 


TABLE  11 

SUMMARY  OF  HUMAN  DYNAMICS  DATA  K)R  EXPLICIT  PURSUIT  DISPLAYS 


INVESTIGATOR 

NUMBER  OF 
TYPES  OF 
RANDCM- 
APPEARBK3 
INPUTS 

INPUT  REMARKS 

ADDITIONAL 

DISTURBANCE 

INPUT 

NUMBER  OF  TYPES  OF 
CONTROLLED  ElfMENTS 

DIRECTLY 
COMPARABLE 
WITH  EXPLICIT 
COMPENSATORY 
DATA 

REMARKS; 

NUMBER  OF  SUBJECTS 

Elkind 
(Ref.  21) 

20 

1"  rms;  (311m  of  sine- 
waves) 

No 

(l),  Yc  =  Kc,  pencll-liko 
manipulator 

Yes 

5 

Waslcko, 
UcRuer,  and 
Magda leno 
(Ref.  127) 

5 

O.25"  to  1"  rms;  all 

3  inputs  were  present 
to  3  of  the  controlled 
elements  (sum  of  sine- 
waves) 

No 

(7),  Yc  *  Kc,  Kc/s,  Kc/s2 
Kc/s(s  -  0-5) 

Kc/s(s  -  1.3) 
Ko/((=’0.25)/(s+5.0)2) 
Lateral  side  stick 

Yes 

\ 

Allen 
and  Jex 
(Ref.  67) 

1 

O.h"  ms,  J.2  rad/sec, 
dominant  frequency 
(sum  of  3inewaves) 

Vos 

(2),  Kc/s,  Kc/s2 

Lateral  side  stick 

Yes 

It 

Obtained  comments 

Reid 

(Ref.  78) 

!» 

Filtered  noise, 

«ds  =  1  -  5,  1.'‘,  2.8  deg, 
U.J  rad/soc 
}  Inputs  had  high- 
frequency  shelves 

Yes 

(2),  K0,  Kc/s 

Low  inertia 

Restraint  free  stick 

Yes 

8 

Disturbance 
magnitude: input 
magnitude 

Ware 

(Ref,  121) 

It 

Filtered  noise 
<i>iQ  «  1.2,  \.U,  2.}, 
and  2.8  rad/sec 

No 

(l),  Yc  3  Ks 

Force  stick 

No 

3 

53 


TABLE  12 

COMPENSATORY  VERSUS  FJRSUIT  DISPLAYS 


COMPENSATORY 


PURSUIT 


Advantages 


1 .  Simple  to  build  Instruments,  and  to  simu¬ 
late.  Most  reliable. 

2.  Expanded  scale  possible;  easy  to  detect 
email  errors. 

3.  Simple  to  Interpret;  only  one  fora  of 
action  and/or  Yp  Is  required. 

U.  Can  be  easily  optimized  and  equalized 
uslrg  pllot/vehlcle  systems  analysis. 

5.  Cnly  type  feasible  for  pure  regulation 
tasks  (e.g.,  suppress  gust  distur¬ 
bances). 


Advantages 

1.  Shows  more  of  available  lnforaatlon;  permits  operator  to  adopt 
separate  control  criteria  for  tracking  versus  regulation; 
enhances  the  ability  to  control  quosl-prcdlctable  Inputs 

2.  Operator  can  develop  a  feedforward  path  (Yp^)  operating  directly 
on  the  Input  to  minimize  the  closed-loop  errors.  Results  in 
reduced  tcff  and  higher  overall  system  bandwidth,  with  lower 
gain  and  fewer  stability  problems  in  the  compensatory  loop. 

3.  Improved  "conforaabllity,"  i.e.,  a  proper  pursuit  Instrument  is 
a  closer  analog  to  the  visual  field,  permitting  easier  VfR  and 
im  transitions. 

U.  Less  control  cx>tion  required  for  higher-order  controlled 
elements 


Disadvantages 

1 .  Cannot  separate  disturbances  from  com¬ 
mands;  leads  to  conservative  closed-loop 
stability  criteria  and  larger  tracking 
errors  for  low-frequency  Inputs. 

2.  Error  may  not  agree  with  secondary  cues 
(e.g.,  normal  accelerations  versus  e, 
while  tracking  a  maneuvering  target). 

3.  Command  patterns  may  be  masked  by 
remnant-induced  errors,  thereby  Impeding 
improvements  due  to  more  effective  signal 
prediction. 


Disadvantages 


1.  Sealing  set  by  largest  Input  command;  may  result  In  errors  being 
too  small  to  use  effectively. 

2.  Harder  to  put  inputs  on  Instruments,  e.g.,  maneuvering  target 
position  In  space  Is  hard  to  derive  except  from  direct  visual 
field. 

3.  No  Improvement  for  regulation  against  disturbances  {for  zero 
Inputs  eases). 

U.  More  difficult  to  Interpret;  Ypt  and  Ype  Involved. 

5.  More  complex  to  build;  less  reliable. 


C.  SUCCESSIVE  ORGANIZATIONS  OF  PERCEPTION 

In  the  preceding  discussion  of  the  dual-input  configuration  we  noted  instances  in  which  the  pilot 
responded  to  the  input  signal  as  if  it  were  a  trigger  to  release  a  stored  or  self-generated  program.  Such 
instances  occurred  for  the  perfectly  predictable  sinewave  or  simple  compound  sinewave  inputs.  The  pilot 
acts  as  if  the  only  signals  he  requires,  once  he  is  fully  familiar  with  the  input  and  controlled  element 
and  knows  when  to  be^in,  are  those  generated  by  the  proprioceptors  involved  in  his  control  movements.  The 
level  of  skill  is  akin  to  that  exhibited  by  a  trained  motorist  applying  controlled  foot  pressure  to  his 
vehicle's  brakes  on  seeing  a  red  signal  light  ahead,  or  overtaking  and  passing  another  car.  In  these 
examples  the  control  skills  constitute  repertories  of  action  which  are  called  up  at  the  appropriate  moment 
and  released.  In  this  manner  a  variety  of  different  criteria  can  be  realized.  The  controller  may  elect  a 
minimum  effort  maneuver,  a  minimum  time  maneuver,  a  minimum  error  maneuver,  or  he  may  apply  his  own  rela¬ 
tive  weights  to  the  foregoing  and  other  considerations.  As  with  any  predictive  control  technique,  this 
type  of  behavior  requires  the  rest  of  the  world  to  maintain  its  position  during  the  maneuver.  If,  for 
example,  after  pulling  out  of  the  lane  and  drawing  alongside  of  the  car  to  be  passed  in  a  preprogrammed 
maneuver  the  passing  vehicle  is  confronted  by  a  change  in  lateral  position  or  a  forward  acceleration  by  the 
vehicle  being  passed,  deep  trouble  could  result.  It  is  the  need  to  watch  the  rest  of  the  world  —  or  main¬ 
tain  an  intermittent  Ype  —  which  can  contribute  apprehension  and  add  to  the  workload  of  otherwise  simple 
tasks.  If  the  Yp0  loop  is  closed  too  tightly  about  lateral  position,  for  example,  the  skilled  passing 
maneuver  is  almost  impossible  to  achieve.  On  the  other  hand,  a  monitoring  mode  is  essential  to  prevent 
disaster  —  however  unlikely  the  contributing  events  may  be.  The  kind  of  intermittent  open-loop  behavior 
which  we  have  been  describing  is  called  precognitive,  and  it  represents  the  highest  level  in  a  progressive 
development  of  motor  skill  by  stages  of  information  exploitation  which  are  Successive  Organizations  of 
Perception  (GOP).  The  theory  has  been  expounded  elsewhere  (Refs.  57,  7b,  125,  128,  and  129)  and  will  only 
be  summarized  here.  Figure  Ui,  where  remnants  and  disturbances  have  been  omitted,  will  be  useful  in  this 
regard . 

•  The  initial  phase  is  the  familiar  compensatory  display.  The  pilot  is  capable  of 
acting  only  in  response  to  the  error  signal.  Some  analogies  to  this  level  of 
skill  are:  a  blind  man  feeling  his  way  with  a  cane  across  stepping  stones;  a 
motorist  proceeding  in  a  dense  fog  by  hugging  the  cartway  guideline;  or  a  very 
drunken  motorist,  having  generated  a  dense  internal  fog,  cautiously  proceeding 
by  edging  along  the  curb.  The  last  is  an  example  of  a  regression  from  a  higher 
to  a  lower  level  of  skill. 

•  The  second  phase  is  an  elaboration  upon  the  dual-channel  configuration  which  we 
have  discussed  in  detail.  The  additional  channels  are  Yp_,  the  proprioceptive 
feedback  which  is  particularly  useful  in  committing  simple  periodic  inputs  to  a 
memory  of  motor  responses,  and  Yp,,,  an  operator  on  the  system  output.  Although 
a  controller  can  conceivably  develop  all  three  of  the  transfer  characteristics, 

Ype,  Ypj^,  Yp^j  any  two  are  adequate  to  describe  the  control  process.  An  expli¬ 
cit  pursuit  display  enables  the  pilot  to  more  readily  advance  to  the  final 
phase  of  skill. 
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•  Th  •  final  or  p-  ^cognitive  phase  can  be 
conveniently,  but  not  usefully  in  terms 
of  operationally  defined  measurements, 
conceived  as  a  series  of  decision  rule 
algorithms  and  stored  programs.  Phase 
three  is  the  graceful  or  adept  maneuver 
mode  of  behavior  —  skilled  gymnastic 
evolutions,  batting  a  ball,  etc.  The 
components  are  a  stored  repertory  and  a 
knowledge  of  which  cues  to  use  in  releas¬ 
ing  this  repertory  in  a  timely  fashion. 
The  open-loop  sequences  may  be  short,  as 
a  timely  bang-bang  control,  or  long  as 
in  generating  a  synchronous  sincwave. 

Hitherto  we  have  emphasized  the  pilot's  predic¬ 
tive  ability  as  the  essential  element  in  achieving 
the  state  of  precognitive  grace  —  for  however 
briefly.  It  should  be  clear  from  the  description 
of  precognitive  control  that  an  additional  element, 
the  dynamics  of  programmed  responses  in  the  pilot's 
repertory,  must  also  be  considered.  These  responses 
are  those  which  may  be  called  up  while  the  pilot  is 
engaged  in  some  form  of  steady- state  behavior  such 
as  continuous  tracking  or  daydreaming.  It  has  been 
found  (Ref.  71*)  that  such  combined  transient  and 
stationary  situations  can  adequately  be  modeled  with 
a  single  input  channel  dual-mode  structure  such  as 
that  shown  in  the  simplified  diagram  of  Fig.  42.  To 
represent  a  multiloop  situation  the  signals  shown  in 
this  block  diagram  could  be  considered  as  vector 
quantities.  The  quasi-linear  steady-state  path  is 
the  one  used  for  tracking  random  inputs  or  distur¬ 
bances.  It  is  the  same  model  described  in  great 
detail  in  previous  chapters.  The  feedforward  ele¬ 
ment  operates  on  the  transient  inputs  provided 
either  from  the  system,  such  as  a  randan-occurring 
step  sequence,  or  induced  by  an  internal  command 
repertory  from  the  pilot. 

The  nature  of  the  switching  and  the  feedfor¬ 
ward  element  is,  in  the  simplest  terms,  such  as  to 
divide  the  total  pilot  behavior  into  temporal 
pluses,  each  having  a  different  system  organiza¬ 
tion.  As  an  elementary  example,  consider  the  typi¬ 
cal  system  step  response  shown  in  Fig.  ij.  In  the 
first  or  time  delay  phase,  nothing  happens.  This 
is  followed  by  a  rapid  response  phase  and  finally 
by  a  more  or  less  oscillatory  error  reduction 
phase.  The  feedforward  element  provides  the  major 
component  of  the  pilot's  output,  c,  during  the 
rapid  response  phase.  The  quasi-linaar  steady- 
state  model  is  predominant  in  the  error  reduction 
phase  (and  in  general  during  the  time  delay  phase 
if  the  system  is  continuously  excited  by  random 
disturbances  or  inputs  other  than  the  step  assumed 
here ) . 

The  output  of  the  feedforward  element  for  a 
skilled  pilot  is  peculiar  to  each  controlled  ele¬ 
ment  lorm.  Some  appreciation  for  this  variation 
with  controlled  element  dynamics  can  be  gained  from 
the  responses  shown  in  Fig.  44.  These  are  sample 
large  maneuver  responses  to  step  commands  for  effec¬ 
tive  controlled  elements  given  by  Kc,  K<;/s,  Kc/s2, 
and  Kc/s5,  proceeding  from  left  to  right,  respec¬ 
tively.  The  most  important  aspect  immediately 
apparent  from  these  data  is  the  pulse-like  bang- 
bang  nature  of  the  stick  deflection  control  move¬ 
ments.  In  fact,  the  pilot's  control,  c,  is  a 
remarkably  good  approximation  to  the  controller 
properties  of  the  nth  order  (n  =0  to  j),  single¬ 
input.  single-output,  time-optimal  control  system 
with  |c(t)|  <  M.  Here,  the  scaler,  M,  may  repre¬ 
sent  either  a  physical  limit  on  the  stick  deflec¬ 
tion  or,  more  likely  in  the  piloted  case,  an 
implicit  restraint  imposed  by  the  pilot  for  the 
given  situation.  Ideal  time-optimal  traces  for 
comparison  with  the  actual  piloted  conditions  are 
shown  in  Fig. 
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c)  Final  Phase  :  PRECOGNITIVE  ( Open-Loop ) 


Figure  41.  The  Three  Main  Phases  in  the  Successive 
Organization  of  Perception  (SOP) 


Figure  42.  Structure  of  the  Dual-Mode  Model 


Figure  4J.  Typical  System  Step  Response 


Figure  U.  Sample  Large  Maneuver  Response  Data  for  Several 


Idealized  Controlled  Elements 
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Figure  L5.  Ideal  Time-Optimal  Response  Characteristics 


A  dual-mode  pilot  maneuvering  model  which  serves  to  explain  all  the  available  data  is  presented  in 
Fig.  tf.  The  nonlinear  error  sensing  blocrts  automatically  route  the  error  signal  through  the  appropriate 
channel  based  upon  whether  the  error  is  greater  or  less  than  some  threshold  magnituae  of  error,  ep.  The  con¬ 
trol  logic  for  each  different  controlled  element  and  as  a  function  of  the  error  state,  e[e  =  col  (e,  e, 
is  given  in  Table  1?  for  time-optimal  response.  Note  that  M,  the  constraint  on  the  control  input,  is  some 
function  of  the  step  input  height,  controlled  element  gain,  and  its  order.  The  decision  logic  model  behaves 
line  a  function  switch  (K5V?)  and  accounts  for  the  initial  increase  in  the  time  delay  (beyond  that  due  to 
quasi-llnear  tracking)  in  response  to  a  step  input. 
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Figure  1*6.  Dual-Mode  Controller  Model 


TABLE  1*, 


CONTROL  LOGIC  FOR  VARIOUS  CONTROLLED  ELEMENTS 


This  same  model  is  appropriate  when  the  step  function  or  other  transient  input  is  initiated  "inside" 
the  pilot.  In  this  instance,  of  cou.se,  the  internal  stimulus  is  indiscernible;  only  the  response  after 
any  time  delay  can  be  seen.  Yet  the  skilled  pilot's  output  in  programmed  maneuvers  exhibits  at  least  the 
timing  of  this  dual-mode,  quasi- time-optimal  character. 

The  decision  time  delay  called  out  in  Fig.  <6  includes  the  decision  interval  during  which  the  pilot 
determines  the  parameters  of  his  preprogrammed  response.  This  interval  may  contain  wild  fluctuations 
under  conditions  of  suddenly  failing  dampers,  autopilots,  or  other  controlled  element  transitions.  An 
extensive  review  of  findings  under  such  circumstances  has  been  carried  out,  (Ref.  150)-  Simple  transitions, 
such  as  ±  changes  in  gain  for  Yc  =  Kc  are  detected,  and  an  appropriate  response  begun  within  about  O.h  sec 


I 


,1 


f>7 

and  the  error  stabilized  at  its  pretransition  value  in  about  1  sec.  For  transitions  from  Ye  =  Kc  to 
Yc  =  Kc/s2,  detection  and  correction  occur  within  about  1  sec,  but  there  may  be  an  extensive  period  devoted 
to  compensating  for  the  dynamics  of  the  new  controlled  element.  Oscillations  from  6-20  sec  cxe  not  uncom¬ 
mon.  These  examples  serve  to  leave  the  reader  with  a  caveat  about  the  desirability  of  operating  at  the 
highest  level  of  shill  when  controlling  an  unstable  plant.  The  pilot  should  monitor  his  control  with  Ype 
at  periods  such  that  any  potentially  disastrous  divergences  can  be  apprehended  and  controller  This  imposes 
a  practical  limit  on  the  duration  that  the  open-loop  mode  in  Fig.  tic  can  be  maintained. 
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CHAPTER  VII 

BRIEF  REVIEW  OF  PHOT  M3DSL  APPLICAIIOIS 


Since  its  first  tentative  and  rudimentary  applications  in  the  early  195O3,  the  use  of  pilot  models  in 
pilot/vehicle  system  analysis  has  grown  exponentially  to  the  point  where  it  is  now  commonly  applied  to  con¬ 
sider  a  very  vide  range  of  problems.  Some  appreciation  of  the  scope  of  application  efforts  is  revealed  in 
the  bibliography  which  constitutes  the  second  half  of  this  chapter-  There,  the  explications  are  classified 
under  five  headings: 

1.  Filot/aircraft  system  single-  and  multi-loop  regulation  and  control. 

2.  Pilot/space-vehicle  regulation  and  control. 

3.  Limits  of  piloted  control. 

U.  Instabilities  of  pilot/aircraft  systems. 

5.  Flight  control  system  display  design. 

This  classification  selection  resulted  in  a  minimum  number  of  repetitive  entries,  although  a  few  still  exist. 
To  make  the  bibliography  as  comprehensive  as  possible  for  the  benefit,  of  specialists,  the  references  pre¬ 
sented  are  ordinarily  the  original  documents.  These,  more  often  than  not,  are  limited  distribution  but 
usually  available  government  or  company  reports.  For  some  of  the  entries  condensed  versions  in  the  form  of 
archival  journal  articles  are  available. 

The  analytical  techniques  used  to  exercise  existing  pilot  models  in  these  applications  can  be  classified 
under  three  headings.  These  are: 

1 .  Conventional  Flight  Control  System  Analysis  —  The  vast  majority  of  application  results 
have  been  accomplished  using  conventional  feedback  system  analysis  techniques  as  special¬ 
ized  for  aircraft  flight  control  purposes.  The  techniques  and  procedures  used  include 
conventional  Bode  and  root  locus  analysis,  coupled  with  Bode  root  loci,  pole-zero  sensi¬ 
tivity,  and  multiloop  analysis  procedures  lumped  under  the  general  heading  of  Unified 
Servo  Analysis  Methods  (USAM) .  A  comprehensive  summary  of  these  techniques  with  many 
application  examples  for  automatic  control  systems  is  provided  in  Ref.  102.  With  these 
methods  pilot/vehicle  analysis  is  accomplished  by  substituting  pilot  models  wherever 
controller  properties  are  required  and  taking  into  account  pilot  ratings  as  well  as 
dynamic  and  system  performance  characteristics. 

2.  Parameter  Optimization  Techniques  —  In  these  procedures,  the  form  of  the  pilot  model  is 
assumed  a  priori,  hopefully  to  one  pertinent  to  the  control  tasks  being  considered  (e.g., 
by  use  cf  the  pilot  model  forms  of  Chapters  IV  and  V  adjusted  appropriate  to  the  controlled 
elements  of  interest).  The  pilot  parameters  are  then  adjusted  via  a  parameter  optimiza¬ 
tion  scheme  to  minimize  some  performance  index.  In  some  cases  this  performance  index  can 
be  pilot  rating  itself  with  the  criterion  then  being  to  adjust  pilot  parameters  such  that 
pilot  rating  is  minimized  (Refs.  96-1OO’). 

3 .  Quadratic  Optimal  Control  —  In  these  formulations,  conventional  optimal  control  theory 
is  modified  to  permit  a  pure  time  delay  and  observation  noise  (remnant)  to  be  given  quan¬ 
tities  along  with  the  plant  characteristics.  A  "reasonable"  performance  criterion  is 
selected  for  minimization,  and  the  results  of  computer-based  optimization  procedures  are 
the  closed-loop  dynamics  and  system  performance  (Refs.  131-1 55) • 

Each  of  these  approaches  has  a  common  basis  in  the  experimental  data  and  mathematical  pilot  models  described 
in  previous  chapters,  although  the  basic  information  about  the  models  may  be  applied  In  quite  different 
forms.  Among  the  cognoscenti  there  are  proponents  of  each  approach,  but  the  authors'  view  is  that  of  the 
middl.  ..an,  i.e.,  all  have  merit  and  tend  to  supplement  each  other  rather  than  compete.  In  fact,  all  three 
procedures  share  in  the  most  significant  single  statement  that  can  be  made  about  the  rise  of  pilot/vehicle 
analysis,  that  it  has  been  established  as  a  fundamental  mode  of  thinking  on  the  part  of  technical  practi¬ 
tioners  in  the  fields  of  aircraft  flying  qualities,  pilot/vehicle  control  system  integration,  display  syn¬ 
thesis,  and  associated  fields. 

Although  the  bibliography  titles  alone  give  a  cross  section  of  applications,  these  may  not  be  as  readily 
appreciated  as  a  more  direct  statement  of  pilot/vehicle  system  problems  which  h  >ve  been  addressed  using 
mathematical  models.  To  illustrate  applications  in  these  ways,  we  have  divided  the  problems  into  three 
categories  for  convenience,  and  have  presented  them  in  outline  form  in  Tables  lh-16  (Ref.  1 38) • 

1.  Flight  Test  Problems  (Table  it) .  The  first  category  is  also  the  smallest  because  it 
lists  problems  initially  encountered  far  down  the  line  in  an  aircraft  development,  i.e., 
in  the  flight  test  phase.  In  all  of  these  cases  the  application  of  the  pilot/vehicle 
analysis  procedures  led  to  a  delineation  of  the  cause  of  the  troubles  and  further  appli¬ 
cation  of  the  procedure  led  to  solutions. 

2.  Design  Problems  (Table  15).  By  far  the  most  widespread  use  of  any  good  predictive  tech¬ 
nique  is  in  the  design  phase.  Here,  specific  instances  are  too  numerous  to  list,  so 
more  general  classifications  are  used  in  Table  13. 
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Simulation  Problems  (Table  16).  In  simulation  experiments  the  theory  is  ideal  for 
program  planning  before  the  experiment,  and  data  interpretation  and  generalisation 
afterwards.  Like  all  good  theories  it  has  the  supreme  attribute  of  focusing  the 
experimental  effort  on  critical  issues.  And  again,  like  all  good  theories,  it  has 
the  characteristic  which  permits  the  fitting  of  experimental  data  into  a  broader, 
more  general,  ccntext.  These  generalizations  not  only  describe  the  use  of  pilot 
mathematical  models  in  simulation  problems  but  also  serve  as  a  suitable  windup  for 
the  entire  discussion. 


TABLE  It 

SOME  PAST  APPLICATIC  :S  OF  PILOT- VEHICLE  SYSTEM  ANALYSES 


SITUATION 

CONTROL  PROBLEM 

CAUSES 

Pilot- Induced  Oscillations 

P»tch(S(ngle  Loop) 

Roll  (Single  Loop) 

Sensitivity; 

Bobwelght/Feel  Spring; 

Loss  of  Pilot  Log , 

Clevotor  Rote  Limiting 

Effect; 

Lateral  Bobwelght 

Weapon  Delivery 

Heoding  Aim  Wonder 
(3  loops) 

Loss  of  RcU  Loop; 

Loterol-  Directional  Multiloop 
Cross  Coupling 

Corner  Landing 

Polh  Con  roMnobllity 
toorres  Rate  of  Sink 
(2  and  3  loops) 

Dynomic  Reversal  in  Path 

Attitude  Control 

PilcMSIngle  Loop) 

Improper  Pilot  /  Stability 
Augments*  Matching 

TABLE  i‘ 


PILOT-VEHICLE-OISPLAY  SYSTEM  ANALYSES: 
DESIGN 


Situation 

ANALYSIS  RESULTS 

0O4iC  A«rf»ame  CM  Pnmo*y 

Control  System 

Predict  multiple, closed-loop  pilot* vehicle 
system  problem  areas  and  assess 
possible  solutions 

P.iot/SiabiMy  Augmentation 
Trodeo.fs 

CoM'dote  stability  augmentation  systems, 
pitot  behoviorond  workload, system 
performance  and  compromises, reliability, 
redundancy,  etc 

SAS  Failure  Effects 

Pilot  oct.ons  ond  resulting  aircraft  excursions 

Competing  Pi'ot  O'tplay  Formots 
Manual  Control 

AFCS  Monitoring 

Information  lequirements,  scan  patterns , 
workload,  assessment  foctors  and  enter, o 

Flight  Director 

Commond  display  lows,  stotus  information 
requirements,  flight  director/pilot/ 
stability  augmentation  tradeoffs 

Carrier  Landing  Aids 

Optimum  FLOtS  control  ond  stabilization 

Cotegor  iei  Dond  Cl  Landing  System 

PfObofcf-ty  of  cpptooch  success,  decision  "stole 
windows",  touchdown  statistics, manuol/ 
automatic  tradeoffs,  guidance  sampling 

Energy-T/im  Monogen.jnt 

Simplified  controts /displays 

TABLE  ]6 


SOME  PAST  APPLICATIONS  OF 
PILOT-VEHICLE-OISPLAY  SYSTEM  ANALYSIS: 
SIMULATION 


SITUATION 

ANALYSIS  RESULTS 

Pre*Exp«rimentol  Analysis 

Predict  critical  areos  and  parameters, 
guidance  for  experimental  design, 
pilot  briefing, questionnaire 

Post-Fxpcnmentol  Analysis 

Interpretation  ond  generalization  of  results 

Competing  Piloting  Techniques 

Pilot  control  procedures, system  performance 
ond  safety  margin  differences. Control 
system  refinements  to  simplify 
piloting  technique 

Motion-Cue  Simulation 

Task -dependent  motion  sensitivity, optimum 
washout  design 
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CHAITER  VIII 

CURRENT  3IATU3  OP  PILOT  MODELER}  AND  3UREHS5R  REQUIREMENTS 


In  the  foregoing  chapters  we  have  emphasized  the  application  of  quasi-linear  models  to  a  large  and 
useful  class  of  pilot  control  situations  while  also  presenting  appropriate  time-dependent  models  ior  tran¬ 
sient  and  mode-switching  behavior. 

As  Chapter  I  indicates,  quasi-linear  models  for  human  controller  behavior  trace  their  history  to  World 
War  II.  They  have  been  questioned  and  criticized  from  their  inception  to  the  present.  let  they  are  still 
with  us  and  continue  to  be  applied  effectively  to  a  large  variety  of  problems.  Our  previous  chapters  pro¬ 
vide  insight  into  why  this  is  the  case.  Consider  the  assumptions  upon  which  applications  of  the  model  rest: 

•  The  forcing  function  is  a  stationary,  random-appearing  process. 

•  The  time-averaged  dynamic  properties  of  the  pilot  are  statistically  stable. 

•  The  controlled  element  characteristics  are  time-invariant. 

•  The  quasi-linear  operator  portion  of  the  model  characterizes  a  usefully  large  portion 
of  the  pilot's  behavior. 

•  Remnant,  when  accounted  for  as  a  pilot-induced  noise,  does  not  "hide"  more  fundamental 
phenomena  which  may  be  important  in  an  overall  pilot/ vehicle  system  context. 

The  first  and  second  assumptions  require  a  minimum  time  duration  of  5-10  sec  for  the  generation  of  mean¬ 
ingful  measurements  and  the  correspondingly  meaningful  application  of  the  model.  This  is  readily  realized 
in  a  very  large  number  of  practically  useful  cases.  The  third  assumption  is  realized  for  the  overwhelming 
majority  of  non-emergency  conditions.  Together  these  assumptions  enable  the  calculation  and  use  of  the 
describing  function  for  the  quasi-linear  model.  It  is  the  fourth  and  fifth  assumptions  which  have  been 
historic  sources  of  criticism. 

Early  critics  voiced  the  uneasy  feeling  that  there  was  a  lack  of  elegance  or  perhaps  analytic  power 
inherent  in  accepting  a  remnant  term.  The  feeling  prevailed  that  imposing  the  Procrustean  bed  of  linearity, 
even  though  situation-specific,  upon  the  pilot  thwarted  the  exposition  of  interesting  and  uniquely  human 
attributes.  A  criticism  that  could  not  be  dealt  with  directly  for  many  years  suggested  that  more  sophisti¬ 
cated  and  crafty  model  making  could  be  expected  to  account  for  significant  portions  of  the  remnant  by  means 
of  nonlinear  operations.  Chapters  IV  and  V  provide  the  basis  for  reducing  this  criticism  to  an  ephemera. 

All  the  available  data  on  the  remnant  make  it  clear  that  for  single-loop  systems  it  is  dominantly  due  to 
random  time  variations  f  the  parameters  in  the  pilot's  describing  function  and  due  to  both  this  and  visual 
scanning  behavior  for  multiloop  systems.  The  majority  of  the  remnant  is  thus  an  irreducible  stochastic 
variability  characteristic  of  the  pilot  which  depends  on  the  four  classes  of  variables  influencing  control 
task  behavior.  The  remnant  contributes  to  system  error  and  can  superpose  small  system  oscillations,  but 
it  is  the  pilot's  transfer  characteristics  which  determine  the  system's  basic  dynamic  stability  and  pro¬ 
vide  an  estimate  of  what  is  usually  the  major  component  of  the  error  signal.  From  our  present  understand¬ 
ing  of  the  components  of  the  quasi-linear  model,  its  long  history  of  applicability  is  not  surprising;  it 
should  work  and  it  does. 

For  most  practical  concerns  related  to  pilot/vehicle  system  analysis,  there  are  no  further  critical 
research  issues  for  this  model.  It  is  in  transition  from  the  domain  of  research  to  that  of  routine  appli¬ 
cations.  There  are,  however,  issues  relevant  to  our  basic  understanding  of  human  control  processes,  in 
general,  and  second-order  effects  in  quasi-linear  control,  which  remain  to  be  addressed: 

•  The  short  time  (less  than  5-10  sec)  behavior  of  the  pilot  cannot  be  described 
adequately  by  constant-coefficient  describing  functions.  Current  indications 
from  remnant  data  are  that  the  pilot's  characteristics  derive  from  a  mathematical 
time-varying  operator  containing  the  terms: 

[KPo  +^K(t)]e-[To+Mt)]s 

where  AK  and  At  are  random  processes.  The  distributions  which  characterize  these 
random  processes  are  not  well  known,  and  their  definition  and  dependence  on  environ¬ 
mental,  pilot-centered,  and  procedural  variables  constitute  areas  for  further  research. 

•  The  inputs  sensed  for  VFR  conditions  are  currently  estimated  on  the  basis  of  con¬ 
trol  needs  (i.e.,  what  feedback  paths  are  necessary  or  desirable  for  the  closed-loop 
system) .  The  actual  quantities  perceived  are  likely  to  be  linear  combinations  of 
these,  with  the  weightings  between  the  independent  inputs  fixed  by  the  geometry  and 
perspective  rather  than  being  independently  adjustable  by  the  pilot.  These  aspects 
of  perception  can  have  profound  effects  on  the  closed-loop  analysis  of  various  maneu¬ 
vers  (e.g.,  approach,  landing,  dive  bombing,  etc.).  There  is  a  need  to  refine  our 
understanding  of  the  actual  real-world  cues  used  under  VFR  conditions  so  that  more 
precise  estimates  of  pilot  dynamics,  scanning  behavior,  and  performance  can  be  made 
for  VFR  conditions  and  for  integrated  displays. 

•  The  dependence  of  the  quasi-linear  feedback  model  in  its  steady-state  form  upon  task 
variables  is  fairly  complete.  For  specific  applications  there  is  a  need  for  codifi¬ 
cation  of  the  influence  of  environmental,  pilot-centered,  and  procedural  (e,  a,  and 
n)  variables  on  the  model. 
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•  Scamsing  behavior  has  been  applied  with  considerable  success  in  some  multiloop 
control  configurations,  but  the  data  base  and  models  are  still  in  a  relatively 
rudimentary  state.  More  complete  data  in  general  is  needed  wherein  cositrol  and 
scanning  behavior  is  measured.  Changes  in  scanning  behavior  and  in  the  induced 
remnant  as  functions  of  e,  it,  and  a  and  the  potential  work  overload  and  desta¬ 
bilizing  influences  need  examination. 

•  Possible  crossfeed  elements  of  the  quasi-linear  model  are  presently  based  more 
on  inference  than  on  directly  applicable  measurements.  The  potential  importance 
of  such  elements  to  a  complete  modeling  of  probably  rare  but  demanding,  highly 
coupled  multiloop  situations,  warrants  the  development  of  a  better  experimental 
base  therefor. 

•  Pilot  opinion  ratings,  despite  their  utility  and  record  of  validation  in  the  past, 
exhibit* an  artistic  quality  in  their  development.  The  physical  continuum  under¬ 
lying  the  ratings  is  multidimensional,  and  the  component  dimensions  have  weight¬ 
ings  which  depend  upon  the  aircraft's  mission.  The  engineering  foundation  for 
this  activity  requires  strengthening. 

•  Although  a  battery  of  extremely  promising  techniques  based  on  critical  task  con¬ 
cepts  is  presently  available  for  assessing  pilot  workload  as  an  excess  control 
capacity,  the  interactions  between  workload  and  pilot  skill  progression  have  not 
been  quantified.  There  exists  anecdotal  material  about  the  minimum  effort  char¬ 
acteristic  of  high  skill,  but  the  conceivable  tradeoffs  between  such  things  as 
effort  extended  in  training  and  potential  subsequent  savings  in  piloting  effort 
are  as  yet  undefined. 

•  Mode  switching  in  response  to  either  abrupt  Yc  changes  or  nonstationary  forcing 
functions  has  been  modeled  on  an  ad  hoc  basis.  Our  understanding  of  the  pilot's 
decision  rules,  switching  criteria,  and  the  mating  of  transient  models  with 
steady-state  models  is  incomplete. 

•  Our  knowledge  of  the  duration  of  intermittent  open-loop  behavior  in  prec  gnitive 
and  near-precognitive  conditions  is  sketchy.  There  are  a  variety  of  research 
questions  in  determining  in  useful  detail  the  parameters  of  subjective  predic¬ 
tability  so  that  tendencies  toward  precognitive  behavior  may  be  either  encouraged 
or  discouraged,  depending  on  the  risks  involved  from  potentially  catastrophic 
transients . 
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